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Abstract 

Six Novel mixed-ligand cyanonitrosyl complexes of chromium (I) of composition, [Cr(NO) 

(CN) 2 (L) 2 (H 2 0)] (where L = N, N-dimethylaniline, N, N-diethylaniline, N-ethyl-N- 
methylaniline, N,N-dimethyl-p-toluidine, N,N-dimethyl-m-toluidine or N,N-dimethyl- 
m-anisidine) have been prepared by the interaction of potassium pentacyanonitrosylchromate 
(I) monohydrate, K 3 [Cr(N0)(CN)5].H 2 0 with the said tertiary amines. The complexes, which 
have been characterized by elemental analyses, magnetic measurements, conductance studies, 
electron spin resonance and infrared spectral studies, contain chromium (I) in a low -spin 
(CrNO) 5 electron configuration. 


Introduction 

A literature survey 1 * 14 of neutral mixed-ligand nitrosyl complexes of monovalent 
chromium reveals few reports on such complexes. Although in recent years there has been a 
great inters! in preparation and characterization of neutral mixed-ligand cyanonitrosyl 
complexes of chromium(I), there is no report on cyanonitrosyl complexes of monovalent 
chromium with tertiary alkylanilines. It was, therefore, thought worth-while to synthesize and 
characterize some neutral mixed-ligand cyanonitrosyl complexes of chromium involving 
{CrNO } 5 electron configuration with tertiary alkylanilines like N, N-dimethylaniline, N, 
N-diethylaniline, N-ethyl-N-methylaniline, N,N-dimethyl-p-toluidine, N,N-dimethyl-m- 
toluidine and NJSI-dimethyl-m-anisidine. 
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Materials and Methods 

N,N-dimethylaniline (Thomas Baker & Company, London), N, N-diethylaniline (The 
British Drug House Ltd., Poole, England), and N-ethyl-N-methylaniline, N,N-dimethyl-p- 
toluidine and N,N-dimethyl-m-toluidine (Aldrich Chemical Company, Inc., Wisconsion, 
U.S.A.) and N,N-dimethyl-m-anisidine (Tokyo Chemical Company Ltd., Japan) were used as 
supplied. 

Preparation of the parent compound : The parent compound, potassium 
pentacyanonitrosyl-chromate(I) monohydrare, K 3 [Cr(N0)(CN) 5 ].H 2 0 was prepared by 
following the procedure of Grifith and Wilkinson 9 . 

Analysis of the constituent elements : Carbon, hydrogen and nitrogen present in the 
investigated complexes were estimated microanalytically. Chromium was estimated by the 
method reported 10 elsewhere. 

Physical methods : The magnetic susceptibility measurements were made at room 
temperature by the Gouy method. The apparatus was calibrated using cobalt mercury 
thiocyanate, Hg[Co(NCS)] 4 . The diamagnetic corrections were computed using Pascal’s 
constant 1546 . LR. spectra (4000-600 cmr 1 ) of the complexes were recorded in Nujol mulls 
supported between NaCl plates on a Beckman IR-20 spectrophotometer. Conductances were 
measured in analytical grade ethanol (EtOH) and dimethylformamide (DMF) using dip type 
cell on Toshniwal conductivity bridge, Japan. E.S.R. spectra of the complexes were recorded 
at room temperature on a Varian E-3 spectrometer using powdered sample at the microwave 
frequency of 9.53 GHz. 

General method of preparation of the complexes : To a filtered aqueous solution of the 
parent compound, potassium pentacyanonitrosylchromate (I) monohydrate (0.01 M, 50 ml), 
an aqueous acetic acid solution 10 ml (1:1) of the corresponding tertiary amine (0.02 M) was 
added with shaking and a coloured solid was precipitated on warming for 15-20 min. over a hot 
plate at 80°C. The resulting mixture was freed from the liberated HCN by passing a current of 
C0 2 through the mixture for a few hours. The precipitate was suction filtered and washed 
several times with dilute acetic acid solution and finally with H 2 0 and dried in vacuo over 
anhydrous CaCl 2 at room temperature to a constant weight. The analytical data are given in 
Table 1. 


Results and Discussion 

All the complexes were air stable coloured solids. They were soluble in DMF, DMSO, 
ethanol, and insoluble in nitrobenzene. The complexes were thermally stable and did not melt 
or decompose up to 280C. They decomposed in dilute acids and alkalis only on heating. All 
the compounds after decomposition in KOH followed by acidifying with acetic acid gave a 
pink colour with few drops of Griess Reagent 15 . This reaction indicates the presence of NO 
group in the complexes. The probable reaction scheme for the Griess reaction is given in 
scheme 1. 

The mixed-ligand cyanonitrosyl complexes (see Table 1 for ligand names) were 
synthesized according to the following equation : 
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K 3 [Cr(N0)(CN) 5 ].H 2 0+ 2L [Cr(N0)(CN) 2 (L) 2 (H 2 0)] 

+3CH 3 COOK + 3HCN(g) + H 2 0 

(where L = NNdman, NNdean, NeNman, NNdmpt, NNdmmt, NNdmma). 

The partial replacement of the cyano groups in the hexa-coordinated complex, 
K 3 [Cr(NO) (CN) 5 ] . H z O by two molecules of ligand arises from the trans effect of the NO 
group. The work of Burgess et al , 17 on the stepwise aquation of the 
pentacyanonitrosylchromate (I) anion, [Cr(NO)(CN) 5 ] 3_ to attain [Cr(N0)(CN) 2 (H 2 0) 3 ] , favours 
the reaction scheme. 


KOH,UO 

[Cr(H0)(C«) 2 ( L ) 2 {H 2 0)] D ec oi-nposition ' NO J + Chromium oxide 




Sulphanilic acid 

(A component of Griess Reagent) 



“(-Naphthyl amine Pink coloured dye 

\ Another component of 
Griess reagent) 


Scheme 1- Probable reaction scheme for the Griess reaction. 
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Table 1 - % yield, analytical, magnetic and E.S.R. data of complexes. 


Compound 

% 

Yield 

%Cr 

Found 

(Calc.) 

%C 

Found 

(Calc.) 

Found 

(Calc.) 

%N 

Found 

(Calc.) 

Meff 

t 

‘ 8 ' 

Decomp 

osition 

Temp. 

(°Q 

[Cr(NO) (0% (NNdman) 2 (H 2 0)] 

50 

12.95 

(13.20) 

54.20 

(54.80) 

5.90 

(6.09) 

17.20 

(17.76) 

1.71 

1.981 

>300 

[Cr(NO) (CN^ (NNdean^O)] 

51 

11.45 

(11.57) 

58.35 

(58.65) 

7.10 

(7.40) 

15.32 

(15.55) 

1.70 

1.980 

>300 

[CiCKO) (CN) 2 (NeNman)2(H 2 0)] 

50 

12.00 

(12.34) 

56.20 

(56.86) 

6.38 

(6.63) 

16.00 

(16.58) 

1.72 

1.983 

>300 

[Cr(NO) (CN) 2 (NNdmpt^HP)] 

49 

11.90 

(12.34) 

56.45 

(56.86) 

6.09 

(6.63) 

16.32 

(16.58) 

1.74 

1.982 

>300 

[Cr(NO) (CN^ (NNdmmt) 2 (H 2 0)] 

49 

12.10 

(12.34) 

55.98 

(56.86) 

5.97 

(6.63) 

15.89 

(16.58) 

1.70 

1.985 

>300 

[Cr(NO) (CN) 2 (NNdmma^CHjO)] 

50 

11.00 

(11-47) 

52.10 

(52.85) 

5.99' 

(6.17) 

14.90 

(15.40) 

1.72 

1.982 

>300 


NNdan = N,N-dimethylaniline; NNdean = N,N-diethylanilme; NeNman = N-ethyl-N-methylaniline; NNdmpt - 
N,N-dimethyl-p-toluidme; NNdmmt = N.N-dimethyl-m -toluidine; NNmma = N,N -dimethyl -m -an is idine. 


Table 2 - Important I.R. spectral bands (cm* 1 ) and their assignments. 


Compound 

[Cr(NO) (CN >2 (NNdman) 2 (H 2 0)] 
[Cr(NO) (CN )2 (NNdean) 2 (H 2 0)] 
[Cr(NO) (CN^ (NeN-man) 2 (H 2 0 )] 
[Cr(NO) (CN >2 (NNdmpt)2(H 2 0)] 

' [Cr(NO) (CN^ (NNdmmt)2(H 2 0)] 
[Cr(NO) (CN) 2 (NNdmma^CHjO)] 


t)(NO) u(C~N) 


1710 (VS) 

2160 (S) 

1705 (VS) 

2160 (S) 

1710 (VS) 

2155 (S) 

1705 (VS) 

2150 (S) 

1705 (VS) 

2160 (S) 

1705 (VS) 

2160 (S) 


\)(C~N) d(OH) 


1380 (S) 

3580 (br) 
3380 (br) 

1385 (S) 

3580 (br) 
3400 (br) 

1380 (S) 

3575 (br) 
3400 (br) 

1385 (S) 

3580 (br) 
3380 (br) 

1380 (S) 

3570 (br) 
3380 (br) 

1375 (S) 

3580 (br) 
3400 (br) 
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Compounds were characterized on the basis of the following results : 

Magnetic and E.SJi. studies : The magnetic and E.S.R. data of the complexes are 
reported in Table 1. The magnetic moments, 1.70 to 1.74 B.M., at room temperature and g 
values, 1.980-1.985 of powdered compounds, which are comparable to the literature 16 - 18 
values, are consistent with a low-spin {CrNO} 5 electron configuration of chromium (I). 

IjR. spectra : The important I.R. spectral bands for the synthesized complexes are reported 
in Table 2. The appearance of a very strong band at 1705 to 1710 cnr 1 and a strong band at 
2150-2160 cm- 1 are assigned to o(NO) + and ti(C=N), respectively, in agreement with the 
result reported 5 elsewhere. The broad bands in the 3570-3580 cm- 1 and 3380-3400 cm- 1 regions 
are due to o(OH) of coordinated water 19 in all the complexes. The absorption bands in 
1590-1600 cm- 1 and 830-850 cm -1 regions are assigned to the coordinated water molecule. A 
comparison of the I.R. spectral bands of the free tertiary anilines and their complexes show 
that the o(C-N) observed at 1350, 1360, 1355, 1355, 1350, .1345 cm- 1 in free NNdman, 
Nndean, NeNman, NNdmpt, NNdmmt and NNdmma, respectively, is shifted to 1380, 1385, 
1380, 1385, 1380, 1375 cm- 1 in their respective complexes. This indicates the bonding of the 
tertiary nitrogen 20 to chromium. 

Conductance measurements : The molar conductances in the range 6.8 to 10.0 ohnr^cm 2 
mole -1 in ethanol and 10.4 to 11.5 ohm -1 cm 2 mole -1 in DMF are in agreement with the non- 
electrolytic nature of these complexes 10 . 


NO 



Fig. 1 - Proposed octahedral structure of [Cr(NO)(CN)2(L) 2 (H 2 0)] where L = NNdman, NNdean, NeNman, 
NNdmpt, NNdmmt or NNdmma. 

The analytical data and physico-chemical studies presented above suggest that the 
complexes may be formulated as [Cr(N0)(CN) 2 (L) 2 (H 2 0)]. Since these complexes show one 
CN stretching band and one NO stretching band, it is reasonable to propose an octahedral 
structure with CN trans to CN, L trans to L, and NO trans to water molecule (Fig. 1). 
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Abstract 

Interaction of Ni(II) with 3-indoleacetate, 3-indolebutyrate and 1-naphthaleneacetate ions was 
studied at DME in aqueous and 50% methanol media at |i = l.OM(KCl) at 25+0.1 °C. The metal 
undergoes reversible diffusion-controlled two electron reduction. The overall stability constants of 
1:1 and 1:2 complexes were determined using DeFord and Hume method. Percentage distribution of 
metal ion in various forms in equilibrium as a function of ligand concentration was calculated. 

Introduction 

In continuation of the polarographic determination of the stability constants of the 
complexes of certain plant auxins and fungicides with trace elements 1 ' 4 we report herein the 
quantitative evaluation of the binding affinity of 3-indoleacetate, 3-indolebutyrate and 
1-naphthaleneacetate ions with Ni(H) in aqueous and 50% methanol media at DME. Percentage 
distribution of the metal in various forms in equilibrium as a function of ligand concentration 
has also been calculated. . 


Materials and Methods 

All the chemicals used were of AR grade. Triple-distilled mercury and double distilled 
water were used. Methanol was purified and distilled before use by the literature method 5 . 

Stock solution (10 2 M) of Ni(II) was prepared in water and standardized 6 . Sodium salts of 
3-indoleacetic, 3-indolebutyric and 1-naphthaleneacetic acids were prepared in aqueous and/or 
50% methanol media and pH of these solutions was maintained with sodium hydroxide. 
Solutions containing metal ions (4xl0~ 4 M) and varying concentrations of the ligand 
(0.00-0.25M) were prepared in water and 50% methanol media at ionic strength 1.0M 
maintained with potassium chloride. 

The capillary characteristics measured in 0.1M potassium chloride (open circuit) and at a 
mercury height of 25.5 cm were #z=1.07 mg s _1 and t = 2.06 s. Purified N 2 gas, presaturated 
with the background solution to be polarographed, was used for deaeration and an inert 
atmosphere was maintained over the solution during electrolysis. Electrolysis was carried out 
using^a thermostated H-cell in conjunction with saturated- calomel electrode. Polarograms of 
4x10 M of Ni(II) were obtained in the presence of different concentrations of the ligand at pH 
10.2 ± 0.1 in aqueous medium and pH 7 ± 0.1 in 50% methanol at 25 ± 0.1°C using EG & 
G Princeton Applied Research Model 174A Polarographic Analyser and Model RE 0089 X-Y 
recorder. While working in aqueous-methanol mixture medium, I.R. compensation was used. 
Gelatin was used as the maxima suppressor. Necessary corrections were made in processing 
the diffusion current data. The resulting polarographic data as a function of ligand 
concentration (C x ) which was calculated from the pK value of the auxin are given in Tables 
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Table 1- E m and F y <[X]) function for Ni (Il)-l-naphthaleneacetate system in 50% methanol at 25±0.rC. 


C* 

flVJ) 

Ejji 

-V(vs. S C E) 

u 

QjA) 

/•ocm) 

FlflX]) 

f 2 ([X]) 

0.00 

1.097 

2.850 

- 

- 

- 

0.02 

1.098 

2.934 

1.11 

5.65 

32.50 

0.04 

1.099 

2.8£4 

1.24 

6.01 

25.25 

0.05 

1.100 

2.839 

1.32 

6.50 

30.00 

0.09 

1.103 

2.835 

1.64 

7.20 

28.00 

0.10 

1.103 

2.800 

1.75 

7.50 

25.00 

0.12 

1.105 

2.785 

1.96 

8.00 

25.00 

0.15 

1.107 

2.750 

2.35 

9.00 

26.06 

0.16 

1.107 

2.700 

2.46 

9.15 

25.93 

0.18 

1.108 

2.645 

2.71 

9.50 

25.00 

0.20 

1.109 

2.600 

3.00 

10 ; 00 

25.06 

0.25 

1.112 

2.550 

3.81 

11.25 

25.00 





Pi= 5.00; 

P 2 = 25.00 


Ni 2+ = 5x10“^!; p = 0.01M (KC1); pH = 7 ±0.1 



Fig. 1- Plot of C x and FjflX]) functions for (a) nickel (U) - 3-indolebutyrate system (- 0 - 0 -), (b) nickel (II)- 
3-indole acetate system (-A-A-), and (c) nickel (13)- 1 -naphthalene acetate system (-n-Li-) in 50% 
methanol medium. 

Results and Discussion 

Hie polarographic reduction of Ni(II) gives a single well defined reversible wave. Current 



Table 2 - E l(2 and Fj {[ X J) function for Ni(H)-3-indolebutyrate system in aqueous and 50% methanol media at 25 ± 0.1 °C. 
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is directly proportional to (h e ff) 1/2 indicating that the reduction is diffusion-controlled. The 
study reveals that the cathodic reduction of Ni(II) in 3-indoleacetate, 3-indolebutyrate and 
1-naphthaleneacetate involves two electron transfer. Decrease in diffusion current of the metal 
ion and shift in its half-wave potential to more negative values as a function of ligand 
concentration suggests complex formation. 

Fj([X ]) functions were calculated at different ligand concentrations 7 and then solved for the 
stability constants by a graphical procedure. Plot of F 0 (\X\), /^([X]) and F 2 ([X]) against C x is 
a smooth curve, a straight line whose slope gives the value of (3 t and a line parallel to the 
abcissa respectively confirming the formation of 1:1 and 1:2 complexes in all the cases, 
depicted for nickel(II)-3-indolebutyrate system in 50% methanol medium as an example, in 
Fig. 1. 

The stability constants of 1:1 and 1:2 complexes of Ni(II) with 3-indoleacetate, 
3-indolebutyrate and l-naphthaleneacetatc in aqueousmethanol medium at pH 7 ± 0.1 are 
Pi = 7.50, p 2 = 40.46; Pi = 5.95, p 2 = 34.96 and p x = 5.00, p 2 = 25.00 respectively. Thus the 
complexes of 3-indoleacetate are stronger than the corresponding complexes of 
3-indolebutyrate which in turn are stronger than the 1-naphthaleneacetate complexes. The low 
values of stability constant of 3-indolebutyrate complexes may be due to the presence of more 
methylene groups in the chain having carboxylic group. The basicity of the ligand and the 
steric effect influence the stability of the complex. It is noted that the complexes in water are 
stronger than those in methanolic medium. (The values for 3-indoleacetate and 
3-indolebutyrate complexes in aqueous medium at pH 10.2 ± 0.1 are p x = 8.90, p 2 = 42.48, 
and Pj = 7.00 and p 2 = 36.25 respectively). The stability increases with the increase in pH of 
the medium. 

Percentage distribution of the metal ion present in various forms in equilibrium as a 
function of the ligand concentration was also calculated and the results are presented in 
Fig. (2&3). 
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Abstract 

The compounds, Pr 2 C 12 H 12 0 12 .6H 2 0 and Dy 2 C 12 H 12 0 12 .8H 2 0 have been prepared and characterised 
by elemental analysis and I.R. spectroscopy. TG and DTA studies show that these lose four and 
six water molecules respectively in the first step, yielding M 2 C 12 H 12 0 l2 .2H 2 0 (M = Pr or Dy) and 
fonn the anhydrous species, M 2 C 12 H 12 0 12 in the second step. There is an indication of formation 
of MO(C0 3 ) 1/2 in the case of Pr compound. M 2 0 3 appears as the final stable product in each case. 
Kinetic parameters like apparent order, energy and entropy' of activation for the dehydration 
process have been arrived at by analysing the data using Piloyan-Novikova, Coats -Redfem and 
Horowitz-Metzger equations, a vs. T curves show that there is no indication of branching or 
surface nucleation preceding the dehydration process which in both the compounds proceeds with 
phase boundary mechanism. 


Introduction 

Tartarates and succinates of rare earths, including those of praseodymium and dysprosium, 
have assumed considerable importance in view of being good starting materials for oxide 
catalysts and also due to their use in several electronic devices. Besides being good materials 
for studying decomposition mechanisms their decomposition products, including 
intermediates, could be of interest in the search for materials of high electrical conductivity. 
Hence, in continuation to our work on rare earth tartarates 1 ' 3 , it was thought worthwhile to 
study the thermal decomposition of Pr(EI) and Dy(III) succinates. 

Materials and Methods 

The compounds Pr 2 C 12 H 12 0 12 .6H 2 0 and Dy2C 12 H 12 0 12 .8H 2 0 were prepared by mixing 
aqueous solutions containing equimolar (0.01 mole) amounts of respective rare earth metal 
chloride (Indian rare earths, 99.9% pure) and succinic acid (LOBA, 99% pure) containing some 
pyridine (Ranbaxy, A.R. 99.5% pure) at room temperature. Presence of pyridine seems to 
control the pH of the solution to the required extent (~4-5). The fine crystals of the 
compounds so obtained were separated by filtration, washed with distilled water and air dried 
without pressing. The percentage of the metal was estimated gravimetrically by precipitation 
as M(OH) 3 and weighing as M 2 0 3 (M = Pr or Dy). The percentage of carbon and hydrogen 
was found by micro-analytical method. 

IR-435 spectrophotometer Shimadzu, (Japan) was used for I.R. spectral studies, using 
KBr pellets in the range 4000400 cm 4 . 

Derivatographic studies of both the compounds were carried out on Paulik-Paulik-Erdey 
MOM derivatograph (Hungary, Budapest) using a-alumina as the reference material up to a 
temperature of 1173K. The thermal runs were recorded at constant heating rate of 10° min- 1 in 
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static air atmosphere using 200 mg sensitivity. As a further check, TG and DTA curves were 
also recorded on Stanton Redcroft 780- series instrument (England) for one of the compounds, 
Pr(III), using alumina as inert material at 5MW sensitivity and 10° min -1 heating rate. The 
results from both the instruments were found wholly compatible. 

Reflectance spectra of Pr 2 C 12 Hi 2 0 12 .6H 2 0 was recorded on Beckmann-DK-2A reflectance 
spectrophotometer (MgO as reference material, speed = 36 nm/min; scale = 10 nm/cm and 
sensitivity = 50). 

ECIL-TDC-316 computer (India) was used to evaluate the various kinetic parameters viz. 
energy of activation, pre-exponential factor, entropy of activation, apparent order of the 
reaction and the mechanism involved in the dehydration process in conformation to integral 
and approximation methods using the Piloyan-Novikova 4 , Coats-Redfern s and Horowitz- 
Metzger 6 equations (non-isothermal). 

Results and Discussions 

On elemental analysis (Table 1) the compounds were found to have the composition 
Pr2C 12 HiA2.6H 2 0 and Dy 2 Ci 2 Hj 2 0j 2 .8H 2 0. 


Table 1 - Values obtained by elemental analysis. 


SL 

no. 

Compound 

Metal % 

C% 


H % 



Found 

Calc. 

Found 

Calc. 

Found 

Calc. 

1 . 

Pr 2 C 12 H 12 0 12 .6H 2 0 

38.07 

38.19 

19.40 

19.52 

3.00 

3.25 

2. 

Dy 2 C 12 H 1 A 2 .8H 2 0 

39.93 

39.78 

17.50 

17.62 

3.50 , 

3.43 


The I.R. spectra of the compounds show a strong and broad band at 3400 cm -1 due to OH 
stretching of water. The broadness of the band is expected due to the intra-molecular hydrogen 
bonding in the compounds. The strong, band -1575 cm -1 can be assigned to t> asym (OCO) and 
the bands -1425 and 1320 cm -1 to u sym (OCO) of the coordinated carboxylate groups 7 . 

Thermal decomposition : The thermal decomposition pattern of both the compounds is 
very similar (Fig. 1-3). The first step of decomposition involves partial dehydration leading to 
the formation of the same’ product viz. dihydrate in both cases. In the case of the Pr 
compound, the first decomposition step starts at 315K and continues up to 433K, the weight 
loss observed (Derivatograph 9.50%; Stanton 9.40%) corresponds to removal of four water 
molecules (calc, loss 9.70%). In the case of Dycompound this decomposition step starts at 
310K and ends at 393K; during this step six molecules of water are removed (observed weight 
loss 13.30% calc, loss 13.20%). This step can be represented as 

M 2 (C 12 H 12 0 12 ) 3 x H 2 0 -> M 2 (C 12 H 12 0 12 ) 3 2H 2 0 + (x - 2)H 2 0 

In the second step of decomposition both the compounds loose the remaining two H 2 0 
molecules leading to the formation of the anhydrous salts. In the case of the Pr compound this 
decomposition step starts at 473K and continues up to 520K showing an observed weight loss 
of 15.0% (Derivatograph), 14.8% (Stanton) (calc. wt. loss 14.6%). In the case of the Dy 
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compound this step starts at 400K and ends at 516K with an observed weight loss of 18.2% 
(calc. 17.6%) 

M 2 (C 12 H 12 0 12 ) 3 2H 2 0 M 2 (C 12 H 12 0 12 ) 3 + 2H 2 0 

The anhydrous salts show remarkable stability over a range of temperature as evidenced by the 
corresponding smooth plateaus in the TG curves. The anhydrous Pr salt remains stable 
between 520 and 698K while the Dy shows stability in a similar temperature range viz. 516 
to 693K. 

In the third step anhydrous praseodymium succinate starts decomposing at 698K. The 
decomposition continues up to 800K giving basic carbonate 8 9 PrO(C0 3 )i ^ as indicated by the 
observed weight loss (Derivatograph 50.0%; Stanton 50.5% calc. 50.6%) up to this stage. 

M 2 (C 12 H 12 0 12 ) 3 — t93/2 ° 2 ■ > 2M0(C0 3 )i + 35C0 2 + 18H 2 0 

This intermediate product remains stable over a relatively small range of temperature (800- 
820K) and decomposes between (820 and 958K) leading to the formation of the oxide M 2 0 3 as 
the final product as shown by the total observed weight loss 55.8% (Derivatograph); 55.4%' 
(Stanton); (55.8% calc.) 



Fig. 3 - TG-DTA curves of Dy2C 12 H 12 0 12 .8H 2 0. 


0 


In the case of dysprosium succinate the basic carbonate intermediate does not appear to be 
stable. The decomposition steps involving the formation of the basic carbonate and its 
subsequent decomposition overlap and the final product formed is the oxide. 
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M 2 {C 12 H 12 0 12 ) 3 — +57/2 °- 2 - > M 2 0 3 + 36C0 2 + 18H 2 0 

This overall step of decomposition starts at 693K and ends at 985K beyond which the oxide 
M 2 0 3 remains stable. Total weight loss up to the formation of M 2 0 3 is observed fo be 54.2% 
(calc. 54.3%). 

DTA curve of praseodymium succinate shows two endotherms with peaks at 368 and 
513K and two exotherms with peaks at 763 and 823K. In the case of dysprosium compound, 
DTA shows two endotherms with peaks at 373 and 493K. The exotherm peak is seen at 
763K. Thus the steps of decomposition indicated by TG correspond to those indicated' by 
DTA. The endotherms correspond to decomposition processes while the exotherms represent 
decomposition accompanied by oxidative processes. 

The reflectance spectra of the original compound (praseodymium succinate) and the 
species got after heating the compound at different temperatures are shown in Fig. 4. Besides 
room temperature (303K), the reflectance spectra were taken for sample heated to 403, 463 and 
523K. A decrease in the intensity of bands is observed near 590, 485, 370 and 445 nm as 
dehydration proceeds with the rise in temperature. 

The I.R. spectrum of the dehydrated product shows all the characteristic bands of the 
starting material except the band at 3400 cm -1 which vanishes completely at 525K (for Dy 
compound) suggesting the formation of the anhydrous salt 


Kinetics of decomposition : The fractional weight loss (a) at different temperatures was 
calculated from TG traces and a vs. T curves plotted for the first decomposition step in the 
case of both the compounds (Fig. 5). The curves do not begin with any apparent induction 
period and thus suggest that die initial stages involved very little or no physical desorption 
and also indicate that no surface nucleation or branching occurs before the decomposition 
starts 10 . The beginning of a-T curves with acceleratory period, followed by a region df 
maximum rate which virtually continues up to the end of process (a=l), is suggestive of there 
being no significant retention or decay period. 


The kinetics of first decomposition step for both the compounds have been studied by 
non-isothermal method using the following equations : 

(i) Coats and Redfem equation 


(a) 



= log 


-ln(l-q) 


E 

2.303RT 



...forn = 1 


( 1 ) 


(b) log = log - 


i-d-q) 


1 -* 


J-n 


2303RT 


+ log 



1 2RTS 
E ) 


..forn^l 


(ii) Piloyan and Novikova equation : 

g( a) a ZR E 

log X- = log— y = log- 

jT 2.303 RT 


..(a = 0.05-.5) 


( 2 ) 


(iii) Horowitz and Metzger equation 






LOG 
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T (K) ► 


Fig. 5 - a-T curves : a, P r 2 ^i 2^ 1 2^ 1 2 - 2 O , b, 




Fig. 7 r Piloyan-Novikova Plots : a, Pr 2 ^i 2 ^i 2 ^i 2 *^ 2 ^* b* 
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£0 

(a) logg(ct) - log[~ln(l-ce)] = - 3Q 5 R T n ? f ° r " = 1 


(6) logg(a) = log 


1 — n 


l-n 


EQ 


l.mRTm 


■j for;) * 1 


( 3 ) 


where T = temperature; R - molar gas constant; £ = energy of activation; Tm = peak 


temperature; 0 = T-Trrr, Z = pre-exponential factor; p = rate of heating and #(<*) ™ J 


da 

f( a)’ 

where /(a) is the appropriate function of a depending upon the mechanism involved). #(a) 
values corresponding to different models were got computerised and the values of g(a)/P v.v. 
1/T plotted for eqn. (1) and (2). The linear fits for different models were investigated. The best 
fit linear plot (Fig. 6) with minimum least square error was selected while applying Coats- 
Redfem equation. The values of slope, intercept, E and Z, were calculated from this very plot 
as also in the case of Piloyan-Novikova plot (Fig. 7) and the mechanism was decided. The 
value of entropy of activation AS* was computed from the equation 
kTm . , 

z = — exp(AS*/fl) (4) 

where K is the Boltzmann constant and h is the Planck constant. * represents the observed 
values and + the values calculated for ideal linear fits respectively. Plots of g(a) v.v. T were 
constructed while applying eqn. (3) and the best fit linear plot (Fig. 8) having minimum least 
square error was selected to find g(a) and hence the mechanism. E was calculated fropi the 
slope of this plot and Z wasidetermined from the equation 
E 


Z = 


RTm 


$exp(E/RTm) 


The value of entropy of activation, AS*, was calculated from eqn. (4). 
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Table 2 presents the computed kinetic results. It is found that the best linear fit for both 


the compounds is obtained for « - 2 or ■?( a ) ~ ^ 
Horowitz-Metzger method. 


1 - (1 - a )2 


by both Coats-Redfem and 


Table 2 - Values of the computed kinetic results. 


SI. 

no. 

Compound 


Equation 

E 

(KJ mol -1 ) 

z 

(S- 1 ) 


s* 

(JK-’mol- 1 ) 



0) 

Coats-Redfem 

17.26 

0.43 x 

10 1 

-253.75 

1. 


(ii) 

Piloyan-Novikova 

36.95 

1.0 x 

10 3 

-188.91 



(iii) 

Horowitz-Metzger 

22.59 

0.54 x 

10 1 

-232.66 



(i) 

Coats-Redfem 

41.27 

4.4 x 

10 3 

-177.23 

2. 

^2^12^12^2*^20 

(ii) 

Piloyan-Novikova 

53.06' 

3.9 x 

10 5 

-139.82 



(iii) 

Horowitz-Metzger 

55.36 

2.7 x 

10 5 

-143.05 


Phase boundary mechanism involves the progress of the decomposition reaction along the 
reactant product interface and the reaction rate at any stage is proportional to the area of the 
interface at that stage. If a solid follows “contracting cylinder” mechanism of decomposition, 
its rate of decomposition is proportional to surface area 

da/dT « 2rcrZ (for a cylinder) 


M 

For a degree of decomposition, a, the volume of the sample is — (1 - a), ( M = mol. 

d - 

wt, d = density), and therefore 

— (1 - a) = nr 2 1 (for a cylinder) 


and r - — -(l-a)i; therefore 

did 


da „ , . f M 

— -«2nr/«2jt — 
dT Idnl 

= K{l-a$ = Kf{a) 


(l-a) . 


and * (p,)= |^)’ 2 [‘- 2(1 - Il)i ] 

Since the use of this value of g(a) in Coats-Redfem and Horowitz-Metzger equations gives the 
best linear fit for both the compounds, phase boundary decomposition model with “contracting 
cylinder” mechanism is suggested for the first step of their decomposition. 
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The values of kinetic parameters obtained by the application of different equations are 
reasonable and in good agreement. The values of AS* suggest that the “transition state” 
involved in the decomposition process has a more ordered structure. 
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Abstract 

A study of the distribution of electric charges of selected soils from West Bengal was made by 
means of potentiometric titration in the presence of varying concentration of indifferent 
electrolyte. The electrochemical behaviour of these soils was found to be similar to that of the 
constant potential systems in which the surface potential is determined by H + and OH“ ions in 
the equilibrium solution, hence charge distribution varied substantially with pH and electrolyte 
concentration. The pH at the zero point of charge (ZPC) (varying from 2.40 to 2.90) occurred 
below the natural pH value of the soil indicating that the soils had a net negative charge under 
natural condition. ^(ZPT-ZPC) values indicated that soils had almost same amounts of permanent 
negative charge. 

Introduction 

Much of the work in soil has been conducted in temperate regions where, in general, soils 
are characterized by surfaces bearing a constant surface charge, whereas soils of the humid 
tropics are characterized by surfaces having a constant surface potential 1 . In such systems, the 
surface potential remains constant and the surface charge varies with changes in electrolyte 
concentration and in the activity of the potential determining H + and OH~ ions 2 . 

The most important parameter to describe a reversible surface charge system is its zero 
point of charge (ZPC). The ZPC is defined as the pH at which the surface charge of the 
system reduces to zero 3 . On the acid side of the ZPC, the system is positively charged and on 
its alkaline side it is negatively charged. Information on the electrochemical properties in Tarai 
and Teesta alluvial regions of West Bengal is rather limited. The aim of this study is to 
examine the ZPC and change of charges in profiles of two soil series in above regions of West 
Bengal due to pH and salt concentrations. 

Materials and Methods 

Nine soil samplesffrom different horizons of two typical pedons located in Tarai and 
Teesta alluvial regions were collected based on the information of soil survey conducted by All 
India Soil and Land Use Survey, Calcutta Centre. General characteristics of the soils are 
reported in Table 1. The soils were air dried under shade, ground and passed through 2 mm 
sieve. The relevant selected properties of the soils are reported in Table 2. Potentiometric 
titrations for determination of the ZPC were done in 0.001, 0.1 and IN sodium chloride 
solution following the method of Van Raij and Peech 1 . 



174 


A. K. DOLUI and S.DEY 


Table 1 - General characteristics of soils. 


Soil series Latitude 

Logitude 

District 

Soil Height 

taxonomy above 

mean sea 
level 
(m) 

Annual Annual 

precipita- mean 

lion temperature 

(nun) (°Q 

Anandapur 26°18'N 

89°41'E 

Cooch 

Typic 


38 

2746,0 

25.00 





Behar 

Udifluvent 




Dalingkote 26*417^ 

88°18'E 

Darjeeling 

Typic 


131 

3483.6 

24.45 






Dystropept ' 




Table 2 - 

Selected properties of soils. 







Horizon 

pH 


specific Organic Cation 

Fe Kd 

A1 KC1 

Particle size 

Texural* 

(cm) 



conduc- carbon exchan- 

(cmol 

(cmol 

distribution (%) 

class 


Bp 

IN 

tance 

(%) ge capa- 

kg’ 1 ) 

kg^ 1 ) 






KQ 

(dStn- 1 ) 

city (c 



Sand 

Silt Clay 






molkg~ 1 ) 









Typic Udifluvent - Anandapur series (Cooch Behar) 



043 

6.00 

4.80 

0.04 

0,62 8.1 

0.020 

0.121 

38.3 

51.7 10.0 

SiL 

13-28 

6.80 

4.90 

0.05 

0.27 5.5 

0.018 

0.111 

45.5 

45.5 9.0 

L 

28-57 

7.75 

4.90 

0.04 

0.42 8.5 

0.019 

0.097 

36.1 

43.0 20.9 

L 

57-79 

7.15 

4.80- 

. 0.05 

0.20 7.3 

0.019 

0.037 

51 . 1 

36.0 12.9 

L 

_ 79450+ 

7.60 

5.00 

0.04 

0.17 6.4 

0.010 

0.038 

68.5 

23.0 8.5 

SL 




Typic Dystropept - Dalingkote series (Darjeeling) 



041 

5.15 

3.90 

0.03 

1.88 8.3 

0.072 

0.360 

63.0 

21.5 15.5 

SL 

11-44 

5.30 

4.00 

0.02 

1.69 9.2 

0.065 

0.356 

51.0 

27.3 21.7 

SCL 

44-77 

5.30 

4.40 

0.02 

1.31 8.1 

0.066 

0.269 

61.0 

12.7 26.3 

SCL 

77430+ 

5.25 

4.20 

0.01 

1.06 6.4 

0.054 

0.295 

64.5 

14.7 20.8 

SCL 


*SCL - Sandy clay loam, SL - Sandy loam, L-Loam, SiL - Silt loam. 


Results and Discussion 

Net electric charge by potentiometric titration : Potentiometric titration curves at different 
ionic strengths for two samples one each belonging to Anandapur series and Dalingkote series 
are shown in Fig. 1. The other samples, in each of these series, behaved in a similar manner. 
The charge-pH curves measured in the various ionic strengths intersected one another at single 
point indicating that these soils belong to the constant potential type of colloids with H + and 
OH~ as the potential determining ions. The ZPC values ranged from 2.40 to 2.90 (Table 3). In 
all cases, the ZPC values occurred below the natural pH values of the soils indicating that soil 
colloids bear net negative charges under natural conditions. This is also in agreement with the 
negative delta pH values (ApH = pH KC1 - pH H:l o) of the soils (Table 3). The ZPC of the 

charge-pH curves is on the acid side of the zero point of titration (ZPT), whereby indicating 
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that soils contain appreciable amounts of permanent negative charge 1 . 

Table 3 - ApH and ZPC values for the two soils. 


Horizon 

(cm) 

ApH=pH Ka -pHH 2 o ZPC Pikers 

dp(ZPT-ZPC) 

(cmolkg~ l ) 

Surface electric potential at 
soil 

pH=59(ZPC-pH)mv at 25°C 



Typic Udifluvent - 

Anandapur series (Cooch Behar) 

0-13 

- 1.20 

2.40 

2.40 

4.94 

- 212.4 

13-28 

- 1.90 

2.45 

2.45 

4.95 

- 256.6 

28-57 

- 1.85 

2.50 

2.40 

4.96 

- 250.7 

59-79 

- 2.35 

. 2.55 

2.25 

4.97 

- 211.4 

79-150+ 

- 2.60 

2.60 

2.40 

4.98 

- 295.0 



Typic Dystropept 

- Dalingkote series (Darjeeling) 


0-11 

- 1.25 

2.45 

1.45 

4.95 

- 159.3 

11-44 

- 1.30 

2.60 

1.40 

4.96 

- 159.3 

44-77 

- 0.90 

2.80 

1.60 

4.97 

- 147.5 

77-130+ 

- 1.05 

2.90 

1.30 

4.98 

- 138.6 


MQi 




Fig. 1 - The net electric charge vs. pH curves of soil as determined by potentiometric titration. 

The values of ZPC (Table 3) obtained from the potentiometric titration curves show 
strikingly the effect of the higher organic matter contei t of the surface horizon in lowering the 
ZPC, when compared with the ZPC of the corresponding subsurface horizons. The ZPC of the 
soils is different for different soils due to the difference in organic matter as well as 
sesquioxides content in soils. 

The distribution of net electric charge on soil colloids (Fig. 1) indicates that the charge is 
influenced by the pH as well as the concentration of the.NaCl. The common point -of 
intersection of potentiometric titration curves of soil samples obtained in the presence of 
varying cocentration of indifferent electrolyte (NaCl) is not affected by the addition of salt, was 
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called the ‘point of neutrality’ by Mattson and Wiklander 4 who interpreted this unique point as 
the pH of the soil at which the exchange acidity due to displacement of H + ions was just 
balanced by the displacement of OH ions by the added salt. Except at ZPC, the net electric 
charge of the soil colloids has been increased with the increasing concentrations of indifferent 
electrolyte (NaCl). Above ZPC, the negative charge of the soils increases with increasing pH 
whereas the positive charge of the soils increases with decreasing pH. 

The [pH m - ZPC] values (Table 3) of Anandapur series vary from 2.25 to 2.45 and 1.30 
to 1.60 in Dalingkote series. It is inferred that Dalingkate soil series is more developed than 
Anandapur series. Because with increasing pedogenic development, the ZPC approaches the 
natural pH of the soil 5 - 6 . The Oi(ZPT - ZPC) values indicate the presence of 4.94 to 4.98 
cmolkg -1 of pH-independent negative charge in the soils. 0) values also reveal that all soil 
samples bear almost same amount of negative pH-independent permanent charge. Surface 
potential of the soils at natural pH varies from -138.6 to -256.6 mV at 25°C. 

From the present study it may be concluded that the charge distributed on the soil colloid 
increases with the increase in the concentration of indifferent electrolyte. Hence, both cation 
exchange capacity and anion exchange capacity measurement of these soil systems would be 
necessary to determine both positive and negative electrical charges of the soil in a dilute 
electrolyte with concentration similar to that encountered in the soil solution. 
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Abstract 

Two fractional integration operators associated with a multivariable /-function have been 
defined. Three theorems are established. Some formal properties of these operators are also 
given. 


Introduction 


In the present paper two fractional integration operators with new representations 
introduced, whose kernels are associated with the multivariable /-function due to Prasad 1 , 
important theorems of these operators are established. Results given earlier by 
authors 2-13 are the special cases of the results proved here. 


The multivariable /-function introduced by Prasad 1 will be defined and represented in 
following manner : 


% •-•*>■] = 7 , . V (r) (,)\ 

p v q 2 \...:p r ,q r ip ,q ,q J 


z l>— * z r 


b 


J--' 


(2nw) 
where w = V-l 


h U 


JO 


.(0 
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a^,Py\ajj,y ,pjj.y (j = l,...,r) ( k = l,...,r) are positive numbers, a‘\ 


b^ (i=l,...,r),a kJ , 


b kj (k=2,...,r) are complex numbers and here m®, n®,p®, <7® (i=l,...,r), /z*, p k , q k ( k=2,...,r ) 
are non-negative integers where 0 < /n® < <g®, 0 < «® < p®, <7* > 0, 0 < n k < p k . Here (i) 
denotes the numbers of dashes. The contours L ; in the complex 5,-plane is of the Mellin- 
Bames type which runs from -woo to +w°° with identations, if necessary, to ensure that all the 
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For further details and asymptotic expansion of the /-function one can refer to the original 
paper Prasad . I function includes, as special cases, the functions studied by Fox , 
Agarwal andBraaksma . 

In what lollows, the multivariable /-function defined by (1) will be represented in the 
contra ^ed notation. 


P2-P2-— :Pr.7r:(A.<7%-- : (y r V r ^ ^ 1 


Z r] 


or simply byI(z lt ...,z T ). 
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Integral Operators 

equations*" 110 ^ 06 ^ fractional “teg 1 " 31 ! 011 operators by means of the following integral 


{Krf^Alfix)^ 11 S'*- 1 j - £) a f(t).l(K 1 lf„..,K r U)dt 
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and mi, m are positive numbers. 


The condition of the validity of these operators are as follows : 
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The Mellin transfonn off (x) will be denoted by M \f (x)] or by F(s). We write s = p' 1 + it, 
where p and t are real. If p > 1,/ (x) « L p { 0,«>), then 

P = 1 . M[f{x)] = F{s) = °\f{x)x S ~ l dx ( 6 ) 

0 


p > 1, M[f{x)] = F(s) = L i. m. J f(x)x s 1 dx 


where l.i.m. denotes the usual limit in the mean of L p -spaces. 

Theorems 

We now enunciate the three main theorems of this paper. 

Theorem 1: If /(*) € L p (0,oo), 1< p< 2[or f(x) e M p { 0,~) and p > 2], 

1 . r b {n) . . r 6 (0 
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den °f/ theclass c ofaI1 Actions M ofL p (0,~) with p >2, which are inverse 
Melirn transforms of functions of L q (- 00,00). 

Theorem 2: Uf(x)e L p ( 0, °°), 1 < p < [or /( x) e M p { 0, ~), p > 2], |arg < - \ iK 
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Proof of Theorem 1 : From (4), it follows that 
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On changing the order of integration which is permissible 18 under the conditions stated 
with the theorem, we obtain 

0 t 

The theorem now follows on evaluating the ^-integral by means of the formula 
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which follows from the Beta function formula and includes Goyal’s formula 19 . 

Theorem 2 can. be established in the same way. Theorem 3 follows immediately on 
interpreting it with the help of (4) and (5). 
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Formal Properties of the Operators 

We give some formal properties of the operators which follow as a consequence of the 
definitions (4) and (5). 

a> x~\K r ) n ’ a A^ lx [f(i ! x)] = 

(u) x -1 (^) 5ia 4x[/(i / *)] = (* r f ,a 4[*~ l /(*)] 

(iii) x^(K r )^ a Al[f(x)] = (ir r f-^J[*'/(x)] 

(iv) *P(z: r ) 5 ' a 4[/W] = (^) 6 ~ p '^p/«] 

(v) If {K r f' a Al[f{x)\ = g(x) 

then = g{cx} 

(vi) andIf(/: r ) 5 ' a 4[/W] = 4»W 
then (AT r ) 8,a B|[/{cx}] = (|){cx} 
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Abstract 

The two-phase homogeneous conducting system has been reduced to a corresponding single 
phase system with the help of two suitably defined quantities for the mixture-uniform density 
p m and permeability coefficient \i m . The equations governing the flow have been integrated 
across the shock surface to give shock relations. The presence of magnetic field affects the- 
propagation of weak discontinuity in such a way that, in place of one, three modes of 
propagation of surface of wfcak discontinuity have been obtained. It has been observed that 
Alfven velocity is independent of velocity of propagation of sound in the mixture while the 
slow and fast wave velocities depend upon sonic velocity. 

Introduction 

Two phase hydrodynamical systems have not been investigated as deeply as single phase 
hydrodynamical system for the want of appropriate mathematical formulation. As a matter of 
fact, single phase system is only an ideal mathematical formulation far away from the natural 
reality. Most of the natural phenomena are represented by two phase/multi-phase system only. 
The study of two phase system is equally important for subsonic and supersonic flows. Some 
of the initial work in this field has been done by Carrier 1 and Campbell and Pitcher 1 . But they 
have started with jump conditions across the shock without going' into the details of 
formulation governing the flow. The work was extended by Mitira 3 , Moore 4 , Rudinger 3 and 
Wizngaarden 6 etc. Upadhyaya and Singh 7 have recently formulated the problem in somewhat 
more general sense and have observed that the results obtained by Campbell and Pitcher 2 for 
homogeneous gas liquid mixture are also true for conducting, non conducting and viscous 
conducting fluids. 

Our two-phase conducting system is a liquid-gas mixture in the presence of magnetic 
field. The former (liquid) is incompressible where as latter (gas) is compressible. In the 
mixture the gas is supposed to be in the form of bubbles. An incompressible fluid mixed with 
gas bubbles behaves like compressible one and thus its fundamental characteristic is changed. 
Formulation of a two phase system is a difficult task. Presence of magnetic field makes it 
more complicated. In Campbell and Pitcher’s formulation 2 relative velocity between the 
phases was not accounted for. Wallis 8 made some improvement in the approach of Campbell 
and Pitcher. Some workers 911 have made significant contribution by improving the 
formulation available earlier. Miura’s 3 formulation can be treated as an improvement on the 
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earlier formulations but the formulation proposed by Upadhyaya and Singh 7 is more 
appropriate and has been used here. 

Formulation 

The bubbles affect and are affected by the motion of liquid. The effect of gas bubbles is 
proportional to the space that they occupy. Let us assume that A be the space occupied by the 
gas in liquid-gas mixture of unit volume; the space occupied by the liquid will be (1-Y). The 
liquid and gas are supposed to be moving with the same velocity. Thus, there is no interaction 
between the phases 2 . Since the magnetic field has no impact on mass, concentrating on liquid 
and gas portions separately and applying the law of conservation of mass, equations of 
continuity for the individual phases can be given as 

p *^^+p*{<i-* )v ,}.,= o cw 

(lb) 

where v t - is the common velocity of liquid and gas and comma followed by Latin index T 
denotes partial differentiation with respect to x it p and p* stand for gas and liquid densities 
respectively. Every term of field equation contains magnetic field intensity (//). The field 
equations can be directly written as 

Hij = 0 (2a) 

+ VjHij - HjV it j + HjVjj = 0 (2b) 

Unlike mass, the magnetic field has impact on moving fluid which is directly related to its 
velocity. Applying laws of conservation of momentum and energy, we get the equations of 
motion and energy for liquid and gas separately as 

(1 -*)p * + (1 - Y)p * VjV itj + {(1 - Y)^ . 

=0 ( 3a ) 


*p|4 Xp +VjVi j + (Xp) . + -±~ p g x\ = 0 


-1 y— s r*-* U |(i-x)p*Jlv 2 + c-r + 4 


: (1 -X)H 2 Vjj - ( 1 - X)pv iti + - X) Vi (H%j - HiHj)} = 0 
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where pe and p^ stand for magnetic permeability of liquid and gas respectively. c v , c* and T are 
respectively the specific heat of gas at constant volume, specific heat of liquid and temperature 
of the media. 

The two media have got different permeabilities and densities. Therefore, we have obtained 
separate equations of continuity, motion and energy for respective phases. To combine the 
equations of continuity, motion and energy, we define two quantities p m and p m respectively 
for density and magnetic permeability of die mixture as 

P/n = Xp + (l-X)p* (5) 

p m = (l-X)p e + Xpg (6) 


We combine the set of eqn. (1), (3) and (4) by using (5) and (6) to give 


dt 


L + v /Pm,i + P/n v i,;-0 


3Vj 1 \ \ H ‘ s uu 

P + PmVjViJ + Pj + — jPml — Sy - HiHj 


= 0 


(7) 

( 8 ) 



(9) 


where c„ 


mc* + c v 
1 + m 


( 10 ) 


c mv represents specific heat of the mixture and m is mass ratio of the liquid to gas such that 


m 


(1-*)P* 

Xp 


and can be treated as constant without affecting the general results 2 . 


( 11 ) 
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Shock Relations 

In order to obtain the shock relations, we intergate the eqn. (2), (7), (8) and (9). For this 
purpose, let us consider a cylindrical fluid volume V , with base £[ and I ^ and curved surface 

S containing the portion X of the shock surface in it. Obviously, S is formed of the stream 
lines of the particles crossing the outer boundary of X! at a particular time and reaching the 

corresponding boundary of Xj. The cylindrical volume of length / is assumed to be stationary 
with tile shock surface. The integration of the representative equation of the set of equations 
governing the flow, i.e., 

fU /- £1 -0 ( 12 ) 

dt ' 


over the volume V , under the conditions / , V , and S tending to zero and Xi and X* tending 
to X yields 12 

JJ[F-n ; ]dX = 0 (13) 

£ 


where [ ] denotes the difference of the quantity on the two surfaces Xj and X^ of the volume. 
The eqn. (13) holds good for arbitrary point in the surface area X, therefore, we have 


M = 0 (14) 

Using the result (14) and integrating (2), (7), (8) and (9), we get 

[H n/ ] = 0 (15) 

[ w -]- M ^] = ° ' (16) 

kV 1 = 0 (17) 


Pirn VuM + [P'} n i -^ H lnl[H t ] = 0 


(IB) 


• * t_rZ 

where p* = p + — - and V t = v,- - Gn„ 
4rt 2 

G denotes the speed of Shock Surface, and 

Pm 


PlmYui 


1 2 P 

— V + c mv T + — — 

2 mv p m 


4tc 




= 0 


(19) 


The suffix nf is used to denote the dot product of the quantity concerned with the unit 
normal vector n,- to the shock surface 
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The eqn. (15) to (19) are the shock relations for liquid-gas mixture in the presence of magnetic 
field. Putting magnetic field intensity H = 0, we get the results for two phase non-conducting 
system which are in full agreement with previous results 7 . 

Propagation of Weak Discontinuity 

Let denote the surface of weak discontinuity moving in the homogeneous liquid-gas 
mixture with velocity G. Obviously the flow parameters are continuous across £(,) i.e., 


[ 23=0 ’ ( 20 ) 

while the first and higher order derivatives are discontinuous across it. Moreover, under the 
proper assumption of the regularity of the surface Z( ( ) and the existence of limiting values of 
the functions and their derivatives as they approach this surface from each side, it is well 
known that 




and 



( 21 ) 


( 22 ) 


where B is the normal derivative of flow parameter. 


For obtaining the jumps of normal derivatives of the flow parameters, we would like to 
use the equation of the state for the mixture in place of equation of energy, as former is 
sufficient for our purpose and latter is cumbersome to handle. Equation of state for the mixture 
under constant entropy can be written as 


P=P(Pm) 


which implies 



V/Af-cJ 


n 

dt 


+ ViP m , 


,>o 


(23) 

(24) 


with the help of (7), (24) gives 


dp 2 

+ V iP,i + c mPm v i,i = 0 


(25) 


Now taking the jumps of (2), (7), (8) and (25) and making use of first order compatibiliy 
conditions 13 , we get 

e„ = 0 (26) 


where 
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V„ l e i -H nJ \ i + Hfr„ / = 0 

(27) 

^/? + Pm^„/ = 0 

(28) 

Pm 3'n/^'i ■+■ ftfl/ + 4 ^ j ^ y tJ;- — H n j 6 ;J = 0 

(29) 

V n j\\. + p m c m \ n j = 0 

(30) 


Before solving (26) to (30), we observe that the situation is not exactly similar to its counter 
part for non-conducting flow 7 . Therefore, we follow a different approach. Eliminating 
unknown quantities other than i from (26) to (30) and substituting 

/, = //,///, we get 

V nlh = a ij^j ( 31 ) 

where 



nff i, 


22 

W 


4 *Pn 


S y + 


^m + 


4itp„ 


2 > 


n i n j 


4n Pm 


( l i n i + n jh) 


(31) implies 

vKa- -*(/} = 0 


which gives 




= o 


(32) 


(33) 


as Xj * 0. 

Taking (x,) parallel to Xj - axis and (/,) in the plane of (x t , Xj) at an angle 8 with x 1 and 
substituting 



(33) yields 
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+ b 2 l 2 nt -C 2 - b 2 ^j (b 2 l nI sin 0 j 0 

jV/„ /S inej {yli-b 2 ^ 0 

0 0 {vl-b 2 l 2 n i] 

which on simplification reduces to 


= 0 


(34) 


(35) 


Eqn. (35) is a third degree equation in V 2 ^, therefore, it will give three values of V 2 ^ and 
consequently three values of velocity of propagation of weak discontinuity. From (35), we 
infer 


2 2 2 
V *I = b l m 


2 1 I 2 

Y„/ = — i c m + b + 


(36a) 

(36b) 

(36c) 


The set of eqn. (36) gives magnitudes of three velocities of propagation of surface of weak 
discontinuity relative to the medium. If medium ahead is at rest ; they yield 


|c* + i> 2 + ^, + f. 4 + 6 2 ci(l-2l„ 2 ( )J 

vl - jU + » 2 - fl + o' + (l - i£,) 


Jl , 2 ,2 

G a - b l, 


nj 



(36a0 

(36b0 

(36c0 


where G a , Gf and G, are respectively magnitudes of Alfven, fast and slow wave velocities. 

A velocity is determined completely only when its .direction and magnitude both are 
determined. It is obvious that respective directions of Alfven, fast and slow velocities are 
given by eigen vectors of symmetric matrix (a,y). 

Determination of Eigen Vectors of Symmetric Matrix ( a,y ) : 
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Let Li be the eigen vector of the symmetric matrix (a,y) with principal value L. We obtain 
from (31) 

vl,Lih = a^Xj = LLjXj (37) 

which implies 

[vh-Ly.X^ 0 (38) 

2 2 
Consequently L = V nh That is, \ is an eigen vector of (a,y) with principal value V n/ . 

Since (a,y) is a symmetric matrix of rank 3, therefore, it possesses three mutually perpendicular 
eigen vectors. These three vectors are corresponding to three values of velocity of propagation 

a f s 

of weak discontinuity given by (36). Let us denote them by L f , L t - and L ; . 

Surface of discontinuity is the surface formed by 0 and <j> parametric curves. Its normal n 
is inclined with r 3 -axis at an angle 0. 

We obviously have 

h = (sin 0 cos 4», sin 0 sin <j>, cos 0) (39) 

in — (cos 0 cos<{>, cos 0 sin <j>, — sin0) (40) 

k = (—sin <f>, cos<j>, 0) ( 41 ) 


where m and k are unit vectors tangential to 0 and <|> parametric curves respectively. 

Now it is easy to show that £ is an eigen vector of the matrix (a,y) with the principal 
value as Alfvan velocity. Thus we have 

L- = k; for i = 1,2,3. ' (42) 


Since X,- is an eigen vector of the matrix (a,y) and eigen vectors are perpendicular to each other, 
we can write 


X f = (cos<)>, sin<|>, Af) 

In view of (43), eqn. (31) yields 

2 2 2 2 
g V n!~ b ~g OT Sin 0 

c m sin 0 cos 0 


(43) 


( 44 ) 
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hence the corresponding eigen vectors are 


L\ = cos<i), sin<(>, — 


V* f -b 2 -c 2 m sm 2 6 


c m sin 0 cos 0 


( 2 2 2 2 
s . V n ,.-b -c m sin 0 

L[ = cosij), ^ 

I c m sin 0 cos 0 


(45) 

(46) 


As it has been already stated, unique solution, for the velocity of propagation of surface of 
weak discontinuity in presence of magnetic field is not possible. That is why we have 
obtained three magnitudes of velocities given by set of eqn. (36) and three different directions 
given by (42), (45) and (46). Mere observation of eqn. (36), (42), (45) and (46) reveals that 

(i) The velocity of propagation depends upon the magnetic strength. 

(ii) The Alfven velocity is independent of velocity of propagation of sound in the mixture 
while the slow and fast wave velocities depend upon sonic velocity. 

(iii) In the absence of magnetic field, the velocity of fast wave reduces to te velocity of 
propagation of weak discontinuity while the remaining two velocities vanish. 

(iv) In the absence of magnetic field the direction of the fast wave coincides with the 
direction of propagation of wave. 
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Abstract 

The effect of suction on the temperature distribution and heat transfer in a uniform circular pipe 
taking small outward normal suction has been studied using perturbation method. The solution 
of energy equation in cylindrical polar coordinates has been obtained under suitable condition. 

It is shown that in the region of no -back flow, the maximum of temperature exists on the axis 
of the pipe, while in the region of back flow the surface of maximum temperature is near the 
wall due to suction. 


Introduction 

The solution for the temperature distribution of a fluid flowing in a circular pipe has been 
given by Goldstein 1 . The heat transfer of a fluid flowing between two parallel plates of which 
one is at rest, the other moving with constant velocity in its own plane is well Known 2 , 

Yuan and Finkelstein 3 have discussed the problem of laminar flow through a pipe of 
porous wall with injection and suction under the assumption that the maximum velocity of 
the Hagen-Poiseuille flow exists at the centre of the mouth of the pipe/ But it is found that 
the pressure gradient along the axis, at mouth of pipe, is not the same as the pressure gradient 
of the Hagen-Poiseuille flow. Verma and Bansal 4 have found the exact solution of the Navier- 
Stokes equations for the flow of a liquid between two parallel plates, one in uniform motion 
and the other at rest with uniform suction at the stationary plate. 

Bansal 5 has considered the steady flow of a viscous incompressible fluid through a circular 
pipe with small outward normal suction under the assumption that the flow is due to pressure 
gradient of the Hagen-Poiseuille flow at the mouth and that the suction at the wall is uniform 
and small. Moreover, he has used a perturbation method to determine the motion instead of 
assuming the form of stream function as done by Yuan and Finkelstein. He has also assumed 
that the radial velocity which vanishes in the Poiseuille flow has a finite magnitude except at 
the axis of the pipe where it vanishes. 

Verma and Bansal 6 have further discussed the effect of suction on the temperature 
distribution of the problem of Bansal 5 in a porous circular pipe with small suction at the 
porous wall. They have determined the solution of energy equation in cylindrical polar 
coordinates with the following assumptions : (i) The wall temperature is constant, (ii) The 
suction parameter is small, (iii) No diffusion takes place and (iv) Fluid is incompressible and 
heat conduction is negligible. 
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In the present, paper the effect of suction on the temperature distribution and heat transfer 
in a uniform circular pipe with small outward normal suction has been considered with the 
same assumptions as made by Verma and Bansal by perturbation method. It has been found 
that in the region of no-back flow, the maximum of temperature exists on the axis of the 
pipe, while in the region of back flow the surface of maximum temperature is near the wall 
due to suction. 


Formulation of the Problem 


Let us consider the steady laminar flow of a viscous incompressible fluid through a 
porous pipe of uniform circular cross-section with small suction at the porous wall. The 
temperature of the wall of the pipe is maintained constant. 


In cylindrical polar coordinates (x,r, 0), the equations of conservation of mass, momentum 
and energies are 


du 1 3 % 

— - + (v,r) = 0 

dx r dr 


( 1 ) 


du du 1 dp 

u — +v— = — — + VI 
ox dr p dx 


( a 2 ... ,2 A 

d u 1 du d u 

dr r dr dx 


( 2 ) 


3v dv 1 dp 

u 1- V — b V 

dx dr p dx 


d 2 v 1 dv d 2 v v 

37 + ^ + a7~7 


(3) 


f 


and 

pc P 

where 

<p = 2 


dT . dT\ 
U ~dx + V ~dr 


f -.2 


= K 


2_^ 


d T 1 dT d T 

v dr 2 r dr dx 2 y 


+ M> 


du) 2 fdv) 2 


dx 


j \* r j 


V 

Cr 


dv du 'j 
dr J 


(4) 


(5) 


x is taken as the axis of the pipe, r is the radial direction, 0 is the azimuthal angle and all the 
quantities are independent of 0. 

The boundary conditions are : 

r = R:v = v 0 , u = 0, T = T 0 , |^=0 

dx 


r = 0 : v = finite, u = finite, — — = 0 

dr 


( 6 ) 
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Since there is uniform suction, Bv/Bx = 0 i.e. v is independent of x and, therefore, it is a 
function of r only. Now, it is evident from eqn. (1) that B 2 u/dx 2 - 0. 

Thus eqn. (1) to (5) become 

Bu Bv v ,~ 

= 0 (7) 

dx dr r 


du du 1 dp 

U T~+ V T~ = — T" + ' u 
dx dr p dx 


d u 1 du 

a7 + ^ 


d v 1 dv v 


- ( dT BT) „ 


9 2 r _19T 3 2 r " 

3r 2 + r Br + dx 2 ’ 


+ 7 + 


Let us introduce the non-dimensional quantities as follows 

_ __ v _ _£ . £ - _ 7 

~ v 0 * Z_ R’ P "pu 2 

___ 'T’ 'T' TIP [±C 

0 = ~ — ~ , — (Reynold number), P r = (Prandtl number), ( 1 1 ) 

a = - (suction parameter) 

where U m is the maximum velocity of the Hagen-Poiseuille flow and T l is the finite 
temperature at the axis when there is no suction such that 


r, = To + - 


Hence eqn. (7) to (10) in non-dimensional form are 


du G f dv v 


3x l 9t| t| 


+ - =0 
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35 ct 35 

u + — V 

dx R e 3q 



R. 


( -2- 
3 u 

3? 


1 35 

ti an 


V 


3v 


R £ 3 p 1 ' d*v 1 3v v 

nan - n ; 


R e p r 


( . 30 a .39 
u — + — v— 
dx R e 3r] 


3 2 e i ae 8 2 e 

3q 2 n SR 33c 2 





'3uf 

3q y 


(14) 


(15) 


(16) 


The boundary conditions (6) with the help of eqn. (11) are 

T| = 1 : v = 1, 5 = 0, 6 = 0, ^ = 0 

3n 

30 

•q = 0 : v = finite, 5 = finite, — = 0 

3q 

Let 

p(x,f[) = p 0 + p'(x,r\) 
u{x,r\) = n 0 + u'(^,R) 
v(n) = v'(T|) 

and 0(3c,Ti) = 0 O + 0’(i,q) 


(17) 


(18) 


where the primed quantities are the perturbations caused by the suction and p 0 , u 0 , 0 O are the 
known quantities for flow when there is no suction (i.e. Hagen - Poiseuille flow in a circular 
pipe), satisfying the equations 


dp o 

3q 


= 0, — ~ = 0, 

dx 


dpQ = ao=: 1 f d 2 “o I 1 d“o^ 
dx R e ^ 3p 2 q 3q ; 


and 


3 2 0 o 1 

o *7" *“ “ 


30r 


3n 


r\ 3q 


+ 4 


duo 

3q 


= 0 


(19) 


Here we have p 0 independent of tj 
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«0 = a(ti 2 -1) where A = ~^- = — and 0 o = b(i-T| 4 ) 


u */ 4 dx 4 u “ T 7 • 

Substituting eqn. (18) in eqn. (13) to (16) and after simplification with the help of eqn. (20), 
we get 

du' a f 3v' v 1 ') 

__ + __l__ + _l = 0 nn 


dx R e l t\ 


, r 3m’ ^ ct ,f duQ du'} dp’ 1 d u' 1 du' 


(“0-1 ^ * 


dx ) R e ^ 3ti 3q J dx R J ^q 2 r( 3r( 


1 3v’ R e dp’ 1 3 v’ 1 3v' v' 


a ^ I 3n n 


. , . ,39’ CT ( 30 o 30'^ 3 0' 130’ 3 0' 

R e P r (“0 + u + TT v = — 2 " *i ^ — 1 T 

1 v ' K l JJ a^ 2 3n 3* 2 


du''] 2 afdv'f ct 2 f v' - ) 2 


44 2 » + TT 


+ 3hr + 


j 

+ B^-12q 2 j + -^£^-4q 3 j + 4(-A.2Tl)'‘ 


3u ,s \ 2 3u 0 du' 

— + 2-—^ — 


3ti I 3rj 3ri 


From the last eqn. of (19) and (20) B = A 2 . 

The boundary conditions (17) with the help of eqn. (18) are 

11 = 1, v’=l, u'= 0,~0, 0'= 0 
3q 

^A1 

T| = 0, v'= finite, u’~ finite, —— = 0 

3q 

The solution of eqn. (21-23) with the help of the boundary conditions (25) is 

«’= -^*[/'(^) + r/( 1 i)"l + F ( 1 i)- v ‘= /Oi) 
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r . . 3 13 7 5 3 1 5 1 7 

f(T|) = — ti — ti + <7 ri ri + — ti ri 

JK U 2 1 2 1 144 1 48 ' 16 1 144 1 


2 f 103 101 3 13 5 5 7 1 9 1 11 

+cr ti Tl + ri ri + ti ri 

9600 4320 864 1728 1920 43200 1 


1 6 14 32 47 1 2 x f 5 32 I 4 


^(Tl) = a — + oA — n — ti + — ti + o — n +~n 

■ ’ R e L 36 8 ' 4 1 72 J R e { 8 4 1 8^ 


2(1 8 1.6 9 4 79 ^ 2 Y 28337 37 2 143 4 

+<7 [256 71 ~ 24 11 + 64 11 ~ 786 J + ^ 86400 + 72 ^ ~ 576^ 


59 6 Is 1 10 

+ TI Tl + Tl 

864 1152 1 7200 1 


f° = b-b 3 ’ 


and so on 

v r\ t? ■* . Y 1 6 1 4 ^2 47 A 

fo . 0 ,F, = - + ^ + 

and so on. 

Now, after substituting the values of u 0 , v 0 , u' and v' from eqn. (20), (26), (27) and (28), the 
eqn. (24) becomes 

= ^ + 1 + 4 [i£! f f <2 + J_ f 2 + 1 ff ■.} 

3n 2 T| 9 t[ 3ic 2 Rj (/ 


\<JX f f f" 


if if-t-t n (*-f- f j-'}j 

Solution of Equation 29 

Let 9-9fa) + *V(R) + * 2 S0n) (30) 

where <pCn), yCn) and £(n) are unknown functions of q to be determined. 

By substituting eqn. (30) into eqn. (29) and equating the coefficient of like powers of x, 
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we get the following set of equations 


5"+^'+^- ) +2aP r rt A,/t-o«/ = 0 

B K e T 1 J T 


11 R 2 l^ri 2 T1 1 


0 R 


Y' 2 a / /' /") 

V" + „ + 4 — 4-U-T- F' 


B 1 


/’ r / r r 

-T-T ^ + W -o- — - — 


ri 1 


-<W +aP r f y + -^/ y - 2f? e P r ^7? e a 0 (q 2 -l) + T 5 = 0 


... 1 », 16a 2 f/' 2 , / 2 //* 


cp"+ — cp + — s- +4T" i +4R.a 0 TlF’ 

11 Rt 1 1 T1 2 T1 J 


il" r\ 


B 4 GP r r?f -oP r N-R e P r |I^ fl0 (q 2 -l) + F|Y + 25 = 0 


The boundary conditions (25) with the help of eqn. (30) become 
T1 = 1: SOn) = 0; Y(T) = 0; <p(il) = 0 

, = o:|i»o ; |r=o ; |£=0 ( 34 ) 

3r| 3 ti 3q 

Now, the solution of eqn. (3 1) can be expressed for small values a of by a power series 
developed near cr = 0 as follows 

$ = $0 + cZ,l + G% +-0 H $n ■ (35) 

where 5„'s are taken to be independent of a. 

By substituting eqn. (35) in eqn. (31) and equating like powers of cr, we get the following set 
of equations 


£o + -i;o = 0 
T 

£i + ^o + 2P r /o + ~j^o = 0 

^2 + “ ^2 + 2P rf 0 + 2P r /l ^0 + ~~ fo £l + ~~fl ^0 ~ P r fo £l ~ P r fl ^0 

b tj 

4 (j r rf_ n 


\rf ii 1 
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The boundary conditions to be satisfied by the $„'s are 

WO) = 0, Wl) = 0 for all n (39) 

Now from eqn. (36) to (38), we get 

$o = 0, $i = 0 and $ 2 = -l(l-q 4 j 

Hence the second order perturbation solution in 0 of (3 1) is 
.2 

$( T l) “ ~ ~ T 1 j (40) 

Now the solution of eqn. (32) can be expressed for small values of 0 by a power series 
developed near 0 = 0 as follows 

V = V 0 + 0 Vi + 0 2 y 2 +. . .0>„ 

where y n ’s are taken to be independent of 0 . 

Substituting this power series in eqn. (32) and equating the coefficients of like powers of 0 , 
we get the following set of equations 

- 1 ■ 

Vo + -Vo = 0 (41) 


" 1 ■ 

V1 + -V1 + 4 


/ 

fo /o /o 

y " n " 1 



= 0 


(42) 


V 2 +^V 2 -/ 5 r/oVi + / , r/o¥i + ^/oV 1 -2R e P r ^ ao(r» 2 - 1) -^(l-q 2 ) 

+ 4a 0 Tlf-4- - — - •=7"] + Trf^T - — - — Vl = 0 

r n 1 J R e l rf n 1 J 1 


Proceeding exactly as before, the second order perturbation solution in 0 of (32) is 
Vfo) = - V) + ~ ^(95 - 108 t| 4 + 16t| 6 - 3q 8 j 

+12a o( 3 - 10q 4 + 8q 6 - q 8 j + a 0 P r |307 - 288q 2 + 36q 4 - 64 q 6 + 9q*| (44) 



TEMPERATURE DISTRIBUTION IN LAMINAR FLOW 


203 


Similarly we proceed to find out the solution of eqn. (33) by substituting the power series in 
cp, as usual, we obtain the second order perturbation solution in <y as 

2 

<p(n) = <po + <*Pi + <P 2 (45) 

where 

<P° = °> <Pi = (95 - 108ri 4 + 16 b 6 - 3b 8 ) + (l3 - 16q 6 + 3 tj 8 ) 

r>2 J2.p 

+ — ~ - (101- 144r( 2 + 36ri 4 + 16ri 6 - 9ri 8 ) 

Xjirr 

and 


1 / 2 4 s'! 

< P2 = — 2 77-99ti 4- 27b -5b 
9R e V J 

R 2 2 . . 

+ ~m ( 1440 )( 72449 " 103050t 1 4 + 39400 1 ) 6 - 9675t) 8 + 936 V° - 60r| 12 J 
-^(23 - 27T] 4 + 4ti 6 ] + ^(77 - IOOtj 6 + 25^ - 2r| 10 ) 

+ lf(i2^oo)( 5589 “ 25200116 + 303751,8 _ 1 1664 V 0+900 t i 12 ) 

+ -^-1(293091 -35640011 2 + 95400b 4 - 2200b 6 
576 V 72007' 

-1 1475B 8 + 1584b 10 - 200b 12 ) 



22 2 , , 

| P r Re a 0 ( 1 

2304 1600 


103764 - 122800b 2 + 16300b 4 + 1600b 6 


8 10 \ 1 ~ 

44300b -3664b + 500b 


The temperature distribution obtained with the help of eqn. (18), (20) and (30) is 
e = (^l_q 4 j + ( p(- n ) + x H r( T1 ) + i 2 ^( T1 ) (46) 

where 4(B). Y(B) and 901) are given by the eqn. (40), (44) and (45) respectively. 

Now by keeping A = -1, we have 5 = 1 and R t a 0 = - 4. Further by putting R t = 1000, 
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0 — >■ 

Temperature distribution in pipe flow without suction. 



.2 .4 .S .8 1.0 


0 


Fig. 2 - Temperature profiles 9 - tj. 
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P T = 5 and a = 0.1 in eqn. (46), we get after simplification 

6 = 1. 727 - 0. 702q 2 - 1. 0049-n 4 + — ^ (- 0. 523 + 0. lOt) 2 + 0. 4 108r| 4 ) 

+ - —— *■ (0. 04 - 0. 04t] 4 ) 

(I000) n > 

The temperature profiles (9 vs. rj) are shown in Fig. 1 for various values of x. 
Similarly by putting B = 1, P r = 5, R t ~ 1000 and x = 3000 in eqn. (46), we get 
0 = 1 - t) 4 + ct(- 6. 52 - 5ri 2 + 1 1. 75rj 4 ) + a 2 (l7. 1 15 + 9. 8 It) 2 - 30. 74iq 4 ) 


(47) 


(48) 


The temperature profiles (0 vs. q) are shown in Fig. 2 for various values of a. 

Discussion 

The temperature profiles in Fig. 1 are divided into two parts : (a) The region of back flow 
and, (b) The region of no-back flow. 

It is observed that the maximum of the temperature distribution exists on the axis of the 
pipe in the region of no-back flew, whereas the surface of the maximum temperature 
distribution is near the wall in the region of back flow. 

The temperature on the axis of the pipe decreases and remains always positive in the 
region of no-back flow, whereas there is a continuous rise in temperature after x = 7000 in the 
region of back flow. 

The effect of the variation of suction parameter <y is shown in Fig. 2. It is seen that the 
temperature decreases with suction parameter a in the region ®f no-back flow. The temperature 
profile without suctipn is also drawn for comparison. 

- 4 

Proceeding to the limit cr = 0, we have 0 = 1 — nq as in the case of Hagen-Poiseuille 
flow. 
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Abstract 

Local similarity solution of the voiticity equation for flow over a sinusoidal wall for both low and 

large Reynolds numbers has been presented without using the perturbation method and the 

linearized boundary conditions. 

Introduction 

The flow problem over a wavy wall is important in connection with the study of the 
effect of inhomogeneities on boundary layer separation, aerodynamic heating and the surface 
distribution of the stresses created by the fluid flow. It is ^lso useful in investigating the 
hydrodynamic instability of a protective layer of fluid injected as a coolent 1 . Inspite of its 
importance, it has attracted the attention of only a few research workers. Bejamin 2 was 
probably the first to consider this problem. His analysis is based on the assumption that the 
basic flow (defined in the absence of waviness) is parallel. Singh and Lumley 3 studied the 
effect of an isolated component of wall roughness on the average velocity distribution. Lessen 
and Gangwani 4 studied the flow stability under boundary-layer approximation. Other studies 
include the Rayleigh problem 5 and the free convection problem 6 . These studies 3 ' 4 ' 6 were aimed 
at seeking a perturbation solution about a basic flow which is either a function of y (co- 
ordinate transverse to the flow direction) only 2 - 3 or an average of the flow quantities over a 
wave length with similarity transformations incorporated 4 or fully developed flow 6 satisfying 
the linearized wall condition 3 ' 4 - 6 . The boundary-layer approximation holds true for large 
Reynolds number (based on the axial distance and the free stream velocity). For low or 
moderately large Reynolds number describing the early flow region, the perturbation solution 
about a y-dependent basic flow will be inappropriate since it has been well established 7 - 8 that 
such basic flow also depends on the axial co-ordinate (x). As such, a perturbation solution will 
involve a basic flow which is function of x and y. 

To the authors knowledge, no attempt has been made earlier to study a steady 
incompressible flow over a wavy wall for both low and moderately large Reynolds numbers. 
To study this problem, one is tempted to seek a perturbation solution about a known flat plate 
solution, which is a function of x* and y*, satisfying the linearized wall condition. These 
complicate the perturbation method, which is an approximate one. Also, the use of linearized 
wall condition provides information regarding flow at the mean position and as such Singh 
and Lumley 3 have suggested that the solution with linearized boundary condition is not valid 



208 


J. DEY and G.NATH 


near the wall. In order to circumvent these difficulties, it would be appropriate to solve the 
Navier-Stokes equations satisfying the actual boundary conditions prevailing on the wavy 
surface. It appears that this is possible via some simple transformations, which permit one to 
obtain exact solution of the Navier-Stokes equations for both low and high Reynolds 
numbers. Such results have been presented in this paper. 

Governing Equations 

We consider a two-dimensional steady laminar incompressible flow over a wavy wall. Let 
u* and v* denote the two velocity components in the x* and y* directions, respectively 
(Fig. la). The wavy wall is defined by 

y^ = e * sin(27cc * A) ( 1 ) 

where e' is the waviness amplitude and X is the wave length. Following dimensionless 
quantities are introduced to obtain the non-dimensional vorticity equation 

x=x*U/ 1> ; y = y*U /‘ o 

u-u/U ; v = v/I7 

P = p*/U 2 ; e = e U/v ( 2 ) 

R x = XU/x> 

where U denotes the uniform free stream velocity, p is the fluid density, v is the kinematic 
viscosity and p* is the dimensional static pressure. The non-dimensional form of the wall is 
y w = e sin( 27 tX//? x ), where R\ is the Reynolds number based on the wave length X. 

The governing non-dimensional vorticity equation is 

+ = ViP (3) 

where suffixes x and y denote partial derivatives with respect to x ^nd y, respectively. 

2 ( 2 2 2 2 \ 

P(=v* - Uy) is the vorticity and V { I = 3 / dx +3 / 3 y I is the two-dimensional Laplacian 
operator. 

The boundary conditions are 
w = v = 0 at y = y w . 

u -* 1 as y -» «> ^ 

It has been assumed that the disturbance dies out as y — > which is of course obvious. 

We define the non-orthogonal co-ordinates (X, F), similar to those of Shankar and Sinha 5 , 
such that 





Flow configurations (a) How configuration considered, (b) and (c) other possible flow 
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X = x;Y = y-e sin(2juc / R 'jJ 


(5) 


which set the range of Y as 0 < T <», i.e., the flow geometry has been transferred to the 
positive Y plane. The vorticity eqn. (3) then becomes 




( 6 ) 


where <p is the non-dimensional stream function and V is given by 


V7 2 3 

V = — + 

dx 


jl + (e 2 jc / R x f COS 2 (2itX / f? x ) j 


3 Y 


-2e(2tc / %) cos(2jtY / 
+e(2n:/i? x ) 2 sin(2jtX'//? x ) 


3 Y dX 


3 Y 


(7) 


The solution of (6) being intractable in the present form, we define the parabolic 
co-ordinates (a, r|) such that 


X + i[y + e sin(2itY / %)] = (cr + iq) 2 / 2 


±sX + 


^X 2 + (F + 6Sin(27tX-/%)) 2 

v 

r\ = |-Y + [^Z 2 + (r4-esin(27t//? Ai )) 2 j 


, 1/2 


( 8 ) 


since B is positive, the jpositive and negative values of c correspond to the positive and 
negative F-plane, respectively. In the parabolic co-ordinates the vorticity equation can be 
solved in either positive or negative half plane alone. Since the sinusoidal function 
encompasses both positive and negative y planes in ( x , y) co-ordinates, one might think that 
the parabolic co-ordinates (cr, r|) are not appropriate in such situation. But then it must be 
remembered that due to the transformations (5), the present problem has been transferred to the 
! positive F-plane (0 <,Y< ~=). As such, the problem of positive or negative y plane does not 
l^rise and the use of these co-ordinates are appropriate. The vorticity eqn. (6) in (cr, T|) 
co-ordinates is 


^ca+^n + <Pa^ti-9ti^o = 0 

where W = (<p co + <p nn ) / ^ 


( 9 ) 
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It is interesting to note that the complicated vorticity eqn. ( 6 ) has been reduced to a very 
simple form ( via a lengthy process involving differentiation and large algebra), similar to that 
of Davis 7 for flat plate. However, in the present case, the lower limit of n(by setting Y = 0) 
being 


no = 


= i-X + 


2 / 

X +(esin(2rcX / R x )j 


, 1/2 


( 10 ) 


is non-zero unlike that for flat plate where Tl o =0. It may be noted that the present approach is a 
departure from the widely used approach of using a linearized wall 3 - 4 ’ 6 . 

The corresponding boundary conditions are 

cp = <Pn = 0 at n = no 

9 ~»crq; 9 ^ -» a as q ^ 


These are dictated by the conditions that u - v = 0 at 7=0 and u— »1 as 7-»«. The boundary 
conditions ( 1 1 ) dictate that the form of 9 be 


9 = afCn) 


( 12 ) 


This is equivalent to say that 9 has been expanded about <J 0 (by setting 7 = 0 in ( 8 )) on the 
surface, similar to that of Davis 7 , and that only the first term of the series has been retained. 
The vorticity eqn. (9) for local similarity at any location becomes 


/ 


f- 


4q 


2 2 

tfo + n 


/ + 


1 

~2 2 

a 0 + n 


(°l-n 2 )f -infY 


0 


(13) 


where prime denotes differentiation with respect to q. 
The corresponding boundary conditions are 


/=/' = 0 at q=q 0 


/-> q ; /' 1 as i]->« (14) 

It may be noted that the eqn. (13) is similar to that for flat plate 7 8 except that the lower limit 
of q is non-zero in the present case and that o 0 is not equal to (2X) W . For large a 0 , i.e, 
representing the boundary-layer approximation, one can obtain the local similarity solutions. 
Unfortunately, the local similarity transformations for the boundary-layer equation are not 
known though the equation can be obtained from (13) for a 0 — > ». It appears that a separate 
search for the local similarity transformations for boundary-layer equation needs to be made. 

The numerical solution of the vorticity eqn. (13) under boundary conditions (14) can be 
obtained provided that q 0 and <r 0 are prescribed at each axial location on the wavy surface. In 
other words, the flow field can be obtained following the wall waviness. 
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Results and Discussion 

Quasi-linearization and interation technique 9 were used to solve (13) for given ti 0 and ct 0 
and for e= 0.3; 10 <X < 41 and R% = 1. This range of X covers low to moderately large 
Reynolds numbers. The solution was built-up with one particular solution and two 
complementary solutions satisfying the two conditions at infinity. The forward integration 
was carried out using Hemming’s predictor-corrector method on DEC 1090 computer. The 
integration step size was optimized to obtain accuracy up to 4 decimal places. A large number 
of computation is involved for unit wave length and unit change in X. This involves very 
large number of computations for a wide range of Reynolds numbers (X). As such, detailed 
computations were carried out for 40 <X < 41. Though computations were carried out for 
X = 10 also, results are presented for 40 < X < 41 for the sake of simplicity. From the flat 
plate results 7 ' 8 , we notice that the variation of shear stress for unit change of Reynolds number 
is very small. As such, in the intermediate range any noticeable variation will occur in third or 
fourth decimal place. 

Representative axial velocity (u) profile has been shown in Fig. 2 for X = 40.6. This has 
been selected because X = 40.6 corresponds to negative value of e sin (2 tcX). Comparison 
with the flat plate velocity profile is also shown in Fig. 2 (To show the velocity difference 
clearly, a part of the velocity profile has been shown. It is understood that the velocity profile 
attains the free-stream value as r \— »°°, not shown in Fig. 2). For other X, representing 
positive value of e sin (2 nX), similar trend is observed. It may be noticed that though T | 0 
values are different the axial velocity is less than the flat plate axial velocity. Reduction in the 
mean velocity due to waviness was also observed by Singh and Lumley 3 , and Lessen and 
Gangwani 4 . In the present case, we find that the local velocity is also reduced due to waviness. 
For Stoke’ s flow through a pipe having sinusoidal corrugation, Phan-Thien 10 observed 
decrease in flow rate and increased pressure drop compared to non-corrugated pipe. The decrease 
in flow rate could be attributed to the decrease in velocity since the average areas are same in 
these two geometries. These suggest that the present analysis is appropriate and that the flow 
. field can be defined for both low and high Reynolds numbers flow over a wavy wall. 

The non-dimensional wall shear stress can be expressed as 

T = x w (X / 2) 1/2 = 2u(v** + uJ,)(X / 2) 1/2 / U 2 

= 2/"('no)(X/2) 1/2 cro(<yo-Tio)/( c7 o + 1 lo) 2 (15) 

where suffixes ;c* and y* denote partial derivatives with respect to x* and y*, respectively, and 
it is shown in Fig. 3. Factor (X/2) 1/2 has been introduced for comparison with flat plate 
geometry. As mentioned earlier* any significant change in the wall shear values occurs in third 
or fourth decimal place (to find the shear stress at any point, 0.5 needs to be added to the 
values given in Fig. 3). Fig. 3 shows that the shear stress attains maximum and then 
decreases to the flat plate value. This pattern is expected for the positive y-plane. In the 
negative y-plane i.e., for negative e sin (2nX) (40.5 < X <41.0) we do not expect negative 
shear stress, as this will mean a flow reversal. Also, physically we find that the fluid contact 
surface in this region (e sin (2tzX) < 0) is more than the flat plate and hence a higher friction 
is expected. Hence, the positive higher wall shear in this region is justified. 
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Fig. 2 - Non-dimensional axial velocity profiles. 
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Fie. 3 - Non-dimensional wall shear stress variation with the axial distance. 
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From the above discussion, we find that the local similarity solution of the vorticity 
equation for wavy wall, under consideration, is possible and that such solution gives more 
accurate results than that of the perturbation method. Also, it is expected that the local 
similarity solution for flow over wavy walls (see Fig. lb and lc) in which the waviness is 
confined to half the wave length should be possible via transformations (5) and (8). 

Conclusions 

To sum up, we find that (i) the local similarity solution of the two-dimensional vorticity 
equation is possible for flow over a wavy wall and for a wide range of Reynolds numbers 
without using the perturbation method, (ii) the local axial velocity is less than that for flat 
plate, (iii) the Wall shear is more than the flat plate geometry, (iv) the local similarity 
solution of the boundary-layer equation is also possible, but the similarity transformations 
have to be found out. 
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Abstract 

The differential equation for bending vibration of a slender beam including the effects of rotatory 
inertia and shear force in a centrifugal force field on Pasternak foundation has been obtained. A 
method based on Gale rk in’s principle has been employed to find the upper bound of fundamental 
mode of vibration. The accuracy of this upper bound is investigated by comparing the results 
when the effect of Pasternak foundation is not taken in consideration. 

Introduction 

The analysis presented in this paper considers vibrations of a slender beam including the 
effects of rotatory inertia and shear force that could represent a turbine blade of simple 
geometry. The beam is attached to a disc of radius r 0 and disc rotates with angular velocity Q , . 
It is assumed that the shear centre of each cross-section coincides with the centre of gravity; 
consequently only transverse bending vibrations will occur. 

For the analysis, it is assumed that S and I which are respectively the area of the cross- 
section of beam and moment of inertia of the cross-section about the axis perpendicular to the 
plane of bending are constant. This, however, is not an essential assumption. If S and I are 
functions of x , the amount of computing work is greater than when S and / are independent of 
x. 


The Differential Equations 

The differential equation 1 for transverse vibration of a beam including the effects of 
rotatory inertia and shear forces is 


-A 

dx 4 


£/ ^ +p5 f^“ p ( : 


1 +- 


3 4 y p 2 / 3 4 y 


K'G) dx 2 dt 2 K'G 3r 4 


( 1 ) 


where E is the Young’s modulus of elasticity, p the density of material, G the modulus of 
rigidity for the material of bar and K' is the ratio of average shear stress to maximum shear 
stress over a cross-section. 

/ 

The differential eqn. (1) on elastic foundation takes the form 


n 4 3 2 

ay 3 y 

EI^ + pS—t 

dx 3 t 


pi 


1 + 


3 4 y 

2 . 2 


K’G) 3xV *'G d[ 


2 r -.4 
P / 3 y 

+ —~ + q = Q 




( 2 ) 
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(3) 


where q is foundation reaction. 

The general form for a Pasternak foundation 2 in two dimensional case may be written as 

/ \ 3 2 y 

q(x,t) - Ky - Gq — % 

dx 

where K is Winkler foundation modulus and G 0 is shear foundation modulus. 

The eqn. (2) represents flexural vibrations of a beam including the effects of rotatory 
-•inertia and shear force on Pasternak foundation, this equation is modified by adding the lateral 

load, — if the centrifugal force effect is to be considered and then it becomes 


dx 

a 4 a 2 

ay d y . 

El— + pS— -p/| 1 + 
dx dt 


a 4 2 r a 4 a 2 a 2 

° y , P 1 ° y . „ . 3 j d m 

K ' G )dx 2 dt 2 K'G dt 4 + y-^ 0 ~-— = 0 ( 4 ) 


dx dx 


where 


d 2 M 

dx 2 


= p sn‘ 


Determination of Natural Frequencies 


(5) 


We now put eqn. (4) in terms of a dimensionless variable by introducing \-xjl. It takes 
the form 


3 4 y r 2 / 2 P 


d^ 4 E/pI/Sl 




l 2 

fi+— 

a 4 > , I 1 l 2 d 2 y 

r l 2 ] 

E / pi 

l AT'gJ 

'd% 2 dt 2 E/p’ns’dt 2 

l*/pj 


-^ 2 4=o 


El’ d? 

The solution of eqn. (6) is of the form 
y = T sin cor 


( 6 ) 


(7) 


where co is the natural frequency of vibration. The function Y satisfies all the boundary 
conditions of the beam which are as follows. 
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At the fixed end i.e. at t, = 0 

(i) deflection due to lateral bending and lateral shear are zero i.e. 7=0. 

(ii) <j> = = 0 where yt, is the lateral bending deflection. At the free end i.e. 

dx dx 

at f = 1 

d 3 2 d 2 y 

(a) bending moment = 0 i.e. — = 0. This implies that pS — j ~ K'GS — 2 = ® 

dx dt dx 


3 3 

3 y 3 y n 
(b) shear force is zero i.e. — - <}> = C . This implies that El — j - pi — j - 0 

dx dx dt dx 


Making the substitution of the solution (7) and introducing £ = x/l, we get the four 
conditions which are 


(0 y = 0 ' 

("> I " 0 


atj = 0 


2 ICG try . 

(iv) 

d% p r d 


at 2; = 1 


( 8 ) 


Substitution of (7) in (6) gives 


d A y O 2 / 2 l 2 


dt? E/p I/S 




0 ) 


2,2 


LG 


E/pt 


l 2 U 


d 2 y co 2 / 2 I 2 ( 




l 


Kl 

H y =75- 

EI El d £ 2 


K'Gjdt, 2 E/p’l/S 
G 0 l 2 d 2 y 


CO 2 r E 77 S' 
E/p’K'G’ l 2 


|y 


= 0 


( 9 ) 


Now we make use of the following substitutions 
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y „ __e_ K _f_ K a 2 / 2 

E/ p’ 1 K'G' Kl I/S’ 3 E/p 
Then the eqn. (9) becomes 


,K 4 = 


Q 2 / 2 _ gp/ 2 

£/p’ 5 £/ 


( 10 ) 


d*y 




4 -k 3 k 2 \ 




+ X(l + Ki)^ 


ft 


■XK^\-&-X\t + K 4 y-K 5 ^ = 0 

K 1 / at. 


dl 




( 11 ) 


For an approximate determination of the fundamental frequency y is chosen as the shape 
function in a series form. This function satisfies the boundary condition (8) and is given by 


y = \ 2 +C^ + D^ 

From the boundary conditions at £=1, we have 


( 12 ) 


N 

2 K'G 1 


/ \ 
6K'G 


12 K'G N 

1 "b 2 2 

+ 

1 2 2 

C + 

1 2 2 

CO p/ ; 


ca p/ J 


i, « p/ J 


D = 0 


(13) 


2 + 


3 + 


6E 

2 ,2 . 

® Pi ) 


C + 


\ 


4 + 


24£ 

~~2 2 
G> Pi J 


D = 0 


(14) 


Solving eqn. (13) and (14) for D, we get 


1 f 1 IE 

D = - + ) + 

6 1 18 12 K' 


„ 2 2 
E 1 <0 l 


1 13 E 

■ + ■ 


G) E / p 54 216 K'G 32 


r £ V 

K'G 


( 2 2 ) 2 
CO f 


£/p 


(15) 


From eqn. (14), we get 


C = + ■< — - — 


1 1 £ 


2,2 

CO / 


1 


7 £ 1 

• + — 


r e \ 


3 L 9 3Ar'GJ£/p 18 54 K'G 18 


K'G 


( 2,2 V* 

CO / 


£/p 


(16) 


On substitution, C and£> become 
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' C = -- + ■[- - - X'jlx +{— K x + -i^i 2 U 2 
3 [9 3 j 118 54 1 8 1 J 

D = - + j- — + -L Ajk + j- iTj - jtflx 2 
6 l 18 12 M I 54 216 1 32 J 


The terms only upto X 2 have been retained in (17). 

If we are satisfied with the trial solution (12), which is sufficient in most of the cases, we 
must form the Stieltjes integrals from Galerkin’s rule i.e. 




,2 _ 

+ (X + LVj - + (-Xfsf 2 + K^ 2 + K^y > ydc, 

Now using (12), the Stieltjes integrals are evaluated as follows 


(M*f 


1= rra + l_fa5 + i 

• J L / 2 W 2 


& i ]\ d y 




: ?{H c+ ii 0+ ? c2+ ! CD+ )j c2 


+ fe + ^C + |D + ^C 2 + |cD +— D ; 
115 3 5 35 8 21 


^3 = j(7 + ^^ = ^{i + C + D + ic 2 +CD + iD 2 


+\\+^c'+^d+-c 2 +-cd+-d 1 

15 6 7 7 5 9 


/< = i^ = ! +2C+ T° + 7 ci+3CD+ T° 2 


/5 = j3- 2 ^=i + Ic 2 +|D 2 +|c + |D+i-CD 


( 19 ) 
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2 2 

Since in these expressions C ,D and CD are occurring, we evaluate them as given below 

= V + f-7" + “ *i - — ^i 2 V 

9 9V3 V l 81 81 1 18 V 


D 2 = — + 1 


36 181 3 2 


1 1 


+ 7:%! - 


l_ 4 

324 4 648 ' 


rKt 


1 Kf 


288 


f 5 Jd_ 1 2 

1 32?" 324 ^ 1+ 72^ 


C£) = -I+fi--iK 1 V + ( 

9 U 8 9 l J \ 

Eqn. (18) takes the form 

h - K l K l{h - h) + {X + XK x - K s )j 4 + (~XK 2 + K x X 2 + Kijjs = 0 


or 


J\ — K 2 K 2 (J 2~ Jj) ~ KkJ 4 + K 4 J 5 + /-|(1 + K X ).I 4 — K 2 J 5I + K^J = 0 (20) 


Since all /’s given by (19) are functions of C and£>, they can be expressed in terms of X with 
the help of expressions (17). Substitution of all J’s in terms of X in (20) and retaining terms 
only upto X 2 one gets an equation of the form 

R — QX 4 - PX — 0 ( 21 ) 

where P, Q and R are constants. 


The solution of eqn. (21) is 


X = 



Q 2 - APR 
~2P 


( 22 ) 


The right hand side of (22) is positive since it may be shown that Q 2 _ 4 PR > 0. 

The smaller of two values of X given by (22) is an upper bound for the frequency of the 
fundamental mode of vibration, while the larger of the two values of X corresponds to the next 
higher mode of vibration. The value of X so obtained from (22) includes the Pasternak 
foundation effect. If this effect is ignored the values of X can be computed by omitting these 
terms contained in eqn. ( 20 ). 


Numerical Example 

A numerical example for vibrations of a slender rotating beam including effects of 
rotatory inertia and shear force is now presented. The frequencies are computed from eqn. (22). 
The physical constants of the blade are taken as follows 
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E/G = |, Cl = 314s -1 , AT'=|,E = 30xl0 6 lbs/in 2 ,/ = lOin, 
p = 0.28 lbs/in 3 , //Q = 10, ^ = 3, ^4=V^5 = 5 > 

These values of -v /^4 and yjK^ were taken from Wang 3 . With these values, the various 
parameters and the constants used in the frequency equation are as follows: 

/ = Sr 2 (r is the radius of gyration). 

K 1 = 4 ,K 2 = 100, K 3 = .0920229, K 4 = K s = 25 

C = - .66667 - 1.22222 X + 1.5370 D = .16667 + .27778 X - 2111 8 X 2 

C 2 = 0.44444 - 1.62163 X - .55556 X 1 , D 2 = 0.027778 + .09259 X - .01543 X 2 

CD = - 0.11111 - .38889 X + .10185 X 2 , = 0.80001 + .33333 X - 1.68520 X 2 

J 2 = 0.09259 - 3.19372 X + 1.58450 X 2 , 7 3 = 0.46790 - 8.08350 X + 4.1 1377X 2 ' 

/ 4 = 0.04762 - 4.63015 X + 1.90872 X 2 , J 5 = 0.06421 - .64777 X + .37736 X 2 

= 0.25684 - 2.59108 X + 1.50944 X 2 , 

^ 6.18290 - 41.62625 X + 28.19240 X 2 , J.-K-JC, (J 2 -J 3 ) -K^+K.J, = 
4.66847 - 82.94784 X - 16.69412 X 2 

Substituting these values in (20), retaining only up to second order terms of X, we get 

4.66847 - 89.13074 X + 25.18897 X 2 = 0 (23) 

The two values of X for this equation are 

Xj = . 05318, Xj = 3.48531 (24) 

which are the upper bounds of the frequencies corresponding to the first two modes of 
vibrations of rotating blade including the effects of Pasternak foundation. If the effect of 
Pasternak foundation is not taken into consideration, the equation in X is 

4.25372 - 50.79006 X + 63 .47306 X 2 = 0 (25) 

which gives the smaller value of X i.e. corresponding to the fundamental ijiode of vibration as 

Xi = .09505 • (26) 

Concluding Remarks 

On comparing the values of X[ from (24) and (26) it can be easily stated that the 
fundamental frequency of rotating beam is considerably decreased, when it is placed on 
Pasternak foundation. 
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Abstract 

The motion of visco-elastic liquid (Kuvshiniski type) with uniform distribution of dust particles 
down an inclined plane under the influence of exponential pressure gradient has been studied. The 
velocities of liquid and dust particles are obtained in elegant forms. The velocity profiles of liquid 
and dust particles have also been plotted to visualize the physical situations of the motion. 

Introduction 

The attention of many research workers in fluid mechanics has been diverted to the study 
of the influence of dust particles on the motion of fluids in the past few years. A model 
equation describing the motion of such mixed system has been given by Saffman 1 . Based upon 
Saffman’s model several authors 2,5 have investigated a number of dusty gas flow problems in 
different situations. 

There is another class of problems of the flow of dusty visco-elastic fluids such as latex 
particles emulsion paints, reinforcing particles in polymer mdts and rock crystals in molten 
lava. The study of these problems and rheological aspects of such flows have not received 
much attention although this has some bearing on the problems of petroleum industry and 
chemicaljmgineering. Several authors 6 * 9 have contributed in this field. Mandal et al. 10 have 
considered unsteady flow of dusty visco-elastic (Kuvshiniski type) liquid between two 
oscillating plates. 

Motivation for studying the problem is likely to have some industrial and chemical 
engineering application on the problems of transport of solid particles suspended in visco- 
elastic fluids through channels. 

Formulation and Solution of the Problem 

In the present investigation, we consider the laminar flow of an unsteady visco-elastic 
Kuvshiniski type liquid 11 with uniform distributiqn of dust particles down an inclined plane of 
inclination 0 to the horizointal. The liquid is bounded by a parallel upper surface at a distance 
h from the plane. 

Let us choose the origin of coordinate system at the bottom of the inclined plane x-axis 
opposite to the direction of the flow and along the greatest slope of the plane and y-axis 
perpendicular to the plane. Since both the dust and liquid particles move along the greatest 
slope of the plane- and the flow is laminar the velocity of both liquid and dust particles can be 
defined by the following relations : 



224 


R.K.S. CHAUDHARY and K.K. SINGH 


= «2 = 0 , «3 = 0 l 

v l = v l(y*0> v 2 = °. V 3 = °J W 


where (m 1( u 2 , u 3 ) and (v lt v 2 , v 3 ) are the velocity components of liquid and dust particles, 
respectively. 


Following Saffman 1 , the equations of motion for the flow of dusty visco-elastic liquid 
(Kuvshiniski) type) are given by 


1 + a 


jn 

dt 


du 

a 


M-;( 


1+ a— 


+ X) 


-\2 

a u x 


KNn 


aw 

a/Jax ” fy 2 p 


l + a— |(v] - «j) - gsin0 


( 2 ) 


1 dp 

- _ + gcos0 = O 


1 

P 


dp 

d2 


= 0 


dvi K . v 

dt m v ’ 


(3) 

(4) 

(5) 


where a is elastic parameters of the liquid particle, \>(=p/p) the kinematic viscosity, p the 
pressure, K the Stoke’ s resistance coefficient, m the mass of dust particles, N 0 the number 
density of dust particles and t the time. 

The initial and boundary conditions are 

t ^ 0 «i = 0 = Vj 

/>0 u 1 = 0 = v 1 aty = 0 


We express the pressure p as 

\ 

p = -pg(x sin 0 + y cos 0) - jc p<j>(t) 
Using eqn. (6) in (2) - (5), we get 


( a ^ a«! . , a 2 «j kn q ( av 

l 1 ° aFj-aT = F(t) + V + + a ^ ~ M i) 


3f m K 1 iJ 


( 6 ) 

(7) 

( 8 ) 


where F(t) = <{>(/) + a<j>'(f) 



FLOW OF A DUSTY VISCO-ELASTIC LIQUID 


225 


Let us choose u, v and F(t) as 

«i (y>t) = u(y)e- x2 ‘ 
vj (y,t) = v(y)e~ x2 ‘ ■ 

F(t) = ce~ x2 ' 


where c and 1 are the teal constants. 

Substituting the Values of u u v x and F(t) in eqn. (7) and (8), we get 

d 2 f 2 • 

-4+a/=o 

dy 


K 


v = 


K-mX 


where 


/ as A u + d 

,2 


M = 


X 
v 




Mm 


ml 


and d - — 
t> 


Here the boundary conditions are : 
f = d when y = 0 

f = -UA^+d when y = h 

The solution of eqn. (9) subject to the boundary conditions (12) is given by 
/ = Bj cos Ay + q sin Ay 


where 


B x = d 

d(l-cosAA)-£M 2 
Cl sin Ah 

The velocities of the liquid and dust particles are expressed as 


( 9 ) 


( 10 ) 

( 11 ) 


( 12 ) 
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«i = 


Vj = 


-d( 1 - cos Ay) sin Aft + jd(l - cos Aft) - UA 2 J sin Ay 


A 2 sin A A 


e ^ 2 ‘ and 


K 


K - rrOc 


-d(l - cos Ay) sin Ah + ■ d{ 1 - cos Ah) - UA 2 j sin Ay 


A z sin Ah 


.-7ft 


(15) 


The discharge of the flux per second for unit width of the plane is 
M 

Q = j o %yu 1 dy 

q _ — 1 J 2Ahl -Adh sin Ah + sin 2 Ah -d(l- cos Ah) cos Ah+UA 2 cos Ah] 

2A 4 sin Ah L 1 J 

+ A 2 dh 2 sin Ah + sin 2 Ah + 2d{l - cos Ah) sin Ah-2UA 2 sin Aft] ( 16 ) 


The velocities of liquid and dust particles at height ft above the inclined plane are 
respectively given by 


V , = 


-d( 1 - cos Aft) + d{\ - cos Aft) - UA 
• 


,2 ,2 

-X ( -X t 

= -Ue 


(17) 


and 


K 


K-mX 


-d( 1 - cos Aft) + d( 1 - cos Aft) - £/A 

P 

Discussion 




KU -X 2 t 
2 ^ 

K-mk 


(18) 


Fig. 1 exhibits visco-elastic effects on the flow field. The elastic parameter a is 
favourable to the velocity for both liquid and dust particles. The velocity profiles are non- 
linear close to the plate there after they become almost linear. From Fig. 2 it is evident that 
near the plate the velocity profiles fend to be linear, of course dusty phase being slower than 
the liquid phase as is expected in the physical situation. 



02H 



Fig* 1 - Velocity profiles for the liquid and the dust particles, 
sin AH r = -^3 / 2, — liquid, — dust. 
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U, , V| -> 

Fig. 2- Velocity profiles for the liquid and the dust particles, A = 2, u,= 1, cos Ah = 1/2, 
sin AH = ^3 /2, —liquid, — dust 
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Abstract 

Two formulae involving series summation and the diagamma function are established with the 
help of a generalized Riemann-Liouville operator of fractional calculus. 

Introduction 

The Riemann-Liouville fractional integral operator of order a of a function/^*) is defined 
and represented as 

<,[/(*)] = ^ J *< \* ~ 0°~/(0 dt, Re(cc) > 0 

= ^[RZ?[f(xj\],0<*+nZl,n = 0,l,2,... (1) 

while the Weyl operator of order a of a function /(*) is given by 

»£»[/«] = ^ JT ( ' " M«)>0 

. = (-l)^[C[/(^)]]' 0<a+n<l,n = 0,1,2... (2) 


A generalization of the Riemann-Liouville fractional integral operator defined by (1) was 
studied recently by Saigo and Raina 1 and represented and defined as follows 

Let a > 0, P and t| be real numbers. 




(3) 


where denotes the familiar Gauss Hypergeometric series, viz.. 


( 4 ) 


For a < 0 
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where 0<a + n^l, n is a positive integer. 

The inverse operator of /“ f ,T1 is given by (see Saigo 2 ) 

The following particular cases are due to Saigo and Raina 1 . 

CVwi - -<r 1 <” rrn 

where E (/(x)] is the fractional integral operator 1 . The inverse operator fori? is given by 

Ki/wir'-so.a/w] 

The two functional relations established in this paper are the following : 

For X > max {0, -a, -a-p, -a-p-T}} , 

V ( a ) B 


(5) 


( 6 ) 

(7) 

( 8 ) 

(9) 


«J a+P, a+B. -11 . ,1 

^«T(a+p+X+n) ^ 2 l a * a+ P +X+ " J 


( 10 ) 


and for X> max {0, a), . 


Y ( tt )n F 

2u \ p+Vq+l 




{ a p)> * ■£ 


~ r?^ f x y 

= X; f - (f) {^+r)-W(-a+X+r)} 


=o rrfz>;i ri 


(ii) 
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Results Required 

(1) If K > max (0, J3— ri) — 1 , then (see Saigo and Raina 1 ) 

7 a,p,n = r(jjr+i)r(-p+n+Ar+i) 


r(-p + K + l)r(a + r\ + K + 1) ' 


(2) From Gradshteyn and Ryzhik 3 , if Re(«+a) > 0 and Re (p+X) > 0, 

J*(jc-0“ +Ot “ 1 / P+X " 1 4 a +P, " n ; 1-^1 dt 

_ a+P+X+n -1 r(P + X)r(a + n) r a+ p 

r(a + fS + X + n) > A J 

Proofs of the Formulae 

To establish the result given in (10), we start with the formula 1 

j- a,-p, a+rTX-l]. _ r(x)r(a+p+x+Ti) p+x-i 
°- x L J r(p+x)r(x+ii) 

valid for X > max (0, -a-p-T|), which with the help of (5), can he put as 

[ X (x-t)~ a+n ~ l 2 Fi[" _a “^ dt 

dx [r(ti-a)L^° v L a+fl xj 

_ r(X)r(a+p+x+ T i) p+x-i 
r(p+x)r(x+ti) 

provided that X > max (0, -a-P-T|) and 0 < -a+n<l, a positive integer. 

Differentiating both sides of (15) with respect to X, we have 


■» f v“+P 


dx^Tf^nj y*- f r a+B_1 ^ 


- a -p,-«-u + ». 1 _i t X-lj dt 
X 


= {log x + g(a,p,X,Ti)} K (a,p,X,Ti) x P+X-1 
which, in view of the definition of operator I in (5), gives 

IoT ^’ a+T1 \* X ~ l lo & x \ = 0°g* + g(ot,p,X,T|)} ^(a.p.X.qJx^^ -1 


for X > max [0, -a-p-q), where we assume that 
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g(a,p, X, q) = ^(X) + T(a + p + X + rj) - ¥(p + X) - ¥(X + q) 


AT(a,p,X,q): 


r(x)r(a+p+x+g) 


r(p + X)r(X + Tl) • 

Next, consider the following Volterra type integral equation 
£(* - 0“ _1 afi[ a+P ’ ^ > 1 " f] fU) dt = r(a)x a+p+x - 1 log . 
which, in view of(3), can be re-written as 

With the help of the inverse operator given by (6), (21) yields 


This yields the solution of the integral equation in (20). Also from (17), we have 

fi*) = {logx + g(a,p,X,q)} AT(a,p,X,q)jt |i+x - 1 (23) 

We next substitute the value of fix) from (23) in (21) to get 

We now evaluate the two fractional integrals occurring on the L.H.S. of(24) We have bv 
definition (3), y 

h = / 0 “i , ‘ n p +X “ 1 log *] = - tf- 1 2 f[ « + P- -n 1 _ i.y^-1 Jog t dt (25) 

and since logr = logx- Y 

"j nx (26) 

Ktm by tem (which is vaM due ,o 

, '-' ogx 4 r p - -m - 4 (27) 

ft — 1 ' ' *- 
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and on further using (12) and (13) in (27), we get 

7 * — l2§jL y x r(p + X)T(a + n) r a+ p ia+K> _ TL I 

K( a,p,X,ri) "»r(a) r(a+p + X+«) [a.«+P+*+» ' l \ (28) 

Again from (12), 

r _ ra.p.iir p+X-l) _ r(P + A,)r(X+tl) 5 i-l „ , „ 

2 1* J-F(X.)r(oc+p+x+-n)' (29) 

Putting the values of 7j and in (24) and then simplifying the resultant expression, we 
arrive at the following result given in (10) 

V p fa+P.a+«, -Ti.il 

" n T(a + p + X + n) ^ 2 L «.«+P+^+» » J 

~ r(X)r(a +P + X+T1) + 'J'(a+P + X + 'n)-'P(p + X)-'P(X + Ti)] (30) 

provided X > max (0, -p, -x\, -a-p-r^). 

To obtain the secdnd result given in eqn. (1 1) above, we start with the formula 4 


which, with the help of (5), can be put as 

d f 1 rx, i i 




-a+X-1 ft \ i 1 

F ( a p)’ x .X 

r(-a + X) p+l 9+1 [ (b q ),~ a+x ’ a J (32) 

where 0<-a + n<l,«isa positive integer. 

aiv^- ff %^ ntiatmg 1)0111 SidCS ° f ( ' 32 ' ) with res P ect t0 ic gives, in view of the definition 
given m p;, 


d n 

1 f 

dx n 

r(-a + n)Jc 


_ 7 -a,a.,T\+a X~\ _ ( a P ) X . 

“ 7 0.x * A /{;- logx 


T(-a+X) 


r(x) iogx m f 


{“p)' X .£ 
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fl{aj) (X), 

+£ (f{# + r)- V (-a + Ur)} 

r "° nfe) r (-a + X) r r! 

J-i 


(33) 


since term by term differentiating is valid due to absolute convergence of the involved series. 
We now again consider the following integral equation 

J*(x- O a ' 1 /(tM=r(oc)x x - 1 ^ 


( a p ) . £ 

.(MV 


log X 


(34) 


The L.H.S. of which, in terms of the /-operator, can be put as ’ a,T| [/(x)] and which on, 
inverting with the help of (6) gives the solution of (34) in the following form 

m=io a x M+a 


X x ~' A 


( a p ) . £ 

(M ’ a 


log* 


(35) 


Now on substituting from eqn. (33) and (26) in eqn. (35), we get the desired result 
mentioned in (11) above. The method given above is essentially due to Kalla 5 . 

Particular cases 


If, in eqn. (10), we take P = 0 and use the result 



r(c)r(c-o-6) 

r(c-a)r(c-6)’ 


Re(c-a-b)>0andc?*0, -1,-2,... 


it reduces to the following result due to Kalla and Ross 6 


r(M-) v 

r w ... 


r(q+n) 
n r(|i + /t) 


= \|/(|i)-\|/(p-a) 


(36) 


where p = a + X, +T(. 

In in eqn. (11), we take p = 2 and q = 1, it reduces to another known result due to Kalla 5 , 
which in turn is a generalization of many known results as indicated therein. It may be pointed 
out that (11) agrees with a result given by Nishimoto and Srivastava 7 . 

By specialising the various parameters appearing in our main results given in (10) and 
(11) above, a number of results can be derived. 
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Abstract 

Two fixed point theorems have been presented over a quasi-metric space. In consequence, an 

analogue of Caristi-fixed point theorem has been obtained. 

Introduction 

In this paper, fixed point theorems over a quasi-metric space have been proved for a 
mapping satisfying certain contraction condition. For quasi-metric topology and allied topics, 
specially in setting of bi-topological. spaces one must see Kelly 1 . In proving theorem 1.1, we 
have used the notion of orbital continuity of the mapping. Theorem 1.1 also leads to Theorem 
1.2 which is an analogue of Caristi-fixed point theorem 2 in a quasi-metric space. Necessity of 
conditions in Theorem 1.2 has been examined by giving suitable examples. Continuity of 
fixed points has also been discussed. 

Definitions 


We start with necessary definitions : 

Definition 1.1 : A quasi-metric p on a nonempty set X is a non-negative real valued 
function pQ onXxX satisfying 

(i) p (x, y) = 0 if and only if x = y(x, y e X) 

(ii) p(x, z) <p(x, y ) +p(y, z); (x, y, z e X). 

Let pQ be a quasi-metric on X, then in a natural way there arises another quasi-metric <?(,) 
on X called conjugate of p(,), defined by q(x, y) = p(y,x) : (x, y e X). Sets of the form 
{y : p(x, y) < r) (r > 0) as x e X give rise to topology x p over X. Similarly, q gives rise to a 
topology on X. The topology on (X, p) is defined by x p . 

Definition 1.2 : A sequence {x n } ofX is said to be p-cauchy if and only if for e > 0, 3 m 
such that 


p(x r , x s ) < £ whenever r>s>m 

A sequence [x n ] in X is said to p-converge to a limit t; in X if and only if p(£, xj -> 0 as 
n -> and in this case [x n ] is called p-convergent. 

Definition 1.3 : (X,p) is said to be complete if and only if every p-cauchy sequence in X 
is p-convergent in X. 
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Definition 1.4 : In a quasi-metric space (x, p), a mapping T : (X, p) — > itself is said to 
be orbital continuous atx 0 e X if for some ue X, p(U, x„.) — > 0 gives p ( T(u ), r(x n; » -> 0. 

where {x„ ; } is a sequence of {x„} such that x y - T(x Q ), x 2 = 7’(xi), ,x n = T( JCn l ), 

T is said to be orbitally continuous if T is so at all points of X. 

Theorems 

Theorem 1.1 : Let (X, p) be a complete quasi-metric space. Let T : X -h> X satisfy 

(i) p(tx, ty) < a [q(x,Tx) + q{y,Ty)] + Pp(x,y) 

where a > 0, p > 0 and Zee + p < 1 

(ii) T is p-orbitally continuous at some point x 0 of X. . 

Then there is a unique fixed point of T in X. 

Proof : Define [x n ) in X where x n = T* x 0 , x n+1 = Tx n , n = 0, 1,2, 

without loss of generality, assume x n * .; B+1 ; then 
Pfa.x i) = p(Tx h Tx 0 ) 

< a[q(x lt Tx x ) + q(x o, Tx 0 )] + p p(x h x 0 ) 

= ap(x 2 , Xj) + aq(x Q , xi) + p q(x 0 , xj) 

or p(x 2 ,Xi)<y^- q{x 0 , Xj) 

p{x 3 ,x 2 ) = p{Px 2 , Tx { ) 

< a[q(x 2 , Tx 2 ) + q{x h Tjcj)] + P p(x 2 , x,) 

= a p(x 3 ,x 2 ) + a q(x h x 2 ) + p q{x h x 2 ) 

or p(x 3 ,.x 2 ) < Q(*b x l) 


and proceeding this way we have 
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= r n q{x o,^) 


( 1 ) 


, a + 3 . 

where r = — < 1. 

1-a 

••• p(xn^x n ) < p(x n+p , x„ +p _i)+ +p[x n + 1> x n) < [r n+p ~ l +......+r*\q(xQ, jq) 

r » 

< - — q(x 0 ,Xi)->0 as n-»°° 

which implies that [x n ] is a p-cauchy sequence in X. By completeness of X, there exists 
U € X such that p-lim x„ = u. By orbital continuity of T, 

p- lim T x„ = Tu i.e. p-lim x n+1 = Tu, i.e. u = Tu 

which shows that u is a fixed point of T. To prove uniqueness, let v e X be another fixed 
point of T. 

Thus, p(u,v) = p(Tu,Tv)<a[q(u,Tu) + q(v,Tvj\+$p(u,-v) 
i.e. p(w, v) < p(«, v); a contradiction. Hence u = v. 


Remark : In case p is a metric = d over X, we observe that condition (ii) of Theorem 1.1 
is redundant Because, from (1) we have 

£ d[t n+x x, r"x)<oo 

n~0 

and, hence Theorem 2 of Bollenbacher and Hicks 3 applies to ensure existence of a lower semi 
continuous (l.s.c.) function <p over X to satisfy 

d(7jc,x)<<p(x)~tp(r(x)), for xeX. 

Thus we arrive at- Caristi-fixed point theorem 2 , an application of which ends with fixed point 
of T: 

However, over a quasi-metric space we now prove Theorem 1.2 which is an analogue of 
Caristi-fixed point theorem. 

Theorem 1.2 : Let (X, p) be a complete quasi-metric space and let T be a mapping of X 
into itself. Let (p : X — » satisfy 

p{Tx, x) < <p(x) - <p(r(x)), xeX 

If T is orbitally continuous at x = x 0 then p-lim T n x = u exists and Tu = u. 


( 2 ) 
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Proof : Define {*„} in X where x n = T n x 0 , i.e. x„+i - Tx a , n = 0, 1,2, , without 

loss of generality let, x n * x n +i. Then 

p(x n + 1, *n) = Ppx x n ) < (p(x„) - <p(x„+i) 

£p(* v+1 , x v ) < (pfa) - < p(jc v +i) ^ <p(*i)- 

V=1 

Hence, 

oo 

*v) < °°- 

V— 1 

This implies that [x n ] is a p- cauchy sequence in X. By completeness of X, there exists u e X 
such that p-lim x n = u. By orbital continuity of T at x 0 , we have p-lim Tx n = u. 

i.e. p-lim x n+ i -Tu 
u = Tu 


which shows that u is a fixed point of T. 

We now examine necessity of conditions of Theorem 1.2. 

Example 1.1 beow shows that one cannot dispense with condition (2) of Theorem 1.2. 

Example 1.1 : Let us consider that complete metric space X = {0} u [1, °°] with usual 
metric of reals; and let T : X -» X be defined by 

T(x) = 0, 1 <* < oo and 7(0) = 1. 

T has no fixed point inspite of it being orbitally continuous at x = 0. However, we show 
that there does not exist such a function q> as to satisfy condition (2) of Theorem 1.2. Because, 
assuming the contrary, we have 

d(T0, 0) = if (1,0) = 1 < <p(0) - <p(l). 


cf(Tl,l) = d(0,l) = lS(p(l)-(p(0). 

This is impossible. 

On the other hand one cannot delete condition of orbital continuity of T from Theorem 
1.2 either. Example 1.2 below supports this contention. 


Example 1.2 : Consider the complete metric space X = 
reals and let T : X — » X be defined by 



u {0} with usual metric of 
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r|T| = -i- „>1 and 2X0) -1. 
n+1 

Here T has no fixed point. It is a routine exercise to see that T is nowhere orbitally 
continuous over X. 

We now construct a function <p : X -»R + to satisfy condition (2) of Theorem 1.2. 

Take <p(0) = 3 and (pf — J = — for n> 1. 

\nj n 


Here, we have 




Further d(0,r(0)) = l = <p(0)-cp(r(0)) 

and d(i,r(i))=i<i=<p(i)-q>(ni)) 


Thus (p satisfies condition (2) of Theorem 1.2. 

Continuity of Fixed Points 

We give Theorems 2.1 and 2.2 below in respect of continuity of fixed points in the 
setting of quasi-metric space. Their proofs run close to those of analogous theorems involving 
a restricted class of mappings in metric spaces as proved by Ray 4 . We include proof of 
Theorem 2.1 only to show the difference; and proof of Theorem 2.2 is left out by making 
reference to Ray 4 . 

Theorem 2.1 : Let (X,p) be a complete quasi-metric space, let T,- : X -*X be a mapping 
with a fixed point U; for i = 0, 1,2, and let {7}}! converge uniformly to T Q . If 

p ( t q ( x ), r 0 (y)) ^ a[F(x,r 0 (^)) + ^(y.Tb^))] +P p(*.>0 
for x, y e X and a is positive and 0 < p < 1; Then { p-converges to u 0 . 


Proof ; Let e > 0 be given. By uniform convergence of (T ; } to T 0 , there is an index N 
such that 



for all x e X and i 2: N. Then 
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p{uo,Ui) = p(T 0 (u 0 ), 

5 p{ T o{ u o) T o( u d) + p{ T o( u i)’ T i( u i)) 

■ S a[p(«o,r 0 (« 0 )) + q{ui, r 0 («i))] + P p(«o.«i) + p{To(ui)> T t (“i)) 
< a p(r 0 (u/), 7}(w,)) + pP(«o,«/) + p(T’ 0 (u 1 -), 2}(« f )) 


or p(u Q ,U ;) : 


1 + a 


p(3o(«/).^(“i)) 


.1-P. 

i.e. p( u o> u i) < 6 whenever i > TV. 


<e 


1 + a 

,1-P.J 


1-P 

1 + a 


Theorem 2.2 : Let (X, /? 0 ) be a complete quasi-metric space and let (p n ) be a sequence of 
quasi-metrics on X converging uniformly to (p 0 )- Let {L,,} be a sequence of mappings 
converging to p 0 -pointwise to a mapping T 0 with fixed point u 0 , and let each T n having fixed 
point u„ satisfy 


Pn{T n x , T n y) < a [p n (x, T n x) + q n (y, T n y)] + p p n (x, y), 


wh^re x, y e X and a is positive and 0 < p < 1. 

Then { u n ) p 0 -converges to Uq. The proof is similar to that of Theorem 3.2 of Ray 4 and is 
left out. 

Theorems 2.1 and 2.2 at once generalise Theorems 3.1 and 3.2 of Ray 4 respectively and 
we remark that uniform convergence of {p n } in Theorem 2.2 cannot be relaxed by pointwise 
convergence. This is evident from the example constructed by Ray 4 in the" setting of metric 
spaces. Also one cannot weaken uniform convergence of (T n ) to pointwise convergence in 
Theorem 2.1. This is supported by Example given below : 


Example : Take 


X = 



U{1}, 


with usual metric of reals. Then X is a complete metric space. 
Define T n : X — > X, as follows : 


7n(0) = 7n(l) — 1, for all n ; and 



1 . f . 

— , if i = n 
2 


= 1, if i * n. 


Then [T n ] converges pointwise to T 0 over X where T 0 : X -+ X, is given by T 0 (x) = 1 for 
all x e X. Here T 0 satisfies requirement of Theorem 2.1 with fixed point = 1. Also T n has 
fixed point = 1/2"; but 1/2" does not converge to fixed point of T 0 . 
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Abstract 

The photoacoustic spectra (PAS) of solid samples of l-chloro-2, 4-dinitro, 2,5-dibromo and 3,4- 
dichloro nitrobenzenes have been recorded and analysed in the light of their solution absorption 
spectra. The PAS of all the samples show three band systems. The longest wavelength band 
system HI, attributed to n-it* transitions, has been found to develop sharply in PAS whereas it 
developed weakly in their U.V. solution phase spectra. 

Introduction 

The photoacoustic spectroscopy is a technique to study spectra of solids and semi- 
solids 1 ' 9 . In this technique, heat is generated in a sample through non-radiative processes 
following the periodic illumination and consequent absorption of light by the sample. This 
technique measures the amount of absorbed optical energy de-excited through non-radiative 
processes and, therefore, the PAS corresponds qualitatively to optical spectrum and this has 
been used in : dentification 10 of the substances. In this work, we report PAS of l-chloro-2, 4- 
dinitro, 2,5-dibromo and 3,4-dichloro nitrobenzenes along with their solution absorption 
spectrum. Analysis of the PAS has been proposed in the light of their solution absorption 
spectrum, because a comparison of PAS with their solution absorption spectra suggests a 
possible method of identifying different electronic band systems in a PAS 11 . 

Materials and Methods 

The samples investigated viz . l-chloro-2, 4-dinitro, 2, 5-dibromo and 3,4-dichloro 
nitrobenzenes were obtained from Koch-Light Laboratories, England. The PAS of all the three 
solid samples were recorded using Princeton applied research model-600 1 spectrophotometer 
with 1 kW xenon arc lamp. The modulation frequency of 40 Hz was used with 
monochromator slit of slit width 2 mm in frequency region 200-500 nm. The scan rate was 
set at 200 nm/min. The electronic absorption spectra of these samples in methanol solution 
were recorded using (i) Carey- 14 spectrophotometer in region 200-800 nm and (ii) Cary-2300 
spectrophotometer in region 200-300 nm. The prominent features in photoacoustic and 
solution spectra are given in Table 1. The electronic origin of band systems I, II and III in 
PAS and the optical electronic spectra of the molecules are given in Table 2. 

The PAS of l-chloro-2, 4-dinitro, 2,5-dibromo and 3,4-dichloro nitrobenzenes show three 
regions of absorption. Each absorption region comprises of a main peak (most intense peak 
of the region) and a number of secondary absorption peaks of varying intensities. 
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Table 1 - Wavelength and prominent features in PAS and optical absorption electronic spectra of 1-chloro- 
2, 4-dinitro, 2,5-dibromo and 3,4-dichloronitrobenzene .. 


Transition 

PAS 

Electronic spectrum 




(A) 

(B) 

1 

2 

3 



l-chloro-2, 4-dinitrobenzene : 
Band system I 



nm 

int. 

nm 

int. 

nm 

int. 

Spectral region : 

200 - 224 

nm 

200 - 222 

nm 

200 - 225 

nm 

Most intense peak : 

220 

ms 

210 ' 

s 

205 

vs 

Absoiption peaks : 

205 

m 

210 

s 

206 

vs 


- 


- 


208 

s 


. 


- 


212 

s 


214 

m 

- 


215 

s 


220 

ms 

- 


218 

s 




- 


222 

s 

Band system II 







Spectral region : 

230 - 280 

nm 

230 - 280 

nm 

230 - 280 

nm 

Most intense peak : 

250 

s 

240 

s 

245 

vs 

Absorption peaks : 

231 

m 

- 


232 

s 


235 

m 

- 


236 

s 




- 


238 

s 


243 

ms 

- 


242 

s 


250 

s 

240 

s 

245' 

vs 


258 

ms 

260 

ms 

258 

s 


264 

ms 

- 


265 

s 


- 


- 


268 

s 


272 

m 

- 


273 

s 


276 

m 

- 


275 

s 

Band system HI 







Spectral region : 

335 - 380 

nm 

330 - 380 

nm 

X 

X 

Most intense peak : 


342 

ms 

350 

w 


Absorption peaks : 

360 

m 

- 





370 

m 

- 





(Contd.) 
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Table 1 (contd.) 

1 2 3 


2 ,5 -dibromonitrobenzene 
Band system I 


Spectral region : 

200 - 227 

nm 

200 - 225 

nm 

200 - 220 

nm 

Most intense peak : 

225 

ms 

220 

s 

210 

s 

Absorption peaks : 

205 

m 

- 


206 

m 


210 

m 

- 


208 

m 


- 


220 

s 

210 

s 


- 


- 


212 

ms 


215 

m 

- 


214 

ms 


225 

ms 

- 




Band system II 





; 


Spectral region : 

230 - 275 

nm 

225 - 280 

nm 

220 - 275 

nm 

Most intense peak : 

265 

vs 

230 

s 

232 

vs 

Absorption peaks : 

232 

s 

230 

vs 

232 

vs 


235 

s 

- 


235 

s 


240 

s 

— 


238 

s 


250 

s 

- 


- 



255 


- 


- 



• 265 

vs 

- 


- 



270 

s 

- 


- 


Band system IH 







Spectral region : 

280 - 370 

nm 

280 - 380 

nm 

- 


Most intense peak : 

300 

vs 

310 

w 

- 


Absorption peaks : 

285 

s 

- 


- 



300 

vs 

310 


- 



325 

s 

- 


- 



345 

s 

- 


- 



360 

s 

- 


- 


3, 4-dichloronitrobeniene 







Band system I 







Spectral region : 

200 - 225 

nm 

200 - 240 

nm 

200 - 240 

nm 

Most intense peak : 

214 

vs 

212 

vs 

212 

s 

Absorption peaks : 



- 


202 

m 


205 

m 

- 


204 

m 




- 


207 

ms 


208 

m 

- 


210 

ms 




212 

vs 

212 

s 


(Contd.) 
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Table 1 (contd.) 


1 

? 



3 




214 

ms 



218 

ms 


220 

m 

- 


•222 

ms 

Band system II 







Spectral region : 

235 - 290 

nm 

245 - 290 

nm 

245 - 290 

nm 

Most intense peak : 

250 

s 

270 

s 

’ 267 

ms 

Absorption peaks : 

248 

ms 

- 


247 

w 


250 

s 

- 


252 

w 


- 


, - 


255 

m 


- 


- 


257 

m 


- 


- 


262 

m 


265 

ms 

270 

s 

267 

ms 


- 


- 


277 

m 


277 

ms 

- 


- 


Band system HI 







Spectral region : 

300 - 360 

nm 

300 - 360 

nm 

- 


Most intense peak : 

328 

ms 

350 

m 

- 


Absorption peaks : 

312 

m 

- 


- 



320 

m 

- 


- 



328 

ms 

- 


- 



335 

m 

- 


- 



345 

m 

350 

m 

- 



355 

m 

- 


- 



Table 2 - 

Electronic origins of bands of l-chloro-2, 4-dinitro 2,5-dibromo and 3, 

4~dichloronitrobenzene 

Transition 

l-chloro-2, 4-dinitrobenzene 

2,5-dibromonitrobenzene 

3 ,4-dichloronitrobenzene 

PAS 

nm 

Electronic 

(A) 

nm 

Spectrum 

(B) 

nm 

PAS 

nm 

Electronic 

(A) 

nm 

Spectrum 

(B) 

nm 

PAS 

nm 

Electronic 
' (A) 
nm 

Spectrum 

<B) 

nm 

System I 

224 

222 

225 

227 

225 

218 

225 

240 

240 

System II 

280 

280 

280 

275 

280 

280 

290 

290 

290 

System HI 

380 

380 

- 

370 

380 

- 

360 


- 


Notations : int. = intensity, nm - nanometer, (A) = solution absorption electronic spectrum recorded on 
Cary-14 spectrophotometer, (B) = Elecronic spectrum recorded on Cary - 2300 spectrophotometer. 








Photo acoustic spectra. 








mo rmro benzene 
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The PAS of l-chloro-2, 4-dinitrobenzene (Fig. 1) shows that the distinction between 
band systems I and II and that between band systems II and HI is clearcut. The PAS signals of 
band system II are stronger as compared to those of systems I and HI. 

In PAS of 2, 5-dibromonitrobenzene (Fig. 1), the distinction between regions II and III is 
much confusing. The bands of region II appear to merge into region III. There is no clear cut 
distinction between systems I and II and between II and III. The band systems II and III are 
found to develop much sharply compared to the bands of system I. 

The PAS of 3, 4-dichloronitrobenzene (Fig. 1) shows that the distinction between band 
systems I and II is clearcut but it is confusing between band systems II and III. The band 
system III has been found to develop much sharply as compared to the band systems I and II. 
The bands of system II appear to merge into region III. The PAS signals of band systems II 
and III are nearly equally intense. PAS signals of band system I are much weaker as compared 
to the bands of systems II and HI. 

The solution absorption spectrum (A) (Fig. 2 & 3) of all the three samples show three 
band systems in good agreement with the three regions of absorptions in their PAS. The 
U.V. solution spectrum (B) shows only first two band systems. The band system III is found 
absent since this spectrophotometer records spectrum in region 200-300 nm. The spectrum 
(B) exhibits discrete vibrational structure of the first two band systems showing a large 
number of intense absorption peaks. 

Analysis : l-chloro-2, 4-dinitrobenzene : The PAS shows that the band system I extends 
the spectral region between 200-224 nm. The main peak lying at 220 nm, the secondary 
absorption peaks have been measured at 205 and 214 nm with medium intensities. The 
solution spectra (A) and (B) show this region between 200-225 nm. The spectrum (A) shows 
only one strong peak at 210 nm whereas the spectrum (B) shows a large number of 
absorption peaks at 201, 206, 208, 212, 215, 218 and 222 nm with strong 'intensities. The 
band observed at 206 nm is the most intense peak of region I. This most intense peak of 
electronic spectrum (B) at 206 nm which appears at 210 nm in electronic spectrum (A) of the 
substance may be correlated with the main peak of its PAS at 220 nm. The main peak of 
PAS of this band system is found to shift to longer wavelength side of the most intense 
peaks in its electronic optical spectra (A) and (B). 

Band system II observed in PAS of this molecule in spectral region 230-280 nm shows 
main peak at 250 nm with strong intensity. The secondary absorption peaks of this region 
have been measured at 231, 235, 243, 258, 264, 272 and 276 nm. Both optical spectra (A) 
and (B) show this spectral region between 230-280 nm which is same as that observed in its 
PAS. The spectrum (A) shows two absorption peaks at 240 and 260 nm of which the former 
is more intense. The spectrum (B) shows absorption peaks at 232, 236, 238, 242, 258, 265, 
268, 273 and 275 nm. The absorption peak at 245 nm is the most intense peak of band 
system n in spectrum (B) which appears at 240 nm in spectrum (A) - may be correlated with 
the main peak of the band system observed in its PAS at 250 nm. The main peak of PAS has 
been observed to shift to longer wavelength side of the most intense peak of optical spectra 
(A) and (B). 

The loiigest wavelength band system III attributed to n-it* transitions has been observed 
in its PAS in spectral region 335-380 nm with main peak lying at 342 nm and secondary 
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absorption peaks at 360 and 370 nm. The optical spectrum (A) shows this system in spectral 
region between 330-380 nm with only one weak absorption maxima at 350 nm. The main 
peak in its PAS observed at 342 nm may be correlated with the absorption maxima observed 
at 350 nm in its optical spectrum (A). It is observed that the main peak in its PAS is shifted 
to lower wavelength side, contrary to the observations of band systems I and II where it 
shifted to longer wavelength side of the most intense peak of the optical spectra (A) and (B). 

Like other substituted benzene derivatives, the first two band systems I and II observed in 
PAS as well as in optical spectra (A) and (B) of this molecule are attributed to tc-tc* 
transitions of phenyl group corresponding to 2600 A transition — 1 A lg and 2000 A 
transition l B lu — x A lg of benzene. The band system III on longest wavelength side is observed 
by electronic transition localized to a great extent in- — N0 2 part of the molecule. This spectra 
may be explained by considering electrons in n-orbitals as well as in %* -orbitals giving rise 
to n-n* transitions. This electronic transition would correspond to an allowed A - A 
transition of substituted benzenes with C s symmetry. This molecule shows low vapour 
pressure at ordinary working condition hence it is very difficult to record its vapour phase 
absorption spectrum. Thus the PAS studies of n-n* transition electronic system of this 
compound may give useful results. 

The spectral features of PAS and optical electronic spectra of 2, 5-dibromo and 3, 4- 
dichloronitrobenzene, showing most intense peak and other absorption peaks of varying 
intensities, are summarised in Table 1. 

Identification of electronic band origin : The identification of electronic origin of bands in 
PAS is, however, a difficult problem. The electronic origins of each band system in PAS of 
l-chloro-2,4-dinitro, 2, 5-dibromo and 3,4-dichloronitrobenzenes alongwith the electronic 
origin of bands in their solution absorption spectra (A) and (B) are given in Table 2. 

The region where the slope of PAS contour changes on longer wavelength side of system 
I around 224, 280 nm in band system II and 380 nm in band system III may be identified as 
the electronic origins of bands of system I, II and III in PAS of l-chloro-2, 4-dinitrobenzene. 

Likewise, the electronic origins of bands of systems I, II and III in PAS of 2, 5-dibromo 
and 3,4-dichloro nitrobenzens have been identified rnd listed in Table 2. 

Results and Discussion 

The PAS and the solution spectra (A) of l-chloro-2, 4-dinitro, 2,5-dibromo and 3,4- 
dichloronitrobenzenes show three absorption regions in good agreement with the three specttal 
regions of their optical solution absorption spectra. Region HI is found absent in solution 
spectra (B) because this records spectrum in region 200-300 nm only. 

The non-radiative transitions of system I in PAS of all the molecules become less and 
less effective as compared to band systems II and III. The photoacoustic signal is directly 
related to the loss of excitation energy by non-radiative processes. Although it is difficult to 
draw conclusion about individual molecules, we notice regularity in the present group of 
molecules. It is found that the non-radiative transitions in region I becomes less and less 
effective compared to those in regions II and III, the order being l-chloro-2,4-dinitro<3,4- 
dichloro- < 2, 5-dibromo-nitrobenzenes. This trend may be correlated with the electro 
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negativity of the substituents. 

Conclusion : The band system III attributed to n-n* transitions is found to develop 
much sharply in their PAS compared to that in their U.V. solution phase spectra. The band 
system III in their PAS shows a large number of absorption bands showing strong 
intensities. Thus PAS studies of n-n* 'transition in l-chloro-2 4-dinitro, 2,5-dibromo and 3, 
4-dichloro-nitrobenzenes, for which it is difficult to record vapour absorption spectrum at 
ordinary working conditions, may give useful informations. 

The three regions of absorptions observed in PAS are found to be in good agreement with 
those of the spectral regions and other features of their solution absorption spectra. This 
shows that the PAS and the solution absorption spectra correspond qualitatively to each other. 
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Abstract 

MINDO/3 semi-empirical molecular orbital technique has been used to calculate bond order indices, 
valence indices, oxidation numbers, ionization potentials, dipole moments etc. of aniline and 
isomeric diamino benzenes. The computed bond orders, the valence indices and ionization 
potentials are found very close to their corresponding classical values and experimental 
observations respectively. The reactivities of different atoms in terms of Jug’s hyper-, normal- 
and sub-valences are discussed in terms of their affinities for covalent bond formation. - 

Introduction 

One of the main objectives of the molecular orbital studies on the electronic structure of 
molecules is to elucidate the nature of the classical chemical bond in terms of quantum 
chemical parameters. Indices like bond order, valence and atomic charges are such classical 
concepts which have been not only used to parameterize and distinguish different bonding 
situations but also they represent useful link between the quantum chemical results and 
classical concepts. Despite the availability of better computation facilities, the quantum 
chemical calculations at ab initio level are expensive 1 even for small and medium size 
molecules. Therefore, semi-empirical methods are still employed to get the structural and 
electronic aspects of organic molecules 2 . Recently modified neglect of differential overlap and 
final version of modified intermediate neglect of differential overlap (MINDO/3) semi- 
empirical methods of Dewar’s group 3 - 4 have proven to be alternate to the ab initio method s . 

The electronic structure of substituted benzenes are of basic importance for a deeper 
understanding of their reactivities and spectral properties. The complete neglect of differential 
overlap (CNDO) and intermediate neglect of differential overlap (INDO) semi-empirieal 
methods have been used to study the electronic structures and geometries of some substituted 
benzenes 6 - 7 . A qualitative prediction of molecular geometry, both at MNDO semi-empirical 
and ab initio STO-3G SCF formation, has been carried out for some fluoro benzenes 8 * 10 . 
Several papers on bond order concept in semi-empirical formalism have been published which 
refer to bond quantities 11 , bond order 12 - 13 , bond indices 14 ot bond overlap 15 but most of them 
actually mean a property other than the valence multiplicity. Furthermore, Jug in his 
subsequent papers 16 * 19 , distinguished the bond order from the bond index 20 * 22 . The former is 
based on the evaluation of eigen values by diagonalising the inter-atomic part of the density 
matrix whereas the later is a function of the square of the density matrix elements. In addition 
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to this, in the same publication, he also provided a definition of valence index which is ab 
initio analogy of Mayer’s definition. 

Since we have continued interest in the electronic structure of substituted benzenes, we 
reported our results on these indices for several benzenes derivatives involving - F, - Cl, 
- OH, - NH 2 and - CH 3 as substituents in ab initio SCF theory 8-10 ' 23-27 , using Mayer’s 28 ' 29 and 
Giambiagi’s 30 definitions for bond orders, valence indices and oxidation numbers, respectively. 
In above publications we have preferred Mayer’s method over Elliott and Richard’s method for 
the computation of bond orders from ab initio molecular wave functions based on Dean and 
Richard’s 31 charge partition technique as the later method needs further parameterization for 
bonds other than C-C. 

Another useful concept for describing and interpreting the extent of electron transfer to or 
from the vicinity of a given molecule in a molecular environment and during the course of 
chemical reaction is the oxidation number. Unfortunately, no reliable theory is yet available 
for direct evaluation of the oxidation number and only qualitative quantum chemical 
approaches based on population analysis 32 or electron density distribution 33-35 are available. 
Recently Iwata 36 obtained analytical expressions for the number and density of electrons in a 
sphere centered at an arbitrary point, provided that the wave-functions are of Gaussian type. 
Using this method Takano et al . 37 ' 38 have systematically studied a series of compounds to 
elucidate the concept of oxidation number but some ambiguities still exist. On the other hand 
Giambiagi et al. 30 have defined oxidation number in terms Qf bond order- indices and found it 
very useful for neutral species. In our earlier paper 26 we have applied Giambiagi’s definition 30 
to calculate oxidation numbers for some benzene derivatives at ab initio level. Seeing the cost 
and time factor involved in the above ab initio calculations and suitability of MINDO/3 semi- 
empirical method, it was thought worthwhile to evaluate above parameters using Jug’s and 
Giambiagi’s definitions with MINDO/3 method and to see whether or not such calculations 
could be helpful in electronic structure studies. 

The main purpose is therefore, to provide reliable and fairly accurate information on the 
electronic structures of some amino-benzenes. In the present paper, we have evaluated bond 
orders, valence indices, exidation numbers, dipole moments and ionization potentials with 
MINDO/3 method for some amino-benzenes. We have already reported 39 the above parameters 
with the same technique for some fluoro derivatives of benzene. The results are discussed in 
the light of available experimental and other theoretical findings. 

Computational Details 

Gopinathan and Jug 21 - 22 suggested that the bond order matrix is sub divided into blocks 
referring to the atoms Afi.C, L in the molecule as 



IV" 

P AB 

p ac ;.... 

i — 

ft. 

p = 

P BA 

P BB 

p BC 



1 ^ 

p 18 

p 1 ^ 



L J (l) 

where P AA refers to charge and hybridization on atom A and P M refers to the bond between 
atoms A and 5. 
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The diatomic portion of P matrix represented by 




0 

P BA 


P M 

0 


( 2 ) 


is now diagonalised to get the characteristic eigenvalues X,- corresponding to bond order 
orbitals (BOOs) and eigenvectors of this matrix. The eigen values are seen to appear in pairs. 
Gopinathan and Jug 21 ' 22 defined the bond order B AS between the atoms A and B and valence 
index V A of atom A in terms of eigenvalues X; as 

+vc 2 

b ab = HLhiAB) (3) 

i 

i -fve 2 

and V A = 'Z'S'XtiAB) (4) 

Bit A i 


respectively. 

A further diagonalization of the intra-atomic part P** of the matrix P reproduces the 
eigenvalues (orbital occupancies) referred to as natural hybrid orbitals (NHOs). These eigen- 
values can be used to define the valence index V A by the equation 


V A = x( 2 «a-<?a) (5) 

a 

a’s being number of atomic orbitals centered on A and q's are their respective occupancies. 

Lastly the same operation applied on the ‘complete’ off diagonal sub-matrix defined by 
equation of Gopinathan and Jug 21 yields the corresponding natural bond orbitals (NBOs) 
eigenvalues £,• and are used to define the valence index of atom A by the equation 


+ve o 


( 6 ) 


These eigenvalues also appear in pairs. The percentage “excess valence” of different atoms 
in the molecule has also been calculated in MINDO/3 method by the expression 


E A = {v r A -v^xmiv r A 


(7) 


where V A is the reference valency of atom A and V A is its valency in the molecule under 
consideration. Following Mayer’s suggestion 40 ' 41 we use the term “excess valence” instead of 
Gopinathan and Jug’s “free valence” in order to avoid the use of latter in two different 
meanings. The definition of normal-, hyper- and sub-valence are given in Gopinathan and 
Jug’s papers 21 - 22 . 
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Giambiagi et al. 20 defined the oxidation number E A of an atom by the relation 

= (\Qa\!Qa) 1*ab ( 8 ) 

A*B 

where Q A is the net charge on A and the sum is carried out over all the atoms with a polarity 
different from that of A. 

The MINDO/3 program is used as such to evaluate net charges, ionization potentials and 
dipole moments and extended to incorporate Jug’s and Giambiagi’s suggestions to calculate 
bond orders, valence numbers and oxidation numbers respectively. 

Each one of the molecule under consideration has been treated as planar with all the ring 
angles of 120°. The successive substituents ground state geometry has been taken from the 
Landolt-B ornstein* 2 table. The C-C, C-H, C-N and N-H bond lengths are 1.397, 1.084, 
1.420 and 1.040 A° respectively and Z. C-N-H is 126.5°. 

Results and Discussions 

Bond orders and valence indices : The calculated bond order and valence indices as defined 
by Gopinathan and Jug 21 - 22 for mono and di amino benzenes are shown in Fig. 1 which shows 
that the magnitudes of the above parameters are very close to the values required classically. It 
is interesting to note that the eigenvalues A,- in semi-empirical theory are obtained just by 
diagonalising the two centre parts of the bond order matrix whereas in ab initio SCF theory 
one has to use the modified density matrix (orthogonalised basis sets) and then diagonalising 
its two centre parts in order to get the value of A,-. The modified density matrix with 
orthogonalised basis sets can be obtained by applying Lowdin’s orthogonalisation technique 43 
on the coefficient matrix with unorthogonalised basis sets. Eqn. (3) is used to evaluate bond 
orders from these eigenvalues. These eigenvalues occur in pairs i.e. ± A ; satisfying the pairing 
theorem 21 - 22 . Thus the eigenvectors corresponding to the positive and negative eigenvalues 
refer to bonding and anti-bonding orbitals. Obviously the smallest positive eigenvalue 
corresponds to the interaction of two basis functions which remain virtually unaffected during 
molecular formation. 

Similar to our prior studies 39 the magnitude of the bond order of C-C bonds attached to 
the substituent bond is found to have smallest value than the magnitudes corresponding to 
other C-C bonds. However, the magnitude of the Ce-Cx bond in o-diamino benzene differs 
considerably from the Q-Ci bond because of the fact that the Q-Q bond is flanked by 
substituent -NH 2 and therefore, reducing its magnitude of bond order. Furthermore if we 
compare the magnitude of bond order of Q-C 2 bond for aniline and p-diaminobenzene, we find 
both of them almost equal (a difference of 0.009) which means that the Ci-C 2 bond is least 
influenced by the presence of another substituent -NH-, at para position. 

Unlike the C-H bonds, the magnitude of bond orders of N-H bonds in all the molecules 
is the same indicating that they are unaffected with the position of the substituents. 

Valence indices computed using bond order matrix from (i) sum of the squares of the 
positive eigenvalues of the BOOs of the appropriate diatomic blocks (eqn.4), (ii) occupancies 
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ANILINE 


0-DIAMLNO BENZENE 



Tu-D! AMINO BENZENE 


F-Dl AMINO BENZENE 


Fig. 1 - Bond order indices of different bonds and valence indices of different atoms shown along the side 
of bonds and on the site of atoms respectively (NBOs, NHOs and BOOs values of valence indices 
have been found coincident) 
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of the NHOs (eqn.5), (iii) sum of the squares of the positive eigenvalues of the NBOs (eqn. 6) 
lead to the same values as summarised in Fig. 1. It is obvious from Fig. 1 that the calculated 
valence index values are very close to those coming from classical concepts. It is to be 
mentioned that a and it separability is maintained as one of the NHOs remains an unchanged 
2p orbital as is seen from the fact that diagonalisation does not affect the occupancy of this 
orbital. Similar results have also been found in earlier studies 9 ' 39 . 

In an attempt to explore the further possibility of covalent bond formation Mayer 28 ' 29 ' 40 ' 41 
proposed an idea of ‘excess valence’ and ‘valence defect’. One can easily see (Table 1) that the 
“excess valence” for all the hyper valent nitrogen atoms and “valence defect” for rest of the 
sub-valent atoms are so insignificant that all the atoms of the molecule can be treated as 
unreactive. The definitions of hyper-, normal and subvalences are given in Gopinathan and 
Jug’s original papers 21 ' 22 . However, the term “free valence” used by them for “excess valence” 
or “valence defect” is unfortunate because the free valence for closed shell systems should be 
equal to zero. This difference in the reference and calculated valence indices can be attributed to 
the delocalization and polarization effects. It should be mentioned that we are dealing with the 
reactivity as a tendency for covalent bond formation, nucleophilic or electrophilic reactivity is 
thus excluded. 


Table 1 -Valency and covalent chemical reactivity of mono and di-substituted amino benzenes 



Aniline 

o-diamino benzene 

m-diamino benzene 

p-diamino benzene 

Atoms 

Free valence 

Atoms 

Free valence 

Atoms 

Free valence 

Atoms 

Free valence 




(%)* 


(%)* 


(%)a 

C, 

+3.13 

p, 

+2.83 

Cl 

+3.38 

C, 

+2.88 

Ci 

+2.35 

C2 

+2.83 

C2 

+3.25 

C* 

+2.20 

Cj 

+2.15 

C, 

+2.23 

Ps 

+3.38 

C3 

+2.20 

c 4 

+2.18 

Q 

+2.10 

C 4 

+2.83 

C 4 

+2.88 

C 5 

+2.15 

Cj 

+2.10 

C5 

+2.30 

C 5 

+2.20 

Ce 

+2.50 

Qs 

+2.23 

Q 

+2.83 

Ct 

+2.20 

n 7 

-3.17 

n 7 

-2.40 

n 7 

-3.43 

N 7 

-2.37 



N« 

-2.40 

n 9 

-3.43 

N10 

-2.37 


*The free valence (%) values of hydrogen atoms in different molecules, not listed in order to reduce the size 
of the table, are found to have either zero or small free valence values-. 

A positive and negative signs of free valence values are indicative of sub and hyper valence respectively. 

Oxidation numbers and net atomic charges : The calculated MINDO/3 oxidation numbers 
and net atomic charges for aniline and isomeric diamino benzenes are reported in Table 2. This 
definition of oxidation numbers 30 , appropriate for neutral species, ensures that oxidation 
numbers in a molecule add up to zero. From the table, it is seen that a small variation in 
electronic densities may switch the polarities and yield very different oxidation numbers. Thus 
it clearly means that in covalent bonding situations, assigning oxidation number are best 
avoided 44 . Since the positive and negative signs of the oxidation numbers are correlated with 
the signs of the net atomic charges and therefore unable to describe the subde balance in a and 
n effects. The assumption of sp 2 hybridisation on the carbon atoms implies a charge 

population described as 2/, 2 p\ y lp y , 2p\, The present calculation do show a similar 
distribution of charges except on the carbon atoms attached to the substituent (table not 
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given). The deficit of electrons in the 2s orbital of the carbon atom, in the molecule under 
study as compared to free carbon representing promotion from 2s to a 2p orbital is explained 
as being due to the energy gain in the molecular formation from sp 2 hybridised orbitals. 


Table 2 - Oxidation numbers and net atomic charges of the atoms* 


Aniline 

o-diamine benzene 

m -diamino benzene 

p-diamino benzene 

Atoms 

Oxidation 

Numbers 

Atoms 

Oxidation 

Numbers 

Atoms 

Oxidation 

Numbers 

Atoms 

Oxidation 

Numbers 

Cl 

+3.707(+224) 

Cl 

+2.548(+130) 

Cl 

+3.786(+271) 

Cl 

+3.721(+165) 

C* 

-3.804(-104) 

C 2 

+2.548(+130) 

C 2 

-3.741 (-212) 

C* 

-2.326(-049) 

C 3 

+3.830(+061) 

C 3 

-2.331 (-056) 

C 3 

+3.786(+271) 

C, 

-2.326(-049) 

C 4 

-2.807(-063) 

C 4 

-0.135(~006) 

C 4 

—3.773 (—1 68) 

C 4 

+3.721 (+165) 

C 5 

+3.880(+061) 

C 5 

-0.132(-006) 

C 5 

+3.857(+l 12) 

C 5 

-2.226(-G49) 

Q 

+3.S04( — 104) 

c* 

-2.331 (-056) 

Q 

-3.773(-108) 

c* 

-2.331 (-049) 

Nj 

-2. 9 68 (-224) 

Nr 

-2.968(-220) 

N 7 

-3. 065 (-223) 

Nj 

-2. 9 47 (-224) 

H 8 

+0.971 (+006) 

N„ 

-2.968(-220) 

H» 

+0,967(+018) 

H 8 

+0.967(+003) 

H, 

-0.966(-010) 

H, 

+0.973 (+004) 

N, 

-2. 979 (-223) 

Hj 

+0.967(+003) 

H 10 

-0.032(-000) 

H 10 

-0.018(-004) 

H 10 

+0.968(+011) 

N 10 

-2.947(-224) 

Hu 

-0.966(-010) 

Hu 

-0.018(-004) 

Hu 

-0.9 67 (-01 4) 

Hu 

+0.967(+003) 

H 12 

+0.971(+006) 

% 

+0.973(+004) 

% 

+0.968(+011) 

Hn 

+0.967 (+003) 

Hy 

+0.970(+078) 

H u 

+0.958(+077) 

Hu 

+1.056(+080) 

Hu 

+0.977(+075) 

H m 

+0.970(+078) 

H» 

+0.973(+040) 

«14 

+0.970(+078) 

H m 

+0.972(+075) 



His 

+0.958(+077) 

H 15 

+0.969(+080) 

Hu 

+0.972(+075) 



H w 

+0.973 (+040) • 

h 16 

+0.970(+078) 

H,* 

+0.972(+075) 


♦Values of net atomic charges are in the parenthesis and are represented in units of 10 3 electrons. 


Ionization potentials and dipole moments : The values of the above parameters for the 
molecules under study are reported in Table 3. The ionization potentials calculated from 
Koopmann”s Theorems 43 are very much consistent with the experiments 46 ' 47 for aniline. The 
symmetry of the HOMO (highest occupied molecular orbital), an indication of the symmetry 


Table 3 - Dipole moments and ionization potentials of mono and di-substituted amino benzenes 


Molecules 

Dipole moments 
(Debye) 

Ionization potential 
(in atomic units) 

Aniline 

0.98 

-0.289 K 


(0.94)* 

(-0.282) c 


(0.39) b 

(-0.295) d 


(1.53) e 


0 -Diamino benzene 

1.50 

-0.269 7t 

m-Diammo benzene 

1.00 

-0.281 7t 

p-Damino benzene 

0.00 

-0.263 K 


* Yadav, J.S. (1973) Ph. D. Thesis, Varanasi, b Singh, O.P. & Yadav, J.S. (Unpublished work), 

0 Ref. 46, d Ref. 47. e Ref. 48. 
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of the ionized species, is found to be of rc symmetry. This is in agreement with the usual 
assumption that the excitation and ionization mechanisms are mainly controlled by electrons 
in conjugated and aromatic molecules. However, the dipole moments differ considerably from 
the experimental 4 ® as well as other theoretical findings (Table 3). 

Conclusions 

Keeping the cost of calculation and computation time involved, it can be concluded that 
the MENDO/3 method may be used for elucidating such useful concepts like bond orders, 
valence indices, and parameters like ionization potentials to other set of molecules. However, 
the dipole moment could not be reproduced by this method and therefore, it is concluded that 
the parameterization incorporated for its calculation in MINDO/3 is not appropriate. On the 
other hand the AMI (Austin Model l) 49 model provides satisfactory results of dipole moment 
which are very close to experimental findings. The evaluation of dipole moments with AM 1 
model for such systems is in process and the results will be published in future 
communications. 
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Abstract 

The I.R. spectrum of thiostrepton-copper(II) complex (Ths-Cu) of 1 : 4 molar ratio indicated the 
involvement of some of the -NH, -C=N, -C=0 and -C-S groups pertaining to amide 
azomethine, ester carbonyl and thrazole rings respectively of the thiostrepton molecule in co- 
ordination to copper(It). The nature of one metal binding site of the antibiotic molecule was 
understood by constructing a model. The subnormal Bohr magneton value of the paramagnetic 
species of Ths-Cu as revealed by magnetic suceptibility measurement suggested an intense metal- 
metal interactions. Optical absorption study of Ths-Cu in different solvents indicated a distorted 
octahedral geometry of the complex with the various donor atoms of the ligand molecule 
occupying the equitorial plane. The structural requirements of the antibiotic for biological 
activity are discussed. 


Introduction 

Interaction of antimicrobial agents with metal ions may result in enhancing 1 , stabilizing 2 , 
or abolishing 3 their biological activity. Study of metal interactions of antimicrobial agents! 
therefore, is envisaged to provide valuable information about their structure-activity 
correlations especially when their metal co-ordinating sites are essential for biological activity. 

The molecular structure of thiostrepton (Ths.) (Fig. I) 4 , a peptide antibiotic produced by 
Streptomyces azureus 5 a soil isolate, reveals that it contains sites at which co-ordination to 
metal ions can occur. During the study of the interaction of Ths. with various bivalent metal 
cations in relation to biological activity, it was observed that the antibiotic was found to 
selectively bind to copper(II) in a stoichiometric ratio of 1 : 4 with the loss of its biological 
activity 4 . The present paper deals with the synthesis and characterization of Ths-Cu complex. 

Materials and Methods 

The following instruments were used; (a) Pye-Unicam SP3 - 300 ratio recording double 
beam spectrophotometer for recording the I.R. spectrum. Polystyrene (1601 cm" 1 ) was used for 
calibration; (b) Beckman DU-50 recording spectrophotometer for electronic absorption spectra; 
(c) Perkin Elmer atomic absorption spectrophotometer for estimating the metal content in 
Ths-Cu; (d) EE & G Princeton Applied Research Model 155 vibrating sample magnetometer 
for measuring the magnetic suceptibility of Ths-Cu at room temperature. 
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Fig. 1 - A schematic view of the molecular conformation of thiostrepton molecule. Deala - Dehydralanine; 
Q Quinaldio acid; Thstn - Thiostreptine; Thz - Thiazole ring; Pip - Piperidine ring. 


Nitrogen and chloride analyses were carried out using the standard microanalytical 
procedures. 

Thiostrepton standard (M r 1650) was a gift sample from the Squibb Institute, New Jersey, 
USA Al other chemicals were of Analar grade and used directly. 

Preparation of thiostrepton-copper(II) complex : The synthesis of Ths-Cu was carried out 
by mixing to a 4 ml 1,4 dioxane solution of Ths. (1 x 10 M), 1 ml of aqueous copper 
chloride solution containing excess copper(II) (6 x 10’ M) arid incubating the mixture at 37°C 
for 24 h. Later, to this solution, 200ml dioxane was added and stirred. A light green coloured 
product was precipitated which was collected by centrifugation. The precipitate was washed 
successively thrice with dioxane and finally once with ether and dried in vacuo over calcium 
chloride. Yield 78%. The purity of the product was checked by thin layer chromatography 
using various solvent systems. 


Results and Discussion 

Some properties of Ths-Cu as compared to that of Ths. are summarized in Table 1. These 
results indicate that the copper(II) complex of Ths. is an altogether a new molecule differing 
from native Ths. in many respects. Analytical data shows Ths. to copper (II) and Ths. to 
chloride molar ratios as 1 : 4 and 1 : 8 respectively. Therefore, the Ths-Cu can be represented 
with the molecular formula as (CuJhsClg). Conductance data (1 x 10~ 3 M, DMF) showed that 
Ths-Cu is a non-electrolyte indicating that the eight chloride ions are in the co-ordination 
sphere. 
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Table 1 - Comparison of Ths-Cu with the thiostrepton standard. 


Property 

Thiostrepton 

Ths-Cu 

Mol. Wt. 

1650 

2176 

Decomp. Temp. (°C) 

246-258 

170-180 

Solubility 

Sol. in DMSO, DMF, pyridine, 

Sol. in DMSO, DMF, pyridine, acetic 

dioxane, CHC1 3 , acetic acid and higher 

acid, methanol and in ImM EDTA (pH 


1 chain alcohols. Insol. in lower 
alcohols. 

7). Insol in higher alcohols. 

Biological activity 
Absorption spectrum 

Active on Gram +ve bacteria 

Ineffective on bacteria 

UV (20% DMSO) 

255nm peak and a shoulder at 280nm 

260nm peak and shoulders at 280 & 
360 nm. 

Visible (pyridine) 

No absorption 

485nm shoulder and 600-800nm 
peaks. 

Specific rotation 
[a]23°D, (DMSO) 

(-) 54“ 

(-) 16° 

Mobility on Tic 

Higher Rf value in non-polar solvent 

Lower Rf value in non-polar solvent 

systems & lower Rf value in polar 

systems & higher Rf value in polar 


solvent systems 

solvent systems 

I R. spectrum 

The spectra of both the samples are different and not superimposable 

Elemental analysis 

N 

16.1% 

Expected - 12.1% 

Observed - 11.9% 

Cu 

- 

4 atoms/molecule thiostrepton 

Cl 

- 

8 atoms/molecule thiostrepton 

Table 2 - Assignments of some of the principal IR absorption bands of thiostrepton. 

Wave numbers 

Approximate frequency description 

(cnr 1 ) 



3300 

Amide A band of -NH & -OH & NH 2 str. 

3000 

Aromatic -C-N str. 


1720 

Ester carbonyl str. 


1680-1600 

Amide 1 of -0=0 str. -NH def. -C-N str. and -C-C bending 

1575-1480 

Amide II of -C-N str. -NH def. 


1200 

-C-0 str. of Sec/Ter. alcohol 


800 

-NH def. 


750 

Amide IV of 0=C-N def. 


587-708 

-C-S str. 


200 

Amide VII tor. 



Str. - stretching frequency; def. - defonnation/bendng frequency; tor. - torque. 


The IR spectrum ofThs. and Ths-Cu were recorded in nujol and in KBr media (Fig. 2). 
Since nujol has intense absorption signals between 2000-3000 cm -1 and also at 1300-1500 
cm -1 frequencies 7 , the spectra of the samples in these regions was read from the KBr media. 

The I.R. spectrum of the native Ths. molecule shows intense absorption bands (cm" 1 ) at 
3300-3500, 2920-2970, 1720, 1600-1680, 1480-1575, 575-710 (Fig. 2) that can be assigned 
to -OH, -NH str; -C=N str; ester -C=0 str; amide I of -C=0 str; and amide II of -C-N str; 
and-C-S str. respectively (Table 2)*. The high intensity of some of these bands is attributed 
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to multiple number of the corresponding functional groups present in the Ths. molecule. 

From the comparison of the I.R. spectrum of metal complex with that of the native 
molecule (Fig. 2), it is evident that the absorption bands pertaining to amide, azomethine, 
ester and amide carbonyl and -C-S groupings of the ligand are utilized in co-ordination to 
copper(II). Since many of these ligand groups are available in the Ths. molecule, the binding 
of the copper (II) to some of them resulted in broadening or reduced intensity of the 
corresponding absorption bands. For e.g. broadening of the absorption band at 3300-3500 cm- 1 
in Ths-Cu can be attributed to some metal bound -NH groups and unbound -OH groups*. The 
involvement of peptide nitrogens of Ths. in co-ordinating to copper(II) was also indicated in 
potentiometric titration study of Ths. in presence of the metal 6 . There is only one ester-C=0 
group in Ths. molecule (1720 cm- 1 ). Since the I.R. spectrum of Ths-Cu shows a downward 
shift of this group from 1720 to 1710 cm- 1 , it indicates the involvement of the ester -C=0 
group in co-ordination. Similarly, the absorption band of Ths. pertaining to -C-S str. 
(575 cm- 1 ) is much reduced in intensity in Ths-Cu indicating the involvement of the sulfur 
atoms of the thiazole ring (Thz) of the antibiotic in co-ordination to the metal ion. The -OH 
groups of the antibiotic are not expected to interact with the metal since, no stable chelates 
could be formed involving these groupings. 

Some of the absorption frequencies in the region 600-900 cm- 1 of Ths. are observed to be 
shifted to lower frequencies in the Ths-Cu, like; 880 cm- 1 absorption band of Ths. appears to 
be shifted to 860 cm- 1 in Ths-Cu. And the 200 cm- 1 absorption band of Ths. pertaining to 
amide VII tor. is reduced in intensity in Ths-Cu. Likewise, the region between 200-500 cm- 1 
of Ths-Cu presents a totally different picture as compared to the native Ths. molecule with the 
appearance of some new structures that can be attributed to M — > N, M — > O and M S 
bonding 8 - 9 . 


From these data it is hence deduced that the nitrogen atoms of the amide and azomethine 
groups, some of the oxygen atoms of the amide and ester carbonyl groups, and sulfur atoms of 
the thiazole ring of Ths. molecule are involved in co-ordinating to the metal ion. 

The solution conformation of Ths. molecule (Fig. 1) reveals that there are several -NH 
groups with the other co-ordinating groups i.e. -C=0, -C=N, -C-S at favourable position to 
achieve efficient chelation with copper(II). Especially, the involvement of nitrogen atom of 
piperidine (pip) and sulfur atom of Thz(4), -NH and NH 2 groups of dehydroalganine (deala) (2) 
and (3) would result in a stable five membered chelate with the possible arrangement shown in 
Fig. 3. A molecular model of this region of Ths. molecule with the bound copper(II) was 
constructed. This approach showed that, with the rotation of the single bond between pip and 
Thz (4) (shows by an arrow mark in Fig. 3), it is also possible to involve the nitrogen atom 
instead of the sulfur atom of the Thz(4) in co-ordination to the metal ion. Since the entire 
‘tail’ portion of the antibiotic molecule (comprising of the two deala groups, and Thz(4)) are 
bound to the pip ring), it is possible that, in solution, there could exist a dynamic 
equilibrium between the nitrogen bound Thz (4) and sulfur bound Fhz (4) complexes. 

It is interesting to observe from Fig. 1 that, binding of one copper(II) ion at the ‘tail’ 
portion of Ths. molecule does not drastically affect the overall conformation of the antibiotic 
molecule because, the rest of the portion of die molecule is unaffected. In such a case, binding 
°f a metal ion at this region would only slightly affect the biological activity of the 
antibiotic. Our observations on biological activity of Ths. in presence of copper (II) at 1:1 
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Fig. 3 - Schematic representation of one possible region of thiostrepton molecule with the bound 
copper(H). R - rest of the portion of the thiostrepton molecule. For purpose of clarity, chlorides in 
the axial position of the metal ion are not shown in the diagram. 



Fig. 4 - Absorption spectra 
pyridine (...) 


of Ths-Cu in the visible region in DMSO/DMF ( — ), methanol ( ) and in 
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molar ratio showed only 12% loss of activity of the antibiotic 6 . Based on this observation, it 
is tempting to assume that, at 1:1 molar ratio of Ths. to copper(II), this region may be the 
site that is initially involved in co-ordinating to the metal thereby forming a stable five 
membered chelate (Fig. 3). Such stable five membered structures involving the other three 
copper(II) ions at the other binding sites in the antibotic molecule could also be possible. 
Therefore, the binding of four copper (II) ions to Ths. molecule is evidently due to the 
availability of four different sites in the antibiotic where the co-ordinating atoms are present. 
All the four regions may not be, however, initially accessible for binding to the four 
copper(II) ions because of unfavourable stereochemical orientations of some of the ligand 
atoms in these regions. Possibly, the initial binding of one or two metal ions at any of the 
favoured regions of the molecule would induce electronic perturbations in the entire ring 
structure of Ths. 10 which would favour the binding by the rest of the copper(II) ions. The 
binding of all the four copper(II) ions to Ths. molecule would, therefore, greatly affect the 
overall conformation of the antibiotic molecule rendering it biologically ineffective. Alteration 
in the I.R. spectrum of Ths-Cu in the fingerprint region (1300-900 cm -1 ) compared to that of 
Ths. indicates a change in the conformation of the antibiotic molecule as a result of binding of 
the four copper(II) ions. 

The electronic absorption spectrum of Ths-Cu compared to that of the native Ths. 
molecule in the U.V. region indicated shifts in both pi — ¥ pi* and n -h> pi* transitions of the 
free Ths. 6 Appearance of a new band at 360 nm in the Ths-Cu molecule suggested a 
bathochromic shift in the transitions jthat are attributed to metal co-ordination. 

i 

The absorption of Ths-Cu in the visible region recorded in different solvents like; 
methanol, dimethylsulfoxide (DMSO) dimethylformamide (DMF) and pyridine shows that, as 
the basicity of the solvent increases from methanol to pyridine, the broadening as well as the 
intensity of the band between 600-800nm increases (Fig. 4). In pyridine, the spectrum shows 
a broad intense multiple band in thd region of 600-800 nm which is resolved into three 
regions, viz. 700, 750 and 780nm being centered at 750nm. These can be attributed to 
2 Blg -> 2 Alg , 2 Blg -» 2b 2j , and 2 Bjg -> 2 Bg transitions respectively that are consistent with a 
distorted octahedral geometry 1 1 . 

Magnetic sucptibility measurement ofThs-Cu indicated a sub-normal Bohr magneton 
value for a single copper(II) ion (0.533). This indicates the existance of a strong metal to 
metal interaction between the copper(II) ions. The Ths-Cu molecule with the four copper(II) 
ions can show subnormal magnetic moment due to the two possible antiferromagnetic 
interactions; (i) the complex molecule can fold in such a way that the four metal ions within 
the molecule can come close to undergo metal to metal interaction trough chloride bridging; or 
(ii) if the complex molecules are assumed to stack one above the other, a copper(II) ion in one 
complex molecule in a distorted octahedral geometry can under spin-spin exchange with 
similar copper(II) ion in another complex molecule through the chloride bridging 12 . 

It was due to these interactions between the four copper(II) ions present well within the 
Ths. molecule that resulted in a single broad unresolved line in the E.S.R. spectrum of Ths- 
Cu in both DMSO and DMF even at liquid nitrogen temperature though, in pyridine, the 
spectrum showed good resolution with superhyperfine splittings (results not presented here). A 
detailed study of E.S.R. investigations are communicated to Current Science, Bangalore. 
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These results suggest that the four copper(II) ions in the Ths. molecule have a distorted 
octahedral geometry with the various N/ O/S donors of the antibiotic molecule occupying the 
equitorial square plane, and with the chloride ions at longer distances in the axial plane. 
However, when the complex was dissolved in pyridine, the chlorides that are weakly 
interacting with the metal ion are being displaced by the solvent molecules. 

In conclusion, it should be emphasized that, it is also possible that the suggested metal 
co-ordinating groups of the Ths. may not be directly involved in its biological activity. This 
is because, metal induced changes in the architecture of an antibiotic molecule at groupings 
other than those involved in metal binding may also render the molecule biologically 
ineffective 3 . Therefore, it is the changes induced in the overall conformation of the Ths. 
molecule as a result of binding of the four copper(II) ions that resulted in elimination of its 
biological activity. 
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Abstract 

Diethazine hydrochloride is proposed as a new reagent for the fapid spectrophotometric 
determination of cerium (TV), arsenic (03) and nitrite. The reagent forms red coloured species with 
cerium (TV) instantaneously at room temperature in phosphoric acid medium. The species has an 
absorption maximum at 516 nm with a sensitivity of 15 ng cm" 2 . The method has been used for 
the analysis of cerium alloys and for the determination of arsenic (D3) and nitrite. 

The colour reaction of diethazine hydrochloride (DH), 10 [2-(diethylamino)ethyl] 
phenothiazine hydrochloride with cerium (IV) lias been studied for the rapid spectrophotometric 
determination of cerium (IV), arsenic (III) and nitrite. The proposed method offers the 
advantages of simplicity, rapidity, sensitivity and stability without the need for heating or 
extraction. 

Stock solutions of cerium (IV) sulphate in 0.5 M sulphuric acid, sodium arsenite 
(prepared from arsenic trioxide) and sodium nitrite were prepared and standardised by literature 
method 1 . A 0.3% (wt/vol) solution of DH was prepared in doubly distilled water and stored in 
an amber bottle in a refrigerator. Solutions of diverse ions of suitable concentrations were 
prepared using analytical grade reagents. Beckman model DB spectrophotometer with matched 
1 cm quartz cell was used for absorbance measurements. 

Determination of cerium (IV) : An aliquot of the stock solution containing 12-500 pg of 
cerium (IV), 15 ml of 10M phosphoric acid and 2 ml of 0.3% DH solution were taken and the 
solution diluted to 25 ml with doubly distilled water. The solution was mixed well and the 
absorbance was measured at 516 nm against a corresponding reagent blank prepared in the 
same manner. The amount of cerium in the sample solution was deduced from the calibration 
curve. 

Determination of arsenic (III) : 1 ml of 5 M sulphuric acid, 0.5 ml of osmium (VIII) 
(1 pg) and 1 ml of 0.002M cerium (IV) solution were transferred to each of eleven 25 ml 
volumetric flasks. An aliquot of the stock solution containing 3.8-87.5 pg of arsenic (III) was 
added only to ten 25 ml flasks and mixed well. An eleventh 25 ml flask was used for a 
simultaneous blank determination. 15 ml of l'OM phosphoric acid and 2 ml of 0.3% DH 
solution were added to all the flasks, diluted to the mark and shaken well and the absorbances 
measured at 516 nm against a reagent blank containing no cerium (IV) and arsenic (III). The 
amount of arsenic (III) in a test solution was then deduced from the calibration graph 
constructed by plotting the concentration of arsenic (III) vs. the difference in the absorbance 
readings between the blank (in eleventh flask) and the sample. 

Determination of nitrite : 1 ml of 0.002M cerium (IV) sulphate and 2 ml of 2M sulphuric 
acid were taken in each of eleven 25 ml volumetric flasks. 2 ml of sodium nitrite solution 
containing 2.5-62.5 pg of nitrite was transferred into each of ten volumetric flasks and the 
eleventh flask was used for blank experiment. The solutions were well mixed, left for 5 min, 
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the colour was developed by adding phosphoric acid and DH solution, the absorbances 
measured and the amount of nitrite in the test solution was determined as described in the 
procedure for the determination of arsenic (III). 

Cerium (IV) readily oxidises DH to red coloured species in sulphuric or phosphoric acid 
medium. The maximum colour development takes place instantaneously at room temperature 
(28 ± 1°C) in 5.5-6.3M phosphoric acid solutions. The red coloured species assumed to be a 
radical cation 1 shows maximum absorption at 516 nm and the absorbance remains constant for 
about 45 min. Phosphoric acid (6M) was chosen because in the presence of sulphuric acid the 
reaction is less sensitive and less stable. The reagent blank and the metal ion under similar 
conditions do not absorb around this wavelength thus promoting excellent analytical 
conditions. A 24-fold molar excess of DH is necessary for maximum colour development. The 
absorbance values are insensitive to temperature range 5-55°C, but beyond 55°C, the 
absorbance gradually decreases. The red species is cationic in nature as revealed by ion- 
exchange experiments. There is no appreciable change in the absorbance if the order of addition 
of reactants is varied. Beer’s law is obeyed over the concentration range 0.5-20 ppm of Ce(IV). 
The optimum concentration range evaluated by Ringbom’s method is 1.2-18.0 ppm. The 
molar absorptivity and Sandell’s sensitivity are 9.7 x 10 3 dm 3 mol -1 cm- 1 and 15 ng cm -2 
respectively. The standard deviation calculated from six determinations in a solution 
containing 4 ppm of Ce(IV) is 0.01 and the relative error is less than 2%. 

It is seen that many diverse ions especially lanthanides do not interfere in the 
determination of Ce(IV). The major advantage of this method is that DH can be used as a 
selective reagent for the determination of Ce(IV) in presence of large quantities of other 
lanthanides in readily attainable oxidation states without the use of masking agents. The 
proposed method was successfully applied to the analysis of misch metal containing 50% Ce, 
25% La, 15% Nd, 5% Fe and 5% mixtures of Pr, Eu, Gd and Er. 

Arsenic (III) and nitrite were indirectly determined. Arsenic (III) was quantitatively 
oxidised instantaneously to arsenic (V) by a known excess of cerium (IV) sulphate (20-80%) 
in 0.5M sulphuric acid medium containing 1 pg of Os (VIII) which did not interfere under the 
experimental conditions. Nitrite was quantitatively oxidised to nitrate in 0.5M sulphuric acid 
in 5 min, by a known excess of cerium (IV) sulphate (20-80%). The unreacted Ce (IV) was 
then determined by the proposed method. In the absence of other reducing substances, the 
reduction in the absorbance of the red colour produced by the fixed amount of cerium (IV) was 
directly proportional to the amount of arsenic (III) and nitrite present. Ce(III), As(V) and 
nitrate formed in the reaction were colourless and did not interfere. 0. 15-3.5 ppm of As(III) and 

0. 1.2.5 ppm of nitrite could be determined. This method can be used for the determination of 
micro-quantities of other substances which are quantitatively oxidised by Ce(IV) in sulphuric 
acid medium to colourless non-interfering products. 
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Abstract 


The spectrophotometric method for determining cobalt(II) using diethyldithiocarbamate has been 
modified by which even quite small quantities of Co(II) can be determined. 


Sodium diethyldithiocarbamate (Na 2 L) is used for the determination of various metal 
ions 1 . Lacoste et al. 2 have reported determination of cobalt (II) using sodium 
diethyldithiocarbamate. A green coloured complex CoL^ (structure I) is formed which absorbs 
at 650 nm in CHC1 3 . 


Ccn 




C-N(C^H 5 )2 


2 


( 1 ) 


We have developed a calibration curve in the range 1-6 mg of Co(H) / ml. For developing 
the colour, 1 ml of a Co(II) solution in the above concentration range was treated with 2 ml of 
2% sodium diethyldithiocarbamate and 10 ml of the acetate buffer (C^COOH-CHjCOQNth) 
of pH 6.3. This aqueous solution was extracted with two lots of 20 ml of CHC1 3 and volume 
was made to 50 ml. It gave Xn, ai at 650 nm on a Beckman spectrophotometer model DU. The 
same procedure was repeated for other concentrations. 

Recently, when we attempted to estimate Co(II) in the concentration range less than 
0.3 mg/ml, it did not give any absorption when the above procedure was adopted. However, 
when instead of 50 ml, a total of 25 ml volume of the solution was made, it gave absorption, 
thus it was decided to modify the above method. 

In the- modified method, which is recommended as a graduate level experiment, a 
calibration curve in the range 0.1 to 1.9 mg Co(II)/ml was made. For this, 1 ml of the 
standard Co(II) solution containing cobalt (II) in the concentration range 0.1 to 1.9 mg/ml was 
treated with 10 ml of buffer and the aqueous solution was extracted using 10 ml, lots of 
CHC1 3 instead of 20 ml (Lacoste’ procedure) and volume made to 25 instead of 50 ml. The 
calibration curve is shown in Fig. 1 and one can see that quite small quantities of Co(II) can 
be determined. The was 640 instead of 650 nm and it may be due to the use of different 
instruments. In the modified method, Shimdazu-UV-visible-spectrophotometer was used. 
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Concentration (mg) > 

Fig. 1 - N, N-diethyldithiocarbamate method for the estimation of Co(II), concentration range 

0.1 - 1.9 mg, X. mftx = 640 nm. 

A standard Co(II) solution was made by EDTA method using xylenol orange as the 
indicator (end point-pink to yellow) 3 . 
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ABSTRACT 

In this note, we shall consider and unify a variety of problems on differential systems. For 
instance, the results cover global existence, boundedness and stability with respect to its approximate 
system, uniqueness and existence of harmonic solutions. 

§ 1, Consider the differential systems 

(i*i) / - p ( y,t ) + g ( y , 0 XO =y a Qo > o) 

(1-2) x' =p (x,t) x(t 0 ) = x 0 

where y, x, p, and q are n-dimensional vectors and where p and q are continuous 
functions on the product space A - I >* Let 1 be the interval o < t < co ; and 
R n be the n-dimensional Euclidean space. Also let R + = [o, + co ]. Let || y || 
denote the norm of the element y. 

The behaviour of the system of differential equations of the form (IT) is fre- 
quently studied by comparing with those of the system (1*2) particularly when the 
perturbation term in (IT) is small, in some cases for sufficiently large t. One of 
the problems of principle interest that arise in this connection is concerned with 
the question of whether the boundedness, stability and other properties of the 
system (IT) are shared with. the system (1*2). The present note takes up this 
question when the perturbations satisfy certain conditions. 

The following definitions are required before proceeding further ■ 

Let y(t) and x{t) be any two solutions of the systems (IT) and (T2) respec- 
tively 

(di) The system (1- 1) or (1-2) is said to be equi-norm-bounded with respect 
to the system (1-2) or (IT), if for each a greater than zero. and t 0 > o, there exists 
a positive function /3 ( t 0> < x) continuous in t 0 for each a satisfying 

|| y(t) - x(t) I] < (S ( t 0 ,a ) 

for all || y 0 - x 0 || < a and t 0 . 

(df) The system (IT) or (T2) is said to be equi-stable with respect to the 
system (1*2) or (IT), if for each e > o' and t 0 > o, there exists a positive function 
rj (t 0 ,i) continuous in t 0 for each 5, satisfying 

liXO - *(0 I! < € 

for all || y 0 — x 0 || sg V (< 0 ,€) and t > t 0 . 
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tiote: If 0 (/ 0 ,«) = jS («), n (t 0 , 0 = (€) i. «. the? are independent of f 0 , equi- 

will be replaced by uniformerly— in the above definitions. 

§ 2. We require the following theorem proved elsewhere [5]. We shall state 
in a suitable form whose proof needs a very little modification of the proof given 
in [5]. 

Theorem 2. —Let the functions m{t) and WO be vectors, -continuous on f 0 ^ t <co 
and satisfy the inequality 

o+ Sup ^ i- ^ £ m(t + h) - m(t) ] A “ b »»(0 + H WO II ) 

where the function <y is continuous and defined on 1 X R + Then 

(2-1) m{t) < M(<) for t a t < co . 

where M(i) is the maximal solution of 

(2*2) z' - «> (f, £(0 + WO ) M(f 0 ) = M 0 

Then we have the following theorem 

Theorem 2.— Let the function «>(<, 2 + #) be continous, non-decreasing in 
■2-fx and defined on 1 X R + ; Suppose 

(i) || y - x + h [p ( ji, O - p (x, /)] 11 < II y - * II 

and 

(ii) q satisfy the condition 

(2*3) II ? (A 0 II < “( il JV II , 0 

Let M(*) be the maximal solution of 

(2*4) z' - (t, z + x) M (i a ) = M 0 

Then, if y(t) and x{t) are any two solutions of (1*1) and (1*2), we have 

MO - x(t) || < M(<) 

Proof: — Let m(t) «= MO — x(l) || 

therefore m(t + h) - U_>(/ -f- h) - x(t + h) II 

we have 

m{t + h) ^\)y - x h [p(y,t) - p{x,t )) 

+ h g(y,t ) II + || h || + II € a h II 
=e II y - * + h [p ( y , t)-p (x, 0] ii + A ii q( y> 0 II 
HMI ®i A || + ll h ii 

Since and tend to zero as h tends to zero, and using the conditions (i) 
and ii), we obtain 

ft“o+ Sup <!-)[ m[t 4 - h) - m(t) J <o (0 II y II ) 
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Since || y II - jj x II ^ || y - x || and by non-decreasing character of <■> in z + x, 
we have 


So+ Su p (t ) [ m(t + A ) “ ]<“<*• m (0 + II *M II ) 

The application of the theorem 1, yields 
(2*5) ||j(0 - <0H ^ M(0 

Remark Suppose x(t) = o, theorem 2 includes theorems 1, 2, 4 in [3J if 
the solutions of (2*4) are assumed to be bounded as . Also the global 
existence of the solutions of (1 # 1) followed immediately by use of the results of 
Wintner [10] if we assume that all the solutions of (2*4) exist for all t ♦ 

Now we are in a position to prove the following : 

Theorem 3: — Suppose that the assumptions of theorem 2 hold. Suppose further 
that (a) the differential equation (2*4) is equi-bounded, then the systems (1*1) and 
(1*2) satisfy the definition (d t ). 

(6) the identically zero solution of (2*4) is equi-stable, then the systems (1*1) 
and (1*2) satisfy the definition (d a ). 

The proof of this theorem is an immediate consequence of theorem 2 and the 
assumptions in view of the definitions. 

Note : The equ ; -will be replaced by uniformly— in the above theorem if the 
functions involved are independent of t 0% 

Remarks : — If p{x,i) =s o, we get the same concepts relative to an arbitrary 
solutions of y =* q(y t i), 

§ 3. In this section we shall consider the existence of harmonic solution of 
(1*1). Suppose that the functions p , q and « are smooth enough to ensure the 
uniqueness of solutions. Let p(y>t\ q{y,t) and w (t $ z + x) be periodic ini with a 
period unity. We are interested in the solutions of (1*1) whose period is also unity 
are called harmonic solutions. We have the following ; 

Theorem 4 Let the function (t,z + x) be continuous, nomdecreasing in 
z ^ x and defined in | x R + . Suppose 

(i) Wy-x + h [p (j, fy—pix, *)] IT C — *|| 

(ii) q satisfy the condition 

u<i{y> 0 IK w ( U ii > 0 

Suppose that the differential equation (2*4) has a harmonic solution. Then, if 
the system (1*2) has a bounded, non-decreasing solution, the system (1*1) has a 
harmonic solution. 

We require the following lemma : 

Lemma: Let E be an inductively ordered set and let F be a transformation 
from E into E such that for any X belonging to E we have FX > X. Then there is 
atieast one point belonging to E such that F (X) = X. 

t 231 ] 



This lemma is due to Zorn, cf [1] 

Proof of theorem I:— The proof of this theorem is similar to the proof of the theo- 
rem 4 given in [4]. 

Remarks : If, further A#, 0 is periodic in t with a period unity, the existence 
of a harmonic solution of (1-2) is assured by the result of [7| and it is also sufficient 
to prove the theorem. 

§ 4, In this section, we shall consider more general uniqueness theorem. 
Let us consider the systems (1*1) and (1*2) again. Suppose t 0 = y Q — x 0 o and 
restrict the interval of t to o ^ t < a. 

Theorems 1 — Let the function u> (t, z 4* #) be continuous, non-dccreasing in 
Z + x and defined for o < t < a, z + x > o. Suppose 

(i) (| v-# + h [p (y, t)-p (x> 01 II < 11.?-* H 

(ii) q satify the condition 

Suppose the only solution M(*) of z! = « {t } Z + x) on o. < t < a such that 

M(o 4 ) = M\o+J « o , 

iS: the tri-vial solution. Then there is at most one solution and the same solution 
for both the system (1*1) and (1*2) on o ^ t < a. 

It is clear that the uniqueness result in [6] can also be defined. accordingly. 

Note ; Similar theorems contained in [4] for two different systems. But we 
have considered here the original system and its approximate system. Moreover 
theorem 1 is a more generalisation of a lemma used in 14] as pointed out 

4“ 15]. . 
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STABILITY CONSTANTS OF PALLADIUM, COPPER, NICKEL, ZINC, 
AND MANGANESE COMPLEXES OF 3-HYD ROXY-3-PHENYL~l-p- 
ACETAMIDO PHENYLTRIAZENE 

By 

D, N. PUROHIT, S. M. DUGAR and N, C. SOGANI 
Chemistry Department , University of Rajasthan, Jaipur , 

[ Received on 2nd December, 1963 ] 

ABSTRACT 

Stabilities of palladium, copper, nickel, zinc, and manganese chelates of 3“hydroxy«3-phenyl-l-p- 
acetami dophenyl triazene have been determined in 70% v/v dioxan-water mixture by employing 
Bjerrum-Calvin pH titration technique. The titiation medium was maintained at a constant ionic 
strength (O' 1M, KC1) and at a temperature, 25 -J- 0*5°. Log K, palladium 23*204, copper 22 485, 
nickel 17*883, zinc 15*736, and manganese 12‘673. 

In recent years Sogani et aL l ~ s have reported hydroxy-triazenes as a very 
useful class of chelating reagents. They possess common functional grouping - 
N (OH) - N == N - . The parent compound of the series, 3-hydroxy- 1, 3-diphenyl- 
triazene, is a highly selective reagent for palladium 1 and copper 2 . Substitution of 
different groups at different positions of the aryl nuclei of the parent compound 
influenced the stabilities of resulting metal chelates. The present investigation 
deals with the determination of stabilities of bivalent metal chelates of para aceta- 
mido substituted hydroxytriazene, namely, 3-hydroxy-3-phenyl-l-p acetamido- 
phenyltriazene, in 70% v/v dioxan- water mixture, by employing Galvin and 
Melchior's method 6 commonly referred as Bjerrum-Calvin pH titration technique 6 . 
Titrations were carried out in duplicate using 40: 1 ratio of ligand to metal ion 
concentration; 


EXPERIMENTAL 

Ligand Solution . — 3-Hydroxy- 3-phenyl- 1-p-acetamidophenyltriazene was pre- 
pared by employing the method given by Sogani et aft Weighed amount of this 
compound was dissolved in 70% v/v dioxan- water mixture and diluted to give 
0*04M solution. 

Standard palladium , copper , nickel , zinc } and manganese solutions . — ‘Analytical grades 
of palladium chloride, copper chloride, nickel sulphate, zinc oxide 5 and manganese 
sulphate were used for preparing standard solutions. After standardising these 
solutions by usual methods, these were diluted to give 0*002M solution. 

Sodium hydroxide.— Approximately OTN carbonate free sodium hydroxide was 
prepared and standardised with analytical grade potassium hydrogen phthalate. 
It was diluted to give 0*02M solution. 

Dioxan . — The B. D. H. dioxan was purified by Weissberger’s method 7 . 

Potenliometfic titrations.— The pH- titrations of the ligands with standard alkali 
solution in absence and in the presence of different metal ions were carried out by 
using Cambridge bench-type pH-meter which gave values accurate to 0*01 pH unit. 
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The pH-meter was standardised before, and checked after, each titration with buffer 
solutions of pH 4*0 and 9*2 Electrode system consisted of glass electrode (pH range 
1 - 13) and reference as saturated calomel-electrode. 

10 ml. of 0*04M ligand solution in 70% v/v dioxan- water mixture were 
pipetted in a titration vessel. Sufficient amount of 0*02M nitric acid was added 
to it to lower the pH to about 2. Final ionic concentration was maintained at 
0*1M by adding 2*5 ml. of 2M KG1. When titrating in presence of metal ions, 
5 ml. of 2 X 10~ 3 M metal solution were added at this stage. The total valume of 
the contents was made to 50 ml. by adding varying amounts of dioxan and water in 
such proportions that it finally became 70% v/v dioxan-water mixture. It was titrated 
against 0*02M sodium hydroxide. The temperature of the solution throughout the 
titration was maintained at 25 ± 0*5°. The pH-meter readings obtained at a 
given alkali addition in duplicate titrations were reproducible with a maxi- 
mum variation of i 0*01 pH unit. The titration curves are given, in Figs. 
1 and 2. 



Calculations . — The horizontal distance between the reference titration curve 
and the curve obtained in presence of a metal ion gives the amount of metal bound 
ligand. This value divided by the total metal ion concentration gives the value 
for n . Ia this way the values of V at different pH values were calculated, 
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At any pH, the value of free ligand ion concentration, [T], was calculated 
from the total ligand concentration, [HT], and its dissociation quotient. This is 
based on the assumption that the amount of chelating agent over metal ion is so 
great that the removal of HT by chelation does not cause any significant change 
in the equilibrium : H T ^ + T“. The pK. value of 3-hydroxy-3-phenyl- 

1-p-acetamidophenyltriazene has been determined spectrophotometricaliy 8 as 
11*663. 



In calculating the values for and [T|, the concentrations were corrected for 
changes in volume, produced by addition of alkali during titrations. In this way, 
a service of values for n and [T], corresponding to different pH values, was obtained. 
The results are given in Table Nos. 1 to 3. The formation curves for metal chelates 
of 3-hydroxy-3-phenyH-p-acetamido-phenyltriazene are shown in Fig. 3. The 



values for log k 1} log k 2 , and log K a v log K) for different metal chelates were read 
directly from the formation curves and are given in Table 4* 
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TABLE 1 


Formation curve data for metal chelates of 3-hydroxy-3-phenyl-l-p-acetami- 

dophenyltriazene 



Palladium 


Copper 

n • 


pH 

: n • 

P IT] 

pH . 

Pi TJ 

1-90 

0-40 

11-866 

2*42 

0-40 

11-349 

1-95 

0-50 

11-818 

2-45 

0'65 

11-321 

2-00 

0-75 

11-770 

2-48 

0-90 

11-293 

2-15 

1-00 

11-627 

2-60 

•HO 

11-181 

2-25 

1-15 

11-531 

2-75 

1-60 

11-036 

2-35 

T40 

11-435 

2-90 

1-70 

10-888 

2-40 

1-50 

11-386 

3-00 

1-90 

10-789 

2-50 

1-65 

11-289 

3-10 

1-95 

10-690 

265 

1-85 

11-146 




2-85 

2-00 

10-947 





TABLE 2 

Formation curve data for metal chelates of 3-hydroxy - 3 - phenyl- 1-p- ace tarn 1- 

do phenyltriazene 



Nickel 



Zinc 


pH 

: n J 

P IT] 

pH : 

n • 

PI T] 

2-90 

0-20 

10-881 

4-70 

0-11 

9-088 

3-25 

0-32 

10-535 

5-12 

0-36 

8-669 

3-85 

0-50 

9-938 

5-30 

0-50 

8-490 

4* 15 

0-70 

9-639 

5-50 

0-72 

8-291 

4-90 

0-95 

8-891 

5-75 

1*14 

8-043 

5-50 

1-25 

8-293 

6-30 

1-42 

7-495 

5-75 

1-35 

8-044 

6-55 

1-50 

7-246 

5-85 

1-50 

7-945 

6*75 

1-62 

7-047 

6-15 

1-68 

7-646 

7-00 

1-80 

6-798 

6-30 

1-92 

7-497 

7-15 

1-92 

6-649 




7-25 

2-00 

6-549 
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TABLE 3 

Formation curve data for metal chelates of 3-hydroxy 3-phenyH -p-acetami- 

do phenyltriazene 




Manganese 




pH 

* n • 

pm 



6-31 

0-10 

7-480 



6-60 

0-20 

7-191 



7*00 

0-50 

6-793 



7-15 

0-61 

6*644 



7-30 

0-75 

6-495 



7-50 

0-93 

6-296 



7*80 

1-30 

5-998 



7-95 

1-48 

5-850 



8*00 

1-62 

5-801 



8-15 

1-88 

5-652 



8-25 

1-93 

5-553 



8-65 

2-00 

5-184 




TABLE 4 



Stability constants of metal chelates of 3-hydroxy-3-phenyl-l-p- 

phenyltriazene 

acetamido- 

Metal ion 

: log k x 

: log k a : 

lug &av 

: log h k% 

Palladium 

11*818 

11-386 

11-531 

23-204 

Copper 

11-335 

11-150 

11-225 

22-485 

Nickel 

9 938 

7-S45 

8-770 

17-883 

Zinc 

8*490 

7-246 

8-120 

15-736 

Manganese 

6-793 

5-880 

6-200 

12-673 
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DISCUSSION 


The consumption of alkali during the course of titration can be due to the 
ligand, the hydrolysis of metal ion, and the ligand protons liberated in complex 
formation. 

The ligand under study is a very weak acid (pK = 11*663). Thus, the ligand 
proton as such is not in a titrable form. Hydrolyses of these metal ions has been 
studied in 70% v/v dioxan-water mixture 9 . The pH for V = T5 in all the cases is 
lower than the pH at which the hydrolysis of these metar ions start. Hence, 
hydrolyses of these metal ions do not interfere in the stability measurements. More- 
over, there was no precipitation during chelation titrations, clearly ruling out the 
possibility of hydrolyses of these metal ions in the presence of large excess of the 
ligand. 

Thus, the consumption of an excess of alkali in chelation titration over the 
simple ligand titration is due to ligand protons liberated during the complex for- 
mation. This may be represented as: 

+ HT 4- OH- ^ MT+ + H a O (1) 

MT+ -b HT 4- OH- ^ MT 2 + H 2 0 (2) 

ii J he ° rder of stabil . it y ot hydroxytriazene metal chelate* has been found as 
palladium > copper > nickel > zinc > manganese. 

r ... T be . autbors express their gratitude to Prof. R. G. Mehrotra for providing 
facilities in the department. ° 
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CHEMICAL EXAMINATION OF SESBANIA GRANDIFLORA (LINN.) PERS. 
LEAVES : PART II. ISOLATION AND STUDY OF A SAPONIN 


By 

R. D. TIWARI and R. K. BAJPAI* 

Department of Chemistry, University of Allahabad, Allahabad 
[ Received .on 7th December, 1963] 

ABSTRACT 

The ethanolic extract of the leaves of Sesbania grandiflora (Linn), Pers. yields a sapoaia m. p. 
235‘240°C which on acid hydrolysis furnished an acid sapogenin, m. p. 300-302 a and galactose, rhamnose 
and glucuronic acid. The sapogeninhas been identified as oleanolic acid (3/3 hydroxy - A 12 oleanene* 
28 oic acid). 

The isolation and chemical examination of an aliphatic alcohol (grandiflorol) 
from the leaves of Sesbania grandiflora (Linn.) Pers., has been described in an 
earlier communication 1 . 

The ethanolic extract of the dried leaves of Sesbania grandiflora after the 
removal of grandiflorol 1 , on concentration and subsequent treatment with an excess 
of chloroform, gave a precipitate of crude saponin which was extracted with 
ether, benzene and acetone successively and dissolved in methanol and precipitated 
by addition of ether, thus affording a light yellow coloured, amorphous and 
hygroscopic compound, sintering at 140° and melting between 180-200°. This 
was crystallised several times from ethanol when a cream coloured amorphous 
hygroscopic compound, melting at 225-230° was obtained which gave all the 
characteristic tests for saponin 3 . 

Further purification was achieved by acetylation followed by deacetylation 3 
and deionisation 4 using Amberlite IR 120 (H) resin. After recrystallisation from 
ethanol a cream coloured amorphous powder melting at 235-240° was obtained. 
Although still hygroscopic, it was pure as only one spot was observed during its 
paper chromatography and paper electrophoresis. 

The saponin was hydrolysed using 4N sulphuric acid whereby an acid 
sapogenin separated out. The liquid part of the hydrolysate after the removal of 
sulphuric acid with barium carbonate and passing over IRA-400 ion exchange 


* Present Address : Government Degree College, Panna (M. P.) 



resin was concentrated and on paper chromatography it was found to contain 
galactose, rhamnose and glucuronic acid. This was also confirmed by paper 
electrophoresis using borate buffer. 

The acid sapogenin after crystallisation from ethanol was treated with 10% 
sodium hydroxide to give a sodium salt which on subsequent decomposition with 
hydrochloric acid gave the product which on crystallisation from methanol was 
obtained in the form of fine needles, m. p. 300-302°. 

On analysis the sapogenin has been found to be G (J0 H 48 O a . It gave yellow 
colour with tetranitromethane and also colour reactions of triterpenes 5 . The 
acidic character of the compound and formation of monoacetate showed that it is 
a hydroxy triterpenoidal acid. The preparation and study of its acetate, methyl 
ester and the methyl ester of the acetate and also its infra-red spectrum have 
confirmed its identity with oleanoiic acid (3 /3-hydroxy -A la -oleanene -28-oic 
acid). 


The leguminosae family is a rich source of triterpenoidal sapogenins^ 7 ’ 8 . 
The plant Sesbania grandiflora (Linn.) Pers« from the leaves of which oleanoiic 
acid has been isolated, belongs to the Papilionaceae sub-group of the Leguminosae 
family. From the two more species of the same genus i. e. from Sesbania aculeata 8 
and Sesbania aegyptioa 9 seeds, oleanoiic acid has been isolated as a sapogenin by 
the hydrolysis of the saponin occuring in them. 


Extraction of the Saponin : 


EXPERIMENTAL 


Ethanolic extract of the dry leaves of Sesbania grandiflora after removal of 
grandiflorol was further concentrated to a syrupy consistency. To this a very 
large excess of chloroform ^ was gradually added while stirring the solution well 
with the help of a magnetic stirrer; a brown gummy mass floating on the surface 
of the solution thus obtained was washed well with cold chloroform and then 
exhaustively extracted with ether and benzene in succession. The light grey 
coloured residue left now was dissolved in a minimum amount of warm methanol 
and ether was added to it. The saponin was obtained as a fiocculent yellow 
coloured powder. The process of dissolving in methanol and precipitating with 
ether was^ repeated several times, when a light yellow coloured, amorphous, 
hygroscopic compound sintering at 140° and melting at 180-203° was obtained. 
It was recrystallised several times from absolute ethanol and obtained in the form 
of a cream yellow coloured compound m. p. 225-230°. It was still hygroscopic 
and retained moisture tenaciously. It gave all the tests for saponin 2 . 


Purification of Saponin by acetylation and de acetylation : 

°f J} e ^ a P°^ n: The saponin was warmed with dry pyridine and acetic 
and Vft flw 18 a * dd f ^ to , lt# The contents were shaken well, the flask was stoppered 
The waw ; ^° Ur k° urs at room temperature and then added to crushed ice. 

well unfh * SO U ^- e a ? et deposited which was filtered at the pump, washed 

with eth^r , dlsso * ve d m minimum amount of chloroform and precipitated 
160 65°. * G Sapoam acetate was cream coloured amorphous powder m. p. 
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Deacetylation of the saponin acetate i 

The saponin acetate was treated with 2% aqueous solution of sodium 
hydroxide, the contents were shaken well for two hours and then warmed for few 
minutes on a water bath. The mixture was allowed to stand at room tempera- 
ture for twenty hours, the saponin acetate was gradually deacetylated to saponin 
and a clear solution was obtained. The saponin solution was passed through a 
column of Amberlite I-R 120 (H) ion exchange resin, after which it was 
concentrated and dissolved in ethanol. Acetone was gradually added to it while 
stirring the solution well; the cream coloured saponin so obtained was crystallised 
from ethanol. The saponin was still hygroscopic m. p. 235-240°, yield = 0T% 
of the airdried leaves. 

Paper chromatography of the saponin : 

„ The aqueous solution of the saponin was applied on a strip of Whatman No. 1 
filter paper. The chromatogram was developed for four hours in w-butanolracetic 
acid-water (4:1:5) using ascending technique 10 . The chromatogram was dried in air 
and sprayed with 25% trichloroacetic acid in ether and dried at 120° for 5 minutes 
when a single spot was observed. 

Paper electrophoresis of the saponin : 

The aqueous solution of the saponin was applied to a strip of whatman No. 3 
filter paper and was subjected to horizontal electrophoresis at 300 volts for six 
hours, using borate buffer (pH 9.4). The paper was dried at room temperature 
and then sprayed with 25% tri-chlor-acetic acid in ether and dried for five minutes 
at 120° when a single spot was revealed. 

Hydrolysis of the saponin : 

The saponin was dissolved in water and to this excess of 4N sulphuric acid 
was added; the mixture was heated on water bath for twenty hours. Frothing was 
reduced by occasionally adding a little amount of ethanol to the reaction mixture. 
The genin obtained as a white precipitate was separated from the hydrolysate by 
filtration under suction. 

Identification of Sugars in the saponin hydrolysate : 

The saponin hydrolysate after separation of the genin was shaken with 
barium carbonate, the precipitated barium sulphate and excess barium carbonate 
were filtered off and washed with hot water. The combined filtrate and the 
washings were concentrated under reduced pressure. The hydrolysate was also 
neutralised by allowing it to pass through Amberlite IRA-400 ion exchange 
resin and the neutralised solution was concentrated under vacuum. 

Paper chromatography of Sugar s\ 

The saponin hydrolysate was chromatographed on whatman No. 1 filterpaper, 
in two different solvents using (a) n-butanol-acetic acid-water (4:1:5) and 
(b) fz-butanol-ethanol-ammonia- water (45:5:1:49) using descending technique. After 
development, the chromatograms were sprayed with aniline-hydrogenphthalate. 
Three spots were shown by each chromatogram. The Rf values of these revealed 
the presence of galactose, rhamnose and glucuronic acid. 
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This was also confirmed, by running a chromatogram of saponin, hydrolysate 
spot along side with authentic galactose, rhamnose and glucuronic acid, the 
sugar spots of the saponin hydrolysate compared well with those of authentic 
sugars. Since glucose and galactose show nearly the same rate of migration 
there was little doubt regarding the identity of galactose. The benzylamine 
derivatives of glucose and galactose exhibit a marked difference in 
their movement 12 . A chromatogram of benzylamine derivative of saponin 
hydrolysate alongside with benzyl amine derivatives of authentic glucose 
and galactose was developed in w-butanol- ethanol- ammonia- water (45:5:1:49), 
Ihe chromatogram was sprayed with 0*01% solution of ninhydrin in butanol. 
The saponin hydrolysate showed a spot corresponding to the benzyl-amine 
derivative of galactose and not of glucose. Absence of ketoses was confirmed., by 
spraying one of the chromatograms of saponin hydrolysate, with naphthoresorcinol- 
trichloro acetic acid when no spot was revealed. 


Paper electrophoresis of saponin hydrolysate : 

A spot of the saponin hydrolysate was applied on a strip of whatman No. 3 
filter paper, along side with this, spots of authentic sugars, galactose, rhamnose 
and glucuronic acid were also applied. The paper was subjected to horizontal 
electrophoresis at 300 volts for six hours using borate buffer. The paper was 
then sprayed with aniline-hydrogen-phthalate reagent and the presence of 
galactose, rhamnose and glucuronic acid was confirmed. 


Isolation of the genin : 

The genin was crystallised from ethanol when it gave fine needles m. p. 
290-292°. The genin was taken in ethanol and treated with 10% aqueous sodium 
hydroxide, the mixture warmed on the water bath and kept for few minutes at 
room temperature. Alcohol was distilled off under reduced pressure. The 
deposited crystalline sodium salt was extracted with ether, but on distilling off 
ether no residue could be obtained. The sodium salt left after extraction with 
ether was dissolved in ethanol and acidified with an excess of dilute hydrochloric 
acid. The precipitated genin was washed well with water and crystallised from 
ethanol giving needles melting at 300-302°. 

It gave yellow colour with tetranitromethane, a red violet colour with 
acetic anhydride and concentrated sulphuric acid (Liebermann-Burchard reaction), 
a purple colour with commercial thionyl chloride (Noller’s reaction) and purple 
colour with chlorosulfonic acid (Brieskorn test) 5 . 

Found : G, 79’09% j* H, 10*53% ; Mol, wt. 458 (neutralisation) ; calculated 
for CsoH^Oa, G, 78*94%, H, 10*52%, Mol. wt. 456. 

30 

|>] D = + 78° (c, 1*00 in chloroform) ; 


Infra-red spectra : 

The peaks at 3420 an - 1 (~OH), 2933 enr 1 (- CH 8l -GH, ), 1701 crrT 1 
(-COOH), 1464 errr 1 (C-GH 3 ) and 1389 1366 1347 cmr 1 , 1325 cm~\ 
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1305 mr\ 1269 cmr 1 (oleanolic acid type) 6 were obtained in the I-R spectrum of the 
compound in KC1. 


Acetylation of the genin', 

The acid genin (800 mg) was mixed with acetic anhydride (4ml) and freshly 
dried pyridine (2 ml). The mixture after being heated for three hours on a water 
bath was allowed to stand overnight at room temperature. The genin acetate 
was crystallised from methanol-chloroform as feathery white needles m. p. 260-262°. 

Q 77*31% H, 10*19%, acetyl group, 8*22%; calculated for C 30 H 47 ~ 
GOGHg G, 77*11%, H, 10*04%, acetyl group 8*64% (lor one-acetyl group. 


At 


Hydrolysis o f the genin acetate . 

The acid genin acetate (200 mg) was refluxed with 10% methanolic 
potassium hydroxide, for three hours, methanol was distilled off and contents 
treated with water, the sparingly soluble salt was filtered, dissolved in little 
amount of ethanol and ethanolic solution was neutralised with dilute hydrochloric 
acid. The precipitate obtained was crystallised from methanol to fine needles 
m ; P» 800-302°. No depression in the melting point was observed when mixed 
with the original acid genin. 


Methyl ester of the acid genin: 

The genin (200 mg) in ether-methanol solution was treated with an excess 
of the ethereal solution of diazomethane. The solvent was evaporated off and the 
residue was crystallised from dilute ethanol yielding fine silky needles m. p. 
195-97°, When treated with stannic chloride in thionyl chloride it gave a light 
pink colour, becoming deep red on standing. Found, G, 79*02%: H, 10*57%; 
calculated for C 31 H 30 O a ; G, 79*14%, H, 10*64%. 


Acetyl methyl ester of the genin \ 

Acetyl methyl ester of the genin was obtained by the diazomethane treatment 
of the acid genin acetate in the usual way. After crystallisation from chloroform- 
methanol it was found to melt at 216-219°. It developed a light yellow colour 
with stannic chloride in thionyl chloride. 

Found; G, 77*59% H, 10*27%; calculated for C 33 H 5a 0 4 ; G, 77*34%; H, 10*16% 

The melting point of the acid genin was not depressed when mixed with 
an authentic sample of oleanolic acid. 


ACKNOWLEDGEMENTS 

Authors are thankful to Dr. I. P. Varshney, Chemistry Department, Muslim 
University, Aligarh for a sample of oleanolic acid for comparison and to the Kanta 
Prasad Research Endowment Trust for a research scholarship to one of them 
(R. K. B.). 


[ 243 ] 



REFERENCES 


1, Tiwari, R. D., Bajpai, R. K, and Khanna, S. S, 5 Arch , der Pharm*, 297, 310 (1964). 

2, Sannie, G., H., Anal . Biochim . Med„ 9, 175, (1948)* 

3- Farmer, R. H. and Laidlaw, R. A., J. Chans Soc 4201 (1955). 

4. Varshney, I. P. and Khan, Mobd. S. Y., J . Sci, indusir . Res. India, 21B, 30, (1962). 

5. Paech, K. and Tracey, M. V., Modern Methoden Der Pflanzenanalysc, Springer Verlag, 
Berlin, Volume 3, p. 64 (1955). 

6. Shatzke, G , Lampeit F. and Tschesche, R., Tetrahedron , 18, 1417 (1962). 

7. Varshney, I.P. and Khan, Mohd. S. Y.i J. Pham , sci., 50, 923 (1961). 

8. Varshney, X. P, and Khan, Mohd. S. Y., J. Sci indusir . Res., India, 21B, 401 (1962), 

9. Farooq, M. O., Varshney, I. P. and Khan, Mohd. Amer, Pham. Assn . Set . Edn., 

48, 466 (1955). 

10. Sarkar, B. and Rastogi, R P. s J . Sci. indusir. Res. India , 19B, 106 (1960). 

11. Jensen, K. B. and Tennoe, K., J. Pharm . Pharmacol 7, 334 (1955). 

12. Barker,* S. A., Bourne, E, J. s Grant, P. N. and Stacey, M., Nature , 177, 1125 (1956). 


[ 244 ] 



FINITE STRAIN IN A ROTATING SHAFT 


By 

Md. MAHFOOZ ALI SIDDIQUI 
Dept . of Mathematics , Osmania University y Hyderabad 
[Received on 13th November, 1963J 


ABSTRACT 

Using non-liner strain components : 


2s.. = 


= u. .+ 


u • 
J > * 


HI « Ul • 

k,j 


and the non-liner stress strain relations : 


T y - ( 1—I i) S y + IJ e.j e a j 

and a, perturbation constant, K, the problem of an isotropic circular shaft rotating with angular 
velocity w, is solved completely, with the radial and longitudinal displacements calculated. 

We consider a solid shaft of circular cross-section of radius a, rotating with 
an angular velocity to. The shaft is taken to be of an isotropic material, and the 
symmetrical manner in which it is strained leads us to assume a radial displace- 
ment and a longitudinal extension only, given by 

a = K « (r) (1,1) 

W = K w (X) 

where u (r) is a function purely of r and w (-2) a function purely of Z , and K is a 
perturbation constant, whose squares and higher powers we propose to neglect 
in the linear theory, where as second order terms in K are retained when applying 
the non-linear theory. 

For the linear theory the components of strain are 

9 ft ==rvU ^ ; e r 6 = e 6X 5=8 $ Xr 



e zz = K •' (*) 


The components of stress are 

rr = \K ( u'+ ~ + W' ) + 2 K n u ' 

$d = \K ( U> +— + w) + 2K i* - 
\ r / r 


( 1 . 2 ) 
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(1.3) 


XX = XK^ u' -f — + W ) + 2 K P a' 
r0 = 9X = Xr = O 

These stresses must satisfy the body-stress equations 

( lr7 4- -L Pr -68) + pr* 1 = 0 ; ■- X = 0 
■0r r OX 

which give 


K(„.+ 4--i) + 

and W* - o 

Tbus« = Ar-P^ (x “_^ } r 3 


X - 2p 


rw* 


( 1 . 4 ) 


(1.5) 


W = B X 

where A and B are constants and the particular solution a = — • of (1.4) is not 

admissible, since « has to be bounded at r = 0, the shaft being complete up to 
the centre. The boundary conditions to be satisfied arc 


r r = 0 for r — a 


( 1 . 6 ) 


the outer surface being stress-free. 

Also, for the resultant longitudial tension to vanish over the plane ends it is 
necessary that 


/ 


r XX dr = 0 


(1.7) 


The conditions (1.6) and (1.7) give 


2 (X + P) A + X B - 


2 X + 3 P 
4 k (X + 2 P) 


p m a a a 


= 0 (1.8) 


2 ^ + <*+ 2 « B -Trm-° (1 - 9 > 

which immediately yield 

K _ p 0,2 Q ‘> (x + P) (X + 6 P) 

8 P (X + 2 P) (3 X + 2 P) 

( 246 ] 


( 1 . 10 ) 



B K= - X P (2 X+ll^H 
8 /a (X + 2/a) 


(Hi) 


Therefore the radial displacement as given by the linear theory is 
u = K u (r) 

_ P m 2 a 2 ( A. + /*) (X -j ~ 6 /A)r pm* 3 (1-12) 

tip- (\ + 2 V) (3A + W) 8 (vf 2 /a) r 


and the displacement in the 2-direction is 
W = K W (X) 


= - p a 2 a* x 


J2 X+_11/a)_ 

8/a (\ + 2f*) 


(1.13) 


When applying the non-linear theory of elasticity, j.c. when we take the 
non-linear strain components (1.10) and the non-linear stress-strain relations 
given by 

T - [X 1 ! (! “ !i)] S y + 2 [/a - (X + /*) Ii] Gjj ' 4/i e aj . s^j . 


we assume 

!i = 1^ + Kl> (f) (1.14) 

w - + Ks (x) 

where u lt and w t are the linear theory solutions, (1.12) and (1.13) 


£ 


ff 


The strain domponents corresponding to (1,14) will be 
a Kui + K a (o' - i U/*) 


= k-^ + k 

r \ r r* / 


« a = KW/ + Ka (5' - i Wx'») 

F r0 " S 0 X = £ Xr= 0 
The stresses will be 


(1.15) 


rr - X [ K ( + JjL + W/ ) + K» 

{<’ , +7+o-K“‘ ,a+J £ +w ‘'’)} ) 

/ \ 2 

- XK 2 ^ + _ 2JL + W x ' ) + 2/i IK < + K* (a' - | Bl '*)] 
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- 2 (\ 4* f*) k. a [ «' ( V + ~ - + w 'i ) ] 

- 4^ k 2 w x /2 (^*^) 

M =x[ K(tV + J^- + W 3 ') + K’ 

{( + -f +*') - *( V a + -£ + Wx«) )] 

- XK»(V+ J^_+ W/) 2 + 2^[ K JJ-+ K^( i - ^- ) ] 


-2(X + M)K 2 [^(V + ^+Wx')] 

- 4/* K a 

1* 


(1.17) 


a -x|^K |v + -^- + Wi'j + K . 2 

J^ + l. + ,' j ^' 2 4- ,^-+ WV» j j> 

- XK» (* + + Wi' j 2 + 2m Wi'. + Pfc' " t Wx' 2 ) 

-2(X + m)K^W/ |< + ^-+Wx'j| (i i8) 


- 4m K* Wj' 2 

The equations of equilibrium are 


3 , 1 

rr + — (rr 
r 


3r 


00) + Pr“ a = 0 



which reduce to 


v" + — ^ v — = - (lOxA + <*XB + 16ft a + 2f*B) 


pw 2 


K(x + 2/t) a 


+ ,„ 49 >,+ ? 3 g_- P V« r3 


16K a (X + 2ft) 3 


(1.19) 


and S" = 0 
Therefore 

v = Cr - (10XA + 4XB + 16f*A + 2ftB) 

, 49X + 83/i 


P 0)^ 


8 (X + 2f*) a K 


24 X 16K a (X + 2? ) E 


f 5 


and S = Dx 

where G and D are constants. 

If we write 


( 1 . 20 ) 


Aj = AK - p.» .» — <) + »> ^ + 6y ) , 
8ft (X + 2ft) (3X + 2ft) 

B _ nv _ x P“ a <J a (2X + Ilf*) 

B ‘ - BK “ ■ x 8, a + 2tr 

j IOxAj + 4XBj^ 4“ lS/^Aj + 2P , B 1 2 

1 8 (X + 2ft)a Hw 

M, 49 X -r _ggf t _ - 2u4 

1 24X 16 (X + 2 ft) 3 


( 1 . 21 ) 


the boundary conditions over the curved surface and the plane ends, given by 
a 

dr = 0 become 


rr 


= 0 for r = a and ) r XX 


/■ 


2 (X + /*) CK a + xDK a = 9 (x+ ft) A x a + -?£_ Bi 2 + 2 (3x + 2ft) AjBj 


pW4 fl2 ^5 (4X + 5ft) Aj + (7X + 3ft) B,| + 0<X + 23ft) 


4 (X + 2ft) ^ 

+ 2 (2x + 3ft) L] a* + 2 (3x + 5ft) M x a 4 


( 1 . 22 ) 


and 


2X CK 2 + (X + 2ft) DK a = 5xA a ® + 7/2 (x + 2ft) B x 2 + 4 (x + 2 ft) A^ 

(1.23) 


P“ 2fl2 _ (5XA, -2ftB x ) + 


4 (X + 2ft) 

4“ 2xBx a a 2XMj a 4 


6 ± 
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From (1.22) and (1.23) 

rn _ 4X 2 + 25XM + 18^ a A a X + f* B .a+. & ± A) (X + M A, B. 
CK VM3X TW~ 1 3T+W + A (3X + 2/*) 1 1 


+ 

+ 


pWa f I (1 5A X + 7Bi) X 2 - 3 (5A t +SB x )XM - 2(25A X + 3B X )M 2 ) 

+-2M.) (3a+ 2M) l 1 


8m (A+2M) (3a+2M) 

p* to 4 a 4 

128m (X+2M) 2 (3X+2/*)" 


|35x 2 + 15Sx/* + 138M a | 


(1.24) 


1 (X+P) (X+6m) 2 T 

T 2M (3x+2m) 

+ (X+2M) (3X+5M) a4 M 
~ 2 m(3x + 2m) 1 


and 


nKa _ 6\ a + 23XM + 12M a fi , _ 2x (4x + 5 m) a , 
2M (3x -}- 2 m) 1 M (3X -f" 2M.) 


2 (4m 2 + 5XM-X 2 ) 
■*" M (3\ + 2m) 


AjB, 


+ ®rir+Wx"+2s( 4 <l5A ‘ +7Bl) x ’ +2 (60A, + 19B ‘ ) «* +5 < 5A ‘+ 2B > a } 


pa to 4 a 4 x (35X+ 62 m) _ 5x 2 -H 5 AM + 6 p2 a% L 
64; (X+2m) 2 (3X+2m) ” m (3X + 2m) 1 


+ 


(3X + 5m) (5X + 2m) fl4 M 
(3X + 2 m) 


Therefore the radial displacement is 


u — 


( p ^ « a ( X + g) (x + 6 m) rK2 

\ 8m (X + 2m) (3x + 2M) + 


where the value of CK 2 is given by (1.24) 
The longitudinal extension is given by 

W = DK 2 - po> 2 o 2 M2 x+_Um) 

8m (X + M) J 

where the value of DK 2 is given by (1 .25) 


(1.25) 

~ (s(\ + 2M) + Ll j" 

, (49x + 83 /a) p \ ** 

^ 24 X 16 (X 4* 2/Aja 
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ABSTRACT 

An examination of the data presented shows that more nutritive grain can be produced in a short 
duration variety of paddy, if the plant is subjected to progressively depleted soil moisture conditions 
during the ripeniDg period, Piotein accumulation is greatly increased if plant is treated with potassium 
in two split doses. With regard to its composition (amino acid), it is not practically affected by any 
treatment. 

INTRODUCTION 

During the last three decades or so f considerable amount of work touching 
the various aspects of the nutritional quality of rice has been carried out in this 
country. Soil moisture plays the most dominating role in plant's life and is 
largely responsible for the accumulation of various nutrients in the grain. But 
unfortunately very little information is available which can throw light on this 
important subject. It is also a matter of considerable significance to find out 
the effect of some specific salt upon the accumulation of various nutrients in the 
grain when it is applied in different soil moisture conditions. 

The authors (1963) have found that in a short duration variety of paddy, 
stoppage of irrigation a few days earlier than harvesting and application of 
potassium in instalments cause greater production of grain per unit area. The 
chief objective of undertaking this study was to see the extent to which these 
treatments affect the nutritive value of the grain. So far as is known no such 
work has been reported which can throw light on this important aspect and the 
present study is an attempt to fill up the gap to some extent. 

EXPERIMENTAL METHOD 

Pre-harvest Operations . 

An early maturing variety of paddy viz., N a2 was selected for the purpose. 
Other details regarding filling of pots with soil, raising of seedlings in the 
nursery, transplanting, weeding, in ter- culturing etc., remained the same as in the 
previous paper (Authors, 1963). 

Plants growing in the pots received liberal but equal and measured amount of 
water till the commencement of the grain maturation period i. e, 9 upto 60 days 
after transplanting. At the beginning of the period, all the pots were brought 
inside the glass-house so as to avoid the interference due to rain, dew etc., 
during the course of the experimentation. 

Pots were now divided in two sets. The first one received liberal but equal 
and measured volume of water and the level 'of standing water was maintained 
2" above the soil surface till the end of the harvest, and was denoted as I n . iu 
case of the other, after keeping 2" of standing water at the commencement - of the 
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period, no further supply of irrigation was provided and the plants were then 
forced to grow under progressively depleted soil moisture condition. This set was 
denoted as l x . 

Plants of each set received potassium in the form of potassium sulphate at 
the rate of 40 lb K 2 per acre in the following manner: 

(1) K 0 — Control, i, e., no application of potassium. 

(2) K x — The whole amount was applied just. after one week of transplantation 

in a single dose. 

(3) K a — The whole lot was split in two equal doses, the first one was applied 

one week after transplanting and the other one month later. 

(4) K 3 — In this case the amount was divided in three equal doses and each 

one was applied at an interval of one month starting from one week 
after transplanting. 

Pots under the following 8 treatment combinations were arranged in a, v 
factorial experimental design with three replications: 

(1) I a K 0 , (2) I a K t (3) In K 2 , (4) I„ K„ (5) I x K 0 , (6) I x K l9 (7) I, K a 
and (8) I x K.. 

Harvesting and Preparation of Grain Sample: 

Harvesting of the crop was done when the grains were fully ripened. It was 
noted that moisture depleted plants took about 5 days more than the normally 
irrigated ones. After the harvest, grains were separated from the straw and 
their husks were removed very carefully. The hulled rice was then ground to ^ 
fine powder in an electrically operated micro-grinding mill and the sample 
was then dried in an electric oven, bottled and sealed. 

Chemical Analyses : 

Help of Loomis and Shull (1937) was taken for the determination of total 
carbohydrate and crude protein in the sample. The value of crude protein was 
obtained by multiplying the total nitrogen with the factor 6‘25. Soxhlet apparatus 
was used to obtain ether solubles. The period of extraction was 8 hours and 
the condensation rate was maintained at 4-5 drops per second. The extract was 
then dried for 30 minutes at 100°C, cooled and weighed. 

For the detection of various amino acids in the sample, paper chromato- 
graphic method as suggested by Block and Weiss (1956) was used. The chromato- 
gram developed from each of the hydrolysates was compared with the standard 
which was prepared under similar conditions. The help of Rp values as given by 

Brimley and Barrett (1956) was also taken to recognise them. 

EXPERIMENTAL FINDINGS 

Data of nutritive value of the rice as presented in the accompanying Table 
clearly show that stopping of irrigation in an early maturing variety a few days 
earlier than harvesting causes more accumulation of total carbohydrate, ether 
soluble and protein than providing continuous irrigation. Higher concentration 
carbohydrate in both the sets was also favoured when potassium was applied in a 
single dose, but this method of application was not found useful for the synthesis 
pf ether solubles and protein. 

L 252 | 



Table Showing the Effect of Soil Moisture Conditions and Potassium 
Application on the Average Chemical Composition 
(percentage of dry matter) of Rice 


Treatment 

Total 

carbohydrate 

Ether 

solubles 

Crude 

protein 

Amino acids detected 

(1) IjK, 

70-76 

3-71 

7-35 

Alanine (Al), aspartic 
acid (Asp), cystine (C), 

glycine (Gl), histidine (H), 
isoleucine (Isol), lysine (Ly), 
methionine (M), proline (Pr) 
and serine (Se). 

(2) I n K x 

79-00 

3-42 

6-04 

Al, arginine (A), Asp, G, Gl, 
H, Isol, Ly, M, Pr, and Se. 

(3) I n K a 

73-53 

3-43 

8-01 

Al, A, Asp, C, Gl, B, Isol, 
Ly, M, Pr, and Se. 

(4) I»K, 

74-26 

2-95 

7-31 

Same as treatment (2) . 

(5) IiK 0 

84-30 

3-73 

7-22 

Do. 

(6) IxK, 

83-33 

3-57 

7-70 

Do. 

(7) I X K 2 

81-52 

3-68 

8-44 

Do. 

(8) I x K a 

72-82 

3-32 

7-96 

Do. 


Further examination of the Table reveals that application of potassium in 
two instalments encouraged greatly to produce grains containing higher content 
of protein. It is also clear from the Table that nearly the same amino acids are 
detected in each treatment and it may be concluded that protein constituents 
are not affected either by moisture conditions or potassium application. 

DISCUSSION 

One of the remarkable findings of the study is that stoppage of fresh supply 
of irrigation water a few days earlier than harvesting of a short duration variety 
of paddy causes more accumulation of total carbohydrate, ether solubles and 
protein in the grain. As has been stated earlier, moisture depleted plants took 
about 5 days more in completing their life-cycle than the normally irrigated ones, 
and therefore, during the extra 5 days the former might have accumulated more 
nutrients than the latter. Low level of soil moisture generally causes better 
aeration and thereby more extension of root system in the direction of available 
moisture and nutrients. This extension of root system might have helped the 
plant in overcoming not only the severe drought injury but also in greater 
absorption of nutrients from untackled area, and thus the nutritive value of ^ the 
grain might have ultimately been increased. Miller and Duley (1925) obtained 
in maize, and also Kido and Yanatori (1961) in rice, more nutritive grains iu 
plants grown on well aerated soil. 

The grains of normally irrigated plants ripened earlier and this might be 
due to better supply of nitrogen to the plant (Russell, 1961). This increased 
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supply of nitrogen might have encouraged new vegetative growth, and the 
nutrients which should have gone for the enrichment of the grain, might have 
been diverted towards the new developing regions thus lowering the nutritive 
value. 

Application of potassium, except in few cases, failed to increase the content 
of either carbohydrate or ether soluble compounds, but its application, _ particularly 
in split doses, definitely encouraged more accumulation of protein in both the 
conditions of soil moisture. As during the later period of growth, potassium is 
secreted from the plant body to the soil (Sayre, 1948), it is likely that split 
application might have met the deficiency at any time and thus reduced the rate 
of protiolysis to amino acids and amides (Wall, 1940 and Richards and Berner, 
1954). This might account for the greater amount of protein in the grain. 

Paddy grains as affected by various treatments have practically the same 
amino acids and it appears that protein constituents are genetically controlled. 

SUMMARY AND CONCLUSIONS 

The present investigation was undertaken to see the extent to which soil 
moisture conditions during the grain maturation period and potassium application 
affect the nutritive value of the grain of a short duration variety of rice. The 
main findings and conclusions may be stated as below: 

(1) Progressively depleted soil moisture is beneficial for the production of 
more nutritive grain. 

(2) More proteinacious grain can be produced by the application of split doses 
of potassium. Application in two instalments have been found to be most useful. 

(3) Grains as affected by different treatments practically have the same 
amino acids and it appears that the protein constituents are genetically controlled. 
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